Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 303,121 Δ1 bp coding (131/615 nt) yagU → DUF1440 family inner membrane acid resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913303,1190C.98.8% 402.5 / ‑0.8 84coding (129/615 nt)yagUDUF1440 family inner membrane acid resistance protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (43/40);  total (44/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

TACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGACAGATCCCAATGCGGCTGTTTA  >  NC_000913/302975‑303265
                                                                                                                                                 |                                                                                                                                                  
tACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTgg                                                                                                                                             >  2:740672/1‑151 (MQ=255)
 aCAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTgg                                                                                                                                             <  1:1349125/150‑1 (MQ=255)
 aCAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTgg                                                                                                                                             <  2:1007160/150‑1 (MQ=255)
  cAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGat                                                                                                                                           <  2:69977/151‑1 (MQ=255)
  cAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGat                                                                                                                                           >  2:1030673/1‑151 (MQ=255)
    aGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGatat                                                                                                                                         <  2:973820/151‑1 (MQ=255)
      cAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGt                                                                                                                                       >  1:723548/1‑151 (MQ=255)
       aGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGtt                                                                                                                                      <  2:992321/151‑1 (MQ=255)
       aGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGtt                                                                                                                                      <  1:1120351/151‑1 (MQ=255)
       aGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGtt                                                                                                                                      <  2:1146788/151‑1 (MQ=255)
          gTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTaa                                                                                                                                   <  2:1192796/151‑1 (MQ=255)
            aGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATg                                                                                                                                 >  1:627299/1‑151 (MQ=255)
             gctatGAATATATGTGAAAAAACTCCACCGAATCTCATACGTTATGGTCTTGCTGTATTCATTGTGCTGATTGCTGGCGTTGTTTCCGCATTCGTGACGAGGGGGGCTGAAGTTTCATTGCCGCCACGTAGCCGGTGGATATGTTTAATg                                                                                                                                 <  1:330798/148‑1 (MQ=255)
              atatGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCa                                                                                                                               <  2:408124/151‑1 (MQ=255)
                  gAATATATTTGAACAAACTCCACCGAACCGCATTCGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCTTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCgt                                                                                                                           <  2:260770/151‑1 (MQ=255)
                  gAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCgt                                                                                                                           <  1:204688/151‑1 (MQ=255)
                  gAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGACTGAAGTTCCATTGCCGCCACGTCGCCGTTGGATATGTTTAATGCAGCgt                                                                                                                           <  1:382392/151‑1 (MQ=255)
                    atatatTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCgtg                                                                                                                          >  2:599704/1‑150 (MQ=255)
                      atatTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTgg                                                                                                                       >  2:296601/1‑151 (MQ=255)
                      atatTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTgg                                                                                                                       >  1:61628/1‑151 (MQ=255)
                         tttGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGccc                                                                                                                    >  2:658093/1‑151 (MQ=255)
                          ttGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGccc                                                                                                                    <  1:1018323/150‑1 (MQ=255)
                           tGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCgg                                                                                                                  >  2:1333394/1‑151 (MQ=255)
                           tGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCg                                                                                                                   <  2:1066886/150‑1 (MQ=255)
                           tGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCg                                                                                                                   <  2:866373/150‑1 (MQ=255)
                            gAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCgg                                                                                                                  <  1:434523/150‑1 (MQ=255)
                             aaCAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGACGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCg                                                                                                                   >  1:309418/1‑148 (MQ=255)
                             aaCAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCgg                                                                                                                  >  1:49448/1‑149 (MQ=255)
                             aaCAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGaa                                                                                                                <  2:303903/151‑1 (MQ=255)
                             aaCAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGaa                                                                                                                <  2:202537/151‑1 (MQ=255)
                              aCAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAAt                                                                                                               <  1:579068/151‑1 (MQ=255)
                                   cTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGCTTAATGCAGCGTGTGGCCCGGAATCAtt                                                                                                           >  2:73615/1‑150 (MQ=255)
                                       aCCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCa                                                                                                      <  2:1302277/151‑1 (MQ=255)
                                        ccGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCa                                                                                                      <  1:47228/150‑1 (MQ=255)
                                        ccGAACCGCAGACGGTCTGGTCTTGCTGCATTCATTGGGCTGATTGCGGGCGGTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATc                                                                                                       >  1:805628/1‑149 (MQ=255)
                                          gAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAggg                                                                                                   >  1:538239/1‑151 (MQ=255)
                                            aCCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCt                                                                                                 <  1:130849/151‑1 (MQ=255)
                                             ccGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTg                                                                                                >  2:224550/1‑151 (MQ=255)
                                              cGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGc                                                                                               <  2:861699/151‑1 (MQ=255)
                                               gCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGACGTGGGGGGCTGAAGTTCCATTGCCGCCCCGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAACCATTAATCAGGTCTGc                                                                                               >  1:829956/1‑150 (MQ=255)
                                                   aCGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGcc                                                                                          >  1:334376/1‑151 (MQ=255)
                                                        aTGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAat                                                                                      <  2:863025/150‑1 (MQ=255)
                                                         tGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAattg                                                                                    >  2:966596/1‑151 (MQ=255)
                                                          ggTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAattga                                                                                   <  1:102977/151‑1 (MQ=255)
                                                               tgctgcATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCt                                                                              >  2:1043284/1‑151 (MQ=255)
                                                                gctgcATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTc                                                                             <  2:683845/151‑1 (MQ=255)
                                                                gctgcATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTc                                                                             <  2:572501/151‑1 (MQ=255)
                                                                     attcattGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTcgc                                                                           <  1:738625/148‑1 (MQ=255)
                                                                     attcattGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTcgc                                                                           >  2:739250/1‑148 (MQ=255)
                                                                     attcattGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTa                                                                        <  1:168494/151‑1 (MQ=255)
                                                                         attGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAAttt                                                                    >  1:1389186/1‑151 (MQ=255)
                                                                         attGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAAttt                                                                    >  1:506535/1‑151 (MQ=255)
                                                                                tGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAAt                                                             >  2:166633/1‑151 (MQ=255)
                                                                                  tttGCTGGCGTTGGTTCAGCACTAGTGCAGTGGGGCGATGACGTTGCATTGCCGCCACGTAGCCCGTGGTTATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCGGCTGGCGAAATTGATTGCTCGCGTAATTTTCTCAAt                                                             >  1:485336/2‑149 (MQ=255)
                                                                                  aTTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATcc                                                           >  2:968443/1‑151 (MQ=255)
                                                                                    tGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCAc                                                         >  2:510608/1‑151 (MQ=255)
                                                                                     gCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCAcc                                                        >  2:423935/1‑151 (MQ=255)
                                                                                     gCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCAcc                                                        <  2:844759/151‑1 (MQ=255)
                                                                                      cTGCCGTTCTTTCCGCATTCGTGAGAGTGGTGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGTCTATGTTTAATTCAGCGTGTGGCCCTTAATCATTCATCAGGGCTGCAGGACAAATTGATTGCTCGCGTAATTTGCTAAATCCAc                                                         <  1:6094/149‑1 (MQ=255)
                                                                                         gCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGtat                                                    >  2:100833/1‑151 (MQ=255)
                                                                                           gttgttTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTatat                                                  >  1:58699/1‑151 (MQ=255)
                                                                                             tgttTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATAttt                                                >  2:236813/1‑151 (MQ=255)
                                                                                              gttTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATAttt                                                <  2:513040/150‑1 (MQ=255)
                                                                                                  ccGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTttctt                                           >  2:884857/1‑151 (MQ=255)
                                                                                                      aTTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCg                                         >  2:13415/1‑149 (MQ=255)
                                                                                                      aTTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCg                                         >  1:1080880/1‑149 (MQ=255)
                                                                                                       ttCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCg                                         >  2:1280311/1‑148 (MQ=255)
                                                                                                        tCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGAc                                     >  2:584849/1‑151 (MQ=255)
                                                                                                            gAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGt                                 <  2:178260/151‑1 (MQ=255)
                                                                                                                gTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTgg                              <  1:752567/151‑1 (MQ=255)
                                                                                                                       gCTGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGa                        >  1:370581/1‑150 (MQ=255)
                                                                                                                         tGAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGAc                       >  1:760775/1‑149 (MQ=255)
                                                                                                                          gAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGAc                       >  2:40773/1‑148 (MQ=255)
                                                                                                                          gAAGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGACa                      <  1:42283/149‑1 (MQ=255)
                                                                                                                           aaGTTCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGACAGAt                   <  2:873048/151‑1 (MQ=255)
                                                                                                                           aaGATCCATTTCCGCCACGCAGCCGGTGGATATGTTTAATGCAACGTGTGGCCCGGAATCATTAATCACGGCTGTAGGCCAAACTGATTGCTCGCGTAATTTTCTCATTACACCGTATATTTTTCTTCGGGACTGGTCGGCGCTg                         >  1:821090/1‑145 (MQ=255)
                                                                                                                               tCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGACAGATccc                >  1:1093582/1‑150 (MQ=255)
                                                                                                                               tCCATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGACAGATccc                >  1:436788/1‑150 (MQ=255)
                                                                                                                                 cATTGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGACAGATCCCAAt             >  2:885506/1‑151 (MQ=255)
                                                                                                                                    tGCCGCCACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGACAGATCCCAATGCg          >  2:1075333/1‑151 (MQ=255)
                                                                                                                                     gccgccACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGACAGATCCCAATGCgg         <  1:740046/151‑1 (MQ=255)
                                                                                                                                          cACGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGACAGATCCCAATGCGGCTGtt    >  2:583708/1‑151 (MQ=255)
                                                                                                                                           aCGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGACAGATCCCAATGCGGCTGttt   <  1:988850/151‑1 (MQ=255)
                                                                                                                                            cGTAGCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGACAGATCCCAATGCGGCTGTTta  <  1:1234621/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGAAGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGACAGATCCCAATGCGGCTGTTTA  >  NC_000913/302975‑303265

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: