Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,929,018 A→T 6.0% 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 6.0% 294.3 / 5.5 99intergenic (‑179/+30)ptrB/yebEprotease II/DUF533 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (39/54);  new base T (2/4);  total (41/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

TAAGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGCGTTCCATAAAAT  >  NC_000913/1928870‑1929163
                                                                                                                                                    |                                                                                                                                                 
taAGGGTGGCAAAACTCACCTGGGATGCAAGCGAAACAGGGGATTCATTGCTTATATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAg                                                                                                                                                  <  2:700752/150‑1 (MQ=255)
  aGGGTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTc                                                                                                                                                <  2:432503/150‑1 (MQ=255)
   gggTGGCAAAACTCACCGGGGATGCAAGCGAAACAGGTGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCa                                                                                                                                               <  1:1345826/150‑1 (MQ=255)
       ggCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACtt                                                                                                                                           <  2:1304153/150‑1 (MQ=255)
        gCAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTc                                                                                                                                          >  1:1568492/1‑150 (MQ=255)
         cAAAACTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCa                                                                                                                                         <  2:1745650/150‑1 (MQ=255)
          aaaaCTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAg                                                                                                                                        >  1:969873/1‑150 (MQ=255)
          aaaaCTCACCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAg                                                                                                                                        <  1:581064/150‑1 (MQ=255)
                 aCCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCTGAAGTTGACTATTTCGGGCTATcc                                                                                                                                 <  1:2098095/150‑1 (MQ=39)
                  ccGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCTGAAGTTGACTATTTCGGGCTATCCa                                                                                                                                <  1:1606012/150‑1 (MQ=39)
                  ccGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCa                                                                                                                                <  1:2098996/150‑1 (MQ=255)
                  ccGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTCTTCAAATGGCTGTGCGATTCTGGATAGCCTGAAGTTGACTATTTCGGGCTATCCa                                                                                                                                <  1:1970016/150‑1 (MQ=38)
                   cGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGGCTGTTCAAATGGCTGTGCGATTCTGGATAGCCTGAAGTTGACTATTTCGGGCTATCCag                                                                                                                               <  2:317765/150‑1 (MQ=38)
                     gggATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCagag                                                                                                                             <  1:1907894/150‑1 (MQ=255)
                         tGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCg                                                                                                                         >  2:919981/1‑150 (MQ=255)
                         tGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCg                                                                                                                         <  1:28394/150‑1 (MQ=255)
                          gCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCCGGCTATCCAGAGAGCgg                                                                                                                        <  1:2078347/150‑1 (MQ=255)
                          gCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCgg                                                                                                                        >  2:1411428/1‑150 (MQ=255)
                          gCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCgg                                                                                                                        <  2:1985523/150‑1 (MQ=255)
                           cAAGCGAAACAGGGGAGTCATTGCTTAGATGTTGACAGGTAATGGCGCCGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTTTGCGATTCTGGATAGCCCGAAATAGTCTACTTCAGTCTAACAATAGAGCTGa                                                                                                                       >  1:323110/1‑150 (MQ=255)
                           cAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGa                                                                                                                       <  1:2078934/150‑1 (MQ=255)
                           cAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGa                                                                                                                       >  1:1394979/1‑150 (MQ=255)
                             aGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCACCTTCAGGCTATCCAGAGAGCGGAat                                                                                                                     >  1:340978/1‑150 (MQ=255)
                             aGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAat                                                                                                                     >  2:1389066/1‑150 (MQ=255)
                              gCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAatt                                                                                                                    >  1:1224033/1‑150 (MQ=255)
                                gAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATAGAATGATATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGAATCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAattat                                                                                                                  >  1:820943/1‑150 (MQ=255)
                                   aCAGGGGAGTCATTGCTTAGATGATGACATGTAATGGCTCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATGCTGGATAGCCCGAAATAGTCCACTTCAGGCTATCCAGAGAGCGGAATTATTcc                                                                                                               <  2:916783/150‑1 (MQ=255)
                                    cAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCg                                                                                                              <  2:1660104/150‑1 (MQ=255)
                                        ggAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCaa                                                                                                          <  1:645098/150‑1 (MQ=255)
                                              aTTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCg                                                                                                    >  1:1824457/1‑150 (MQ=255)
                                                 gCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGttt                                                                                                 >  1:1722661/1‑150 (MQ=255)
                                                  cTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGtttt                                                                                                >  1:514435/1‑150 (MQ=255)
                                                  cTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGtttt                                                                                                >  1:593191/1‑150 (MQ=255)
                                                     aGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGc                                                                                             <  2:2012897/150‑1 (MQ=255)
                                                           tgaCAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCg                                                                                       <  2:1225170/150‑1 (MQ=255)
                                                           tgaCAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCg                                                                                       <  2:953364/150‑1 (MQ=255)
                                                            gaCAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCga                                                                                      <  1:775853/150‑1 (MQ=255)
                                                            gaCAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCga                                                                                      <  2:1212695/150‑1 (MQ=255)
                                                             aCAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCgag                                                                                     <  1:1905725/150‑1 (MQ=255)
                                                             aCAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCgag                                                                                     <  1:1981134/150‑1 (MQ=255)
                                                              cAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCgaga                                                                                    >  2:2111380/1‑150 (MQ=255)
                                                               aGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGAt                                                                                   >  1:587707/1‑150 (MQ=255)
                                                               aGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTACCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTGGCTGTTCGAGAt                                                                                   <  1:1054728/150‑1 (MQ=255)
                                                                ggTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATc                                                                                  >  1:1080767/1‑150 (MQ=255)
                                                                      ggCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTc                                                                            <  1:2057208/150‑1 (MQ=255)
                                                                      ggCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTc                                                                            >  2:7838/1‑150 (MQ=255)
                                                                      ggCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTCCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTc                                                                            <  2:1851927/150‑1 (MQ=255)
                                                                       gcgcGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCa                                                                           <  2:486178/150‑1 (MQ=255)
                                                                       gcgcGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCa                                                                           <  2:1183124/150‑1 (MQ=255)
                                                                        cgcgGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCaa                                                                          >  1:857511/1‑150 (MQ=255)
                                                                        cgcgGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAACGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCaa                                                                          >  1:1305853/1‑150 (MQ=255)
                                                                            gATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGcc                                                                      <  2:1007982/150‑1 (MQ=255)
                                                                              tatCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGt                                                                    <  2:595705/150‑1 (MQ=255)
                                                                                tCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCa                                                                  <  2:1824457/150‑1 (MQ=255)
                                                                                 cGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCAc                                                                 >  2:1753814/1‑150 (MQ=255)
                                                                                 cGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCAc                                                                 >  2:1501847/1‑150 (MQ=255)
                                                                                  gAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACg                                                                >  2:1143402/1‑150 (MQ=255)
                                                                                    aTGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGaa                                                              <  2:996926/150‑1 (MQ=255)
                                                                                     tGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCTCGTTCAATGCCGTCACGAAc                                                             >  2:1373294/1‑150 (MQ=255)
                                                                                        tATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATc                                                          >  2:752819/1‑150 (MQ=255)
                                                                                         aTGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGCTCGCGTTCAGTGCCGTCACGAACATcc                                                         >  2:575719/1‑150 (MQ=255)
                                                                                          tGCAGACAGAGAAAATCAGCCTGTTCTAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCt                                                        <  2:2097035/150‑1 (MQ=255)
                                                                                           gCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTg                                                       <  2:1221582/150‑1 (MQ=255)
                                                                                            cagacagaGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTgg                                                      <  2:1378357/150‑1 (MQ=255)
                                                                                            cagacagaGAAAATCAGACTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCCTTTTTGCAGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTgg                                                      >  1:375849/1‑150 (MQ=255)
                                                                                              gacagaGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAAATATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTgggg                                                    <  1:67060/150‑1 (MQ=255)
                                                                                                cagaGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGa                                                  <  1:845000/150‑1 (MQ=255)
                                                                                                 agagAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGACAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGAt                                                 >  2:1863784/1‑150 (MQ=255)
                                                                                                  gagaAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGAtt                                                >  1:988315/1‑150 (MQ=255)
                                                                                                  gagaAAAACAGCCTGTTCAAATGGCTGTGAGATTCTGGATAGACCGAAATAGTCAACTTCAGGCTATACAGAGAGCGGAATTATTCCTCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGAtt                                                >  2:939117/1‑150 (MQ=255)
                                                                                                   agaAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGAttt                                               >  2:1274128/1‑150 (MQ=255)
                                                                                                    gaAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTTCTGTTCGAGATCGCGTTCAATGCCGTCACGAAAATCCTGGGGGAtttt                                              >  1:1766566/1‑150 (MQ=255)
                                                                                                     aaaaTCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTc                                             <  2:2091994/150‑1 (MQ=255)
                                                                                                     aaaaTCAGCCTGTACAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTACGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTc                                             <  1:986254/150‑1 (MQ=255)
                                                                                                         tCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAgcg                                         <  1:1039955/150‑1 (MQ=255)
                                                                                                          cAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAgcgc                                        >  2:813024/1‑150 (MQ=255)
                                                                                                           aGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGcgcg                                       <  2:1879405/150‑1 (MQ=255)
                                                                                                                    aaaTGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGt                              >  1:713017/1‑150 (MQ=255)
                                                                                                                     aaTGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTg                             <  1:1610820/150‑1 (MQ=255)
                                                                                                                       tGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCa                           <  2:1394979/150‑1 (MQ=255)
                                                                                                                         gCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCAtt                         >  1:1940546/1‑150 (MQ=255)
                                                                                                                          cTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTc                        <  2:969873/150‑1 (MQ=255)
                                                                                                                          cTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTc                        <  2:514435/150‑1 (MQ=255)
                                                                                                                          cTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTc                        >  1:1977727/1‑150 (MQ=255)
                                                                                                                           tgtgCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCa                       <  2:1568492/150‑1 (MQ=255)
                                                                                                                            gtgCGATTCTGGATAGCCCGAAATAGTCAACTGCAGGCTATCCAGAGAGCGGAATTATTCCGACAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACTAACATCCTGGGGGATTTTCAGCTCGTCACCCAGTGCATTCAg                      >  2:1995783/1‑150 (MQ=255)
                                                                                                                             tgCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAgg                     <  2:1873802/150‑1 (MQ=255)
                                                                                                                             tgCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAgg                     <  2:347775/150‑1 (MQ=255)
                                                                                                                               cGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCTGTGCATTCAGGTa                   <  1:976781/150‑1 (MQ=255)
                                                                                                                                gATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTaa                  >  2:1531387/1‑150 (MQ=255)
                                                                                                                                   tCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTg               <  2:1061955/150‑1 (MQ=255)
                                                                                                                                    cTGGATAGCCTGAAGTTGACTATTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCAGTCAGGTAACTGc              >  2:835464/1‑150 (MQ=255)
                                                                                                                                    cTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGc              <  2:2093396/150‑1 (MQ=255)
                                                                                                                                     tGGATAGCCTGAAGTTGACTATTTCGGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGCg             >  2:896018/1‑150 (MQ=39)
                                                                                                                                     tGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATTCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGCg             <  2:1080767/150‑1 (MQ=255)
                                                                                                                                       gATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTTCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGCGtt           <  1:1411428/150‑1 (MQ=255)
                                                                                                                                            cccGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGCGTTCCATa      <  2:2085389/150‑1 (MQ=255)
                                                                                                                                               gAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATAGCGTTCAATGCAGTCCCGAACATCATGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGCGTTCCATaaaa   >  1:147981/1‑150 (MQ=255)
                                                                                                                                                aaaTAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGCGTTCCATAAAAt  <  1:919981/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: