Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,929,028 A→G 8.1% intergenic (‑189/+20) 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,0280AG 8.1% 246.9 / 9.4 86intergenic (‑189/+20)ptrB/yebEprotease II/DUF533 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (23/56);  new base G (4/3);  total (27/59)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

CCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGCGTTCCATAAAATGGTCAATATCAA  >  NC_000913/1928888‑1929175
                                                                                                                                            |                                                                                                                                                   
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:1169787/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)
              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                                                                                             >  1:106039/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)
                                                                                                                                   gTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGCGTTCCATAAAATGGtca         <  2:1513928/150‑1 (MQ=255)
                                                                                                                                   gTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGCGTTCCATAAAATGGtca         >  1:815148/1‑150 (MQ=255)
                                                                                                                                    tCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGCGTTCCATAAAATGGtcaa        >  2:1520114/1‑150 (MQ=255)
                                                                                                                                       aCTTCAGGCTATCCAGAGAGCGGAATTAATCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGCGTTCCATAAAATGGtcaatat     <  1:1427281/150‑1 (MQ=255)
                                                                                                                                        cTTCAGGCTATCCAGAGAGCGGCATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATTGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGCGTTCCATAAAATGGtcaatatc    >  2:726043/1‑150 (MQ=255)
                                                                                                                                        attCAGGCTATCCAGAGAGAGGAATTATTCCGCCAAATTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGCGTTCCATAAAATGGtcaatatc    <  2:837267/149‑1 (MQ=255)
                                                                                                                                         ttgaGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGCGTTCCATAAAATGGtcaatatca   <  2:126585/147‑1 (MQ=255)
                                                                                                                                          tCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGCGTTCCATAAAATGGtcaatatcaa  <  1:590032/150‑1 (MQ=255)
                                                                                                                                          tCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGCGTTCCATAAAATGGtcaatatcaa  >  1:1035024/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
CCGGGGATGCAAGCGAAACAGGGGAGTCATTGCTTAGATGATGACAGGTAATGGCGCGGATATCGAATGTTATGCAGACAGAGAAAATCAGCCTGTTCAAATGGCTGTGCGATTCTGGATAGCCCGAAATAGTCAACTTCAGGCTATCCAGAGAGCGGAATTATTCCGCCAAAGTGCGTTTTTGCTGTTCGAGATCGCGTTCAATGCCGTCACGAACATCCTGGGGGATTTTCAGCGCGTCACCCAGTGCATTCAGGTAACTGCGTTCCATAAAATGGTCAATATCAA  >  NC_000913/1928888‑1929175

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: