Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,929,012 C→T 9.3% intergenic (‑173/+36) ptrB ← / ← yebE protease II/DUF533 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,929,0120CT 9.3% 251.7 / 9.0 85intergenic (‑173/+36)ptrB/yebEprotease II/DUF533 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (23/54);  new base T (4/4);  total (27/58)
Fisher's exact test for biased strand distribution p-value = 2.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGCGTTCCATA  >  NC_000913/1928863‑1929159
                                                                                                                                                     |                                                                                                                                                   
cATCTTATAAGGGTGGCAACACTCACCGGGGATGCAAGCGAAACCGGGCAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGccc                                                                                                                                                     <  1:288552/150‑1 (MQ=255)
cATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGccc                                                                                                                                                     <  2:1527685/150‑1 (MQ=255)
cATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGccc                                                                                                                                                     <  1:869745/150‑1 (MQ=255)
    ttATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGaaa                                                                                                                                                 >  1:1317462/1‑150 (MQ=255)
    ttATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGaaa                                                                                                                                                 <  1:1789170/150‑1 (MQ=255)
       taAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTGTGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAg                                                                                                                                              <  1:1594311/150‑1 (MQ=255)
       taAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAg                                                                                                                                              <  2:1386308/150‑1 (MQ=255)
       taAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAg                                                                                                                                              <  1:533781/150‑1 (MQ=255)
        aaGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGTGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCATCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATatt                                                                                                                                             <  1:1353281/150‑3 (MQ=255)
          gggTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCa                                                                                                                                           <  1:445754/150‑1 (MQ=255)
               gCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTc                                                                                                                                      >  1:1706309/1‑150 (MQ=255)
               gCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTc                                                                                                                                      >  1:1699393/1‑150 (MQ=255)
               gCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTc                                                                                                                                      >  2:590032/1‑150 (MQ=255)
                         ccGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCa                                                                                                                            >  1:397184/1‑150 (MQ=255)
                          cGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCTGAAGTTGACTATTTCGGGCTATCCag                                                                                                                           <  2:688766/150‑1 (MQ=38)
                           ggggATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCTGAAGTTGACTATTTCGGGCTATCCaga                                                                                                                          <  2:6368/150‑1 (MQ=38)
                            gggATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCagag                                                                                                                         <  2:1582541/150‑1 (MQ=255)
                                 gCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCgg                                                                                                                    >  2:676687/1‑150 (MQ=255)
                                 gCAAGCGAAACAGGCGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAAGGTTATGCAGACAGAGAACATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCTGAAGTTGACTATTTCGGGCTATCCagaatcgca                                                                                                                    <  1:273309/150‑7 (MQ=37)
                                  cAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGa                                                                                                                   >  2:815103/1‑150 (MQ=255)
                                  cAAGCGAAACAGGGGAGTCATTGCTTAGATGATGAAAGTTAATGGCGAGGATATCGAATGTTATGCAGACAGAGAAAATCAGACTGTTCAAATGGCTGTGCGATTCTGCATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGAGGa                                                                                                                   >  1:550589/1‑150 (MQ=255)
                                   aaGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGCGCGGaa                                                                                                                  >  1:837267/1‑150 (MQ=255)
                                    aGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAat                                                                                                                 >  2:724361/1‑150 (MQ=255)
                                    aGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAat                                                                                                                 <  1:1680475/150‑1 (MQ=255)
                                    aGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAat                                                                                                                 <  1:306372/150‑1 (MQ=255)
                                     gCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAatt                                                                                                                <  1:1690504/150‑1 (MQ=255)
                                     gCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAatt                                                                                                                <  1:973348/150‑1 (MQ=255)
                                     gCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAatt                                                                                                                <  2:264021/150‑1 (MQ=255)
                                     gCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAatt                                                                                                                <  1:1169787/150‑1 (MQ=255)
                                       gAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAattat                                                                                                              >  1:1513928/1‑150 (MQ=255)
                                           cAGGTGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCg                                                                                                          <  2:36362/150‑1 (MQ=255)
                                              gggAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCATCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCTGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCa                                                                                                       <  1:476318/150‑1 (MQ=255)
                                                 cgTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGCTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAg                                                                                                    <  1:398338/149‑1 (MQ=255)
                                                   tCATTGCTTAGATGATGACAGGTAATGGCGAGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTg                                                                                                  <  2:955237/150‑1 (MQ=255)
                                                     aTTTCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCg                                                                                                <  2:1280656/150‑1 (MQ=255)
                                                     aTTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCg                                                                                                >  2:1347448/1‑150 (MQ=255)
                                                           tAGATGATGTCAGGTAATGTAGAGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTTTTCATATTTCTGTGCGTTTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTg                                                                                          <  2:215198/150‑1 (MQ=255)
                                                            tgatgatGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCCGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGc                                                                                         <  2:1186290/149‑1 (MQ=255)
                                                            aGATGATGACAGGTAATGGCGCGGATATCGAATGTTATTCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTTGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGc                                                                                         <  2:1330964/150‑1 (MQ=255)
                                                            aGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGc                                                                                         <  2:1347189/150‑1 (MQ=255)
                                                               tgatgaCAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGt                                                                                      >  2:1399801/1‑150 (MQ=255)
                                                                gatgaCAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAACTGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGtt                                                                                     <  2:1698866/150‑1 (MQ=255)
                                                                gatgaCAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATATCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGtt                                                                                     <  2:1692038/150‑1 (MQ=255)
                                                                    aCAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCgag                                                                                 <  2:1291522/150‑1 (MQ=255)
                                                                     cAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCgaga                                                                                <  1:1693397/150‑1 (MQ=255)
                                                                        gTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATcg                                                                             >  2:504096/1‑150 (MQ=255)
                                                                        gTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATcg                                                                             <  2:848935/150‑1 (MQ=255)
                                                                        gTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATcg                                                                             >  1:106039/1‑150 (MQ=255)
                                                                          aaTGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATcgcg                                                                           <  1:376140/150‑1 (MQ=255)
                                                                            tGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGtt                                                                         <  1:43967/150‑1 (MQ=255)
                                                                            tGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAAAAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGtt                                                                         <  2:50492/150‑1 (MQ=255)
                                                                                gcgGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAAt                                                                     <  2:1826207/150‑1 (MQ=255)
                                                                                   gATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGcc                                                                  <  1:1688386/150‑1 (MQ=255)
                                                                                      atCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTc                                                               >  1:1172602/1‑150 (MQ=255)
                                                                                      atCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTc                                                               <  1:1625843/150‑1 (MQ=255)
                                                                                          aaTGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTGAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGa                                                           <  2:885939/150‑1 (MQ=255)
                                                                                          aaTGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGa                                                           <  2:740121/150‑1 (MQ=255)
                                                                                            tGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAAc                                                         <  1:697789/150‑1 (MQ=255)
                                                                                             tttATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACa                                                        <  1:1609866/149‑1 (MQ=255)
                                                                                                   cagacagaGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTgg                                                  <  1:535071/150‑1 (MQ=255)
                                                                                                    agacagaGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTggg                                                 <  2:1526651/150‑1 (MQ=255)
                                                                                                     gacagaGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGTAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGGTCAATGCCGGCACGAACATCCTgggg                                                <  2:931112/150‑1 (MQ=255)
                                                                                                      acagaGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTggggg                                               <  2:23804/150‑1 (MQ=255)
                                                                                                       cagaGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGa                                              >  2:187629/1‑150 (MQ=255)
                                                                                                        agagAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGTATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGAGTTCAATGCCGTCACGAACATCCTGGGGGAt                                             >  2:1631026/1‑150 (MQ=255)
                                                                                                            aaaaTCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGCGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTc                                         <  2:1617376/150‑1 (MQ=255)
                                                                                                            aaaaTCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTc                                         <  2:1764982/150‑1 (MQ=255)
                                                                                                              aaTCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAg                                       <  2:1270539/150‑1 (MQ=255)
                                                                                                               aTCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAgc                                      <  1:772435/150‑1 (MQ=255)
                                                                                                                tCAGCCTGTGCAAATGGCGGTGCGATGCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTCTTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAgcg                                     <  1:938536/150‑1 (MQ=255)
                                                                                                                 cAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAgcgc                                    >  1:1572151/1‑150 (MQ=255)
                                                                                                                  cgCCTGTTCAAATGGCTGTGCGAATCTTGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGTGAGCGGAATTATTCCGGCAAAGTGCGTTTTTGCTGTTCGATATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAgcgcc                                   <  1:1045080/149‑2 (MQ=255)
                                                                                                                  aGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGcgcg                                   <  2:46483/150‑1 (MQ=255)
                                                                                                                   gCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGt                                  >  2:739483/1‑150 (MQ=255)
                                                                                                                       gTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCAcc                              <  1:13670/150‑1 (MQ=255)
                                                                                                                       gTTCAAATGGCTGTGCGATGCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCAcc                              <  1:1502413/150‑1 (MQ=255)
                                                                                                                            aaTGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTg                         <  1:1085626/150‑1 (MQ=255)
                                                                                                                             aTGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGc                        <  1:504096/150‑1 (MQ=255)
                                                                                                                                  tgtgCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCa                   <  1:332386/150‑1 (MQ=255)
                                                                                                                                    tgCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAgg                 <  1:1347448/150‑1 (MQ=255)
                                                                                                                                    tgCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAgg                 >  2:1015342/1‑150 (MQ=255)
                                                                                                                                      cGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCAACCAGTGCATTCAGGTa               >  1:256296/1‑150 (MQ=255)
                                                                                                                                         ctCTGGATAGCCTGAAGTTGACTATTTCGGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACAAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACt            >  1:730806/2‑150 (MQ=38)
                                                                                                                                           cTGGATAGCCTGAAGTTGACTATTTCGGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGc          >  2:283975/1‑150 (MQ=38)
                                                                                                                                           cTGGATAGCCTGAAGTTGACTATTTCGGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGAGTCACCCAGTGCATTCAGGTACCTGc          >  1:103298/1‑150 (MQ=38)
                                                                                                                                            tGGATAGCCTGAAGTTGACTATTTCGGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCCCCCAGTGCATTCAGGTAACTGCg         >  1:920708/1‑150 (MQ=38)
                                                                                                                                                 aGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGCGTTCCa    <  2:1825725/150‑1 (MQ=255)
                                                                                                                                                  gCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTACCGGCGTTCCAt   >  2:870055/1‑150 (MQ=255)
                                                                                                                                                   cccGAAATAGTCAACTTCAGGATATCCAGAGATCGGAATTATTCCTCCAAAGTGCGTTTTTTCTGTTCGATATCGCGTTCATTGCGGTCACGAACATCCTTGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGCGTTCCATa  <  1:635066/150‑1 (MQ=255)
                                                                                                                                                   cccGAAATAGGCACCTTCAGGCTATCCAGAGAGCGGACGTATTCCGCCAAAGTGCGGTTTTGCGGGGCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCCGCGCGTCCCCCAGTGCAACCAGGTAACTGCGTTCCATa  <  1:942263/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CATCTTATAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGCGTTCCATA  >  NC_000913/1928863‑1929159

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: