Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,423,231 T→A 5.5% intergenic (+37/+192) yhdZ → / ← rrfF putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,2310TA 5.5% 363.5 / 5.7 120intergenic (+37/+192)yhdZ/rrfFputative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base T (47/66);  new base A (4/3);  total (51/69)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 122 total aligned reads displayed.

GGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATA  >  NC_000913/3423083‑3423378
                                                                                                                                                    |                                                                                                                                                    
ggTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTc                                                                                                                                                     <  2:1450295/150‑1 (MQ=255)
ggTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCCCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTc                                                                                                                                                     >  2:1612870/1‑150 (MQ=255)
 gTAATTTTTATGGATCGTGGGGAAATAGTGGATCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCg                                                                                                                                                    <  2:768624/150‑1 (MQ=255)
 gTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCg                                                                                                                                                    <  1:454785/150‑1 (MQ=255)
 gTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCg                                                                                                                                                    >  1:594198/1‑150 (MQ=255)
  tAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCgg                                                                                                                                                   <  1:1404217/150‑1 (MQ=255)
  tAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCgg                                                                                                                                                   >  2:1179254/1‑150 (MQ=255)
    aTTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGAt                                                                                                                                                 <  2:1448806/150‑1 (MQ=255)
     tttttATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATg                                                                                                                                                >  1:768021/1‑150 (MQ=255)
       tttATGGATCGTTGGGAAATAGTTGAGCAAGCTGCACCTTATTAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGTCCATCCTTTCGGATGGc                                                                                                                                              <  2:1327669/150‑1 (MQ=255)
         tATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCGAAAGGATGGCCt                                                                                                                                           <  2:39773/150‑1 (MQ=255)
          aTGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCtt                                                                                                                                          <  2:1361869/150‑1 (MQ=255)
          aTGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCGAAAGGATGGCCtt                                                                                                                                          <  2:353100/150‑1 (MQ=255)
          aTGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCGAAAGGATGGCCtt                                                                                                                                          <  1:1434065/150‑1 (MQ=255)
             gATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATCTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCg                                                                                                                                       >  2:898432/1‑150 (MQ=255)
                 gTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCttg                                                                                                                                   <  2:1644812/150‑1 (MQ=255)
                 gTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCttg                                                                                                                                   <  1:533665/150‑1 (MQ=255)
                  tGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCttga                                                                                                                                  <  2:354869/150‑1 (MQ=255)
                    gggAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATg                                                                                                                                                >  2:100372/1‑135 (MQ=255)
                    gggAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATg                                                                                                                                                <  1:100372/135‑1 (MQ=255)
                     ggAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCttgattt                                                                                                                               <  2:1481563/150‑1 (MQ=255)
                      gAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCttgatttg                                                                                                                              >  1:1255181/1‑150 (MQ=255)
                       aaaTAGTGGAGCAAGCTGCACCTGATGTATCTTTTGCGCATCCTAAATCAGAGCGTTCGAGGGCATTTTTATCGCAGGTAATCCATTAATTGATTGTTAGTGCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCttgatttga                                                                                                                             <  1:140861/150‑1 (MQ=255)
                         aTAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATg                                                                                                                           <  1:240406/150‑1 (MQ=255)
                           aGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGta                                                                                                                         >  2:1369197/1‑149 (MQ=255)
                           aGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTc                                                                                                                         <  2:653815/150‑1 (MQ=255)
                              ggAGCCAGCTGCCCCTGATGAAGTGTTTGCGCCTCCTACCTCAGCGCGTACGCGGGCATTTTTATCGCAGGTACTCCCTTACTTGAATGGTAGTTCGAACAGCACAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTgg                                                                                                                      <  1:1074293/150‑1 (MQ=255)
                              ggAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAACAGCATAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTgg                                                                                                                      <  1:1074294/150‑1 (MQ=255)
                               gAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAGAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGc                                                                                                                     <  1:422350/150‑1 (MQ=255)
                               gAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGc                                                                                                                     >  2:1272669/1‑150 (MQ=255)
                               gAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGc                                                                                                                     <  1:922893/150‑1 (MQ=255)
                                aGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCa                                                                                                                    <  2:1134639/150‑1 (MQ=255)
                                   aaGCTTCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGtt                                                                                                                 <  2:545347/150‑1 (MQ=255)
                                     gCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTa                                                                                                               <  1:1669671/150‑1 (MQ=255)
                                      cTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTAt                                                                                                              >  1:94975/1‑150 (MQ=255)
                                       tGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTATCGGATGGCCTATCGCTTGATTTGATGTCTGGCAGTTTATg                                                                                                             >  1:110845/1‑150 (MQ=255)
                                              gatgaATTTTTTGCCCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATggcgggcg                                                                                                      <  1:1612870/150‑1 (MQ=255)
                                               atgaATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGt                                                                                                     <  1:828575/150‑1 (MQ=255)
                                                tgaATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTc                                                                                                    <  1:1582516/150‑1 (MQ=255)
                                                 taaTTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTcc                                                                                                   <  2:1323465/149‑1 (MQ=255)
                                                     tttttGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGcc                                                                                               >  1:1462020/1‑150 (MQ=255)
                                                        ttGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGc                                                                                            >  1:1435615/1‑150 (MQ=255)
                                                                ccTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTcc                                                                                    >  1:3557/1‑150 (MQ=255)
                                                                ccTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTcc                                                                                    <  1:776449/150‑1 (MQ=255)
                                                                 cTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCACAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCg                                                                                   <  1:285926/150‑1 (MQ=255)
                                                                 cTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCg                                                                                   >  1:1657640/1‑150 (MQ=255)
                                                                        agagCGTACGCGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGtt                                                                            <  2:594198/150‑1 (MQ=255)
                                                                              tACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGTCCTTTAGATTTATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGAGCCGTTGCTACg                                                                      >  1:362715/1‑150 (MQ=25)
                                                                                   gggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACg                                                                 >  2:662696/1‑150 (MQ=255)
                                                                                    ggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATATCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGt                                                                <  2:600429/150‑1 (MQ=255)
                                                                                      cATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGAGGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCGAAGCAACggcc                                                              <  2:1180801/150‑4 (MQ=25)
                                                                                       aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCa                                                             <  2:410595/150‑1 (MQ=255)
                                                                                       aTTTTTATCGAAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCa                                                             <  2:466504/150‑1 (MQ=255)
                                                                                          tttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAAt                                                          <  1:1212400/150‑1 (MQ=255)
                                                                                           ttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCAGCCTTTCCGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATc                                                         <  2:268272/150‑1 (MQ=255)
                                                                                             aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCg                                                       >  2:956896/1‑150 (MQ=255)
                                                                                             aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCg                                                       <  1:596564/150‑1 (MQ=255)
                                                                                                   ggTAATCCCTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcg                                                 <  1:1196862/150‑1 (MQ=255)
                                                                                                   ggTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcg                                                 <  1:841212/150‑1 (MQ=255)
                                                                                                       aTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCcggcgg                                             >  1:787334/1‑150 (MQ=255)
                                                                                                          cATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCTCCGGGCCGTTGCTTCGCACCGTTCAAATCCGCTCCCGGCGGAtt                                          >  1:925329/1‑150 (MQ=255)
                                                                                                          cATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCGTGCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt                                          >  2:1472657/1‑150 (MQ=255)
                                                                                                          cATTAATTGAATGTTAGTTAGAAAATCAAAAAGGGCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCATTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCTGCGGAtt                                          <  2:110845/150‑1 (MQ=18)
                                                                                                           aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGAttt                                         >  2:1531237/1‑150 (MQ=255)
                                                                                                              aaTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTc                                      >  1:780158/1‑150 (MQ=255)
                                                                                                               aTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTcc                                     <  2:1010650/150‑1 (MQ=255)
                                                                                                                 tGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGTTCCCGGCGGATTTGTCCTa                                   >  2:143115/1‑150 (MQ=255)
                                                                                                                 tGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTa                                   >  1:1248971/1‑150 (MQ=255)
                                                                                                                  gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTAc                                  >  1:1103664/1‑150 (MQ=255)
                                                                                                                   aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACt                                 >  2:1074302/1‑150 (MQ=255)
                                                                                                                   aaTGGTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGGCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACt                                 <  2:1372987/150‑1 (MQ=255)
                                                                                                                     tGTTAGTTCGAAAAGCAACAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCg                               <  2:282607/150‑1 (MQ=255)
                                                                                                                     tGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCg                               <  2:94975/150‑1 (MQ=255)
                                                                                                                      gTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCgg                              <  2:1435615/150‑1 (MQ=255)
                                                                                                                      gTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCgg                              >  1:303369/1‑150 (MQ=255)
                                                                                                                       ttAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCggg                                                                                 <  2:1230970/98‑1 (MQ=255)
                                                                                                                       ttAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCggg                                                                                 >  1:1230970/1‑98 (MQ=255)
                                                                                                                         aGTTCGAAAAGCAAAAAGGCCATCCTTTGGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGgag                           <  2:1496144/150‑1 (MQ=255)
                                                                                                                         aGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGgag                           <  2:138138/150‑1 (MQ=255)
                                                                                                                         aGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGgag                           <  2:1678390/150‑1 (MQ=255)
                                                                                                                          gTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGgaga                          >  2:1671362/1‑150 (MQ=255)
                                                                                                                           ttCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGagag                         <  2:1103664/150‑1 (MQ=255)
                                                                                                                             cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAgtg                       >  2:12224/1‑150 (MQ=255)
                                                                                                                                aaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTc                    <  1:898432/150‑1 (MQ=255)
                                                                                                                                  aGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCAc                  >  2:1019815/1‑150 (MQ=255)
                                                                                                                                   agcgaaaGGCCATCCGAAAGGATGGCCTTTTTGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCAc                  >  2:1333803/5‑150 (MQ=11)
                                                                                                                                     aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCgc               >  1:1486540/1‑149 (MQ=255)
                                                                                                                                       aaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGaca             <  2:831071/150‑1 (MQ=255)
                                                                                                                                        aaGGCCATCCGAAAGGATGGCCTTTTTGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGaca             >  2:158063/1‑150 (MQ=11)
                                                                                                                                         aGGCCATCCGAAAGGATGGCCTTTTTGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGGTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGGCCTACTCGGGAGAGTGTGCCCCGccaa            >  1:408253/1‑150 (MQ=11)
                                                                                                                                          ggCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaac          >  2:1413769/1‑150 (MQ=255)
                                                                                                                                           gCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaaca         >  2:1325410/1‑150 (MQ=255)
                                                                                                                                           gCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaaca         >  1:1478065/1‑150 (MQ=255)
                                                                                                                                            ccATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCAcc                 >  2:927321/1‑141 (MQ=255)
                                                                                                                                              aTCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAacaaca      >  2:76892/1‑150 (MQ=255)
                                                                                                                                              aTCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAacaaca      <  2:445053/150‑1 (MQ=255)
                                                                                                                                                ccTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGa    <  2:787334/150‑1 (MQ=255)
                                                                                                                                                ccTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGa    <  1:987347/150‑1 (MQ=255)
                                                                                                                                                  tttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATa  <  2:1235997/150‑1 (MQ=255)
                                                                                                                                                  tttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATa  >  2:515922/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGTAATTTTTATGGATCGTGGGGAAATAGTGGAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATA  >  NC_000913/3423083‑3423378

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