Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,929,018 A→T 7.8% intergenic (‑179/+30) 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,0180AT 7.8% 241.1 / 6.3 85intergenic (‑179/+30)ptrB/yebEprotease II/DUF533 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (22/56);  new base T (4/3);  total (26/59)
Fisher's exact test for biased strand distribution p-value = 1.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGCGTTCCATAAAAT  >  NC_000913/1928870‑1929163
                                                                                                                                                    |                                                                                                                                                 
taAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTGTGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAg                                                                                                                                                  <  1:1594311/150‑1 (MQ=255)
taAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAg                                                                                                                                                  <  1:533781/150‑1 (MQ=255)
taAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAg                                                                                                                                                  <  2:1386308/150‑1 (MQ=255)
 aaGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGTGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCATCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATatt                                                                                                                                                 <  1:1353281/150‑3 (MQ=255)
   gggTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCa                                                                                                                                               <  1:445754/150‑1 (MQ=255)
        gCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTc                                                                                                                                          >  2:590032/1‑150 (MQ=255)
        gCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTc                                                                                                                                          >  1:1699393/1‑150 (MQ=255)
        gCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTc                                                                                                                                          >  1:1706309/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                                                                                                                     <  1:306372/150‑1 (MQ=255)
                             aGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAat                                                                                                                     >  2:724361/1‑150 (MQ=255)
                             aGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAat                                                                                                                     <  1:1680475/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                                                                                 >  1:106039/1‑150 (MQ=255)
                                                                 gTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATcg                                                                                 >  2:504096/1‑150 (MQ=255)
                                                                 gTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATcg                                                                                 <  2:848935/150‑1 (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)
                                                                                                                                                aaaTAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGCGTTCCATAAAAt  <  2:395298/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGCGTTCCATAAAAT  >  NC_000913/1928870‑1929163

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: