Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,770,390 A→G 6.2% noncoding (19/106 nt) rprA → sRNA antisense activator for RpoS translation, Hfq‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,770,3900AG 6.2% 336.7 / 9.2 113noncoding (19/106 nt)rprAsRNA antisense activator for RpoS translation, Hfq‑dependent
Reads supporting (aligned to +/- strand):  ref base A (63/43);  new base G (4/3);  total (67/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.14e-01

Only 100 of 115 total aligned reads displayed.

TTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGTTGCACAAGATACT  >  NC_000913/1770243‑1770538
                                                                                                                                                   |                                                                                                                                                    
ttgctCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATAtt                                                                                                                                                    <  2:207954/150‑1 (MQ=255)
  gctCTTCCACATATTTTCCTGCTTAACGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAATATGTTGa                                                                                                                                                  <  1:361464/150‑1 (MQ=255)
   ctctTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGAt                                                                                                                                                 >  2:1055761/1‑150 (MQ=255)
   ctctTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGAt                                                                                                                                                 >  1:682263/1‑150 (MQ=255)
    tctTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGAtt                                                                                                                                                <  1:158954/150‑1 (MQ=255)
     ctTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGAttt                                                                                                                                               <  2:645179/150‑1 (MQ=255)
      ttCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTta                                                                                                                                              <  2:697724/150‑1 (MQ=255)
        ccccATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTtata                                                                                                                                            >  1:1192755/1‑150 (MQ=255)
         cccATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATaa                                                                                                                                           >  2:1370250/1‑150 (MQ=255)
          ccATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAATATGTTGATTTATaac                                                                                                                                          <  1:452079/150‑2 (MQ=255)
           cATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAATATGTTGATTTATaacc                                                                                                                                         <  2:412661/150‑3 (MQ=255)
           cATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGGCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGc                                                                                                                                         <  1:689362/150‑1 (MQ=255)
           cATATTTACATCCTTACCGCCATTCAGCTCGTAGTACCTGTCGCAAATTCTTAACAGTTTTTAAACTAATGTCACGAATAAGATAGACGCATAAAGTCCGTATGCCTACTATTAGCTCACGCTTATAAATAAACATATTTATTTAAAAGc                                                                                                                                         >  2:880059/1‑150 (MQ=255)
              atTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATg                                                                                                                                      <  2:1075045/150‑1 (MQ=255)
                tttCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGa                                                                                                                                    <  2:665963/150‑1 (MQ=255)
                  tCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGaaa                                                                                                                                  <  1:1622046/150‑1 (MQ=255)
                   ccTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAAt                                                                                                                                 >  1:1546608/1‑150 (MQ=255)
                   ccTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAAt                                                                                                                                 >  1:1595778/1‑150 (MQ=255)
                    cTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATc                                                                                                                                >  1:67893/1‑150 (MQ=255)
                    cTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATc                                                                                                                                >  2:304561/1‑150 (MQ=255)
                    cTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGCTTATAAATCAACATATTGATTTTTAAGCATGGAAATc                                                                                                                                >  2:299890/1‑150 (MQ=255)
                           ccGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAGAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAg                                                                                                                         >  2:633545/1‑150 (MQ=255)
                           ccGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAg                                                                                                                         >  1:706673/1‑150 (MQ=255)
                           ccGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAg                                                                                                                         >  2:1474715/1‑150 (MQ=255)
                             gCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTg                                                                                                                       <  1:932459/150‑1 (MQ=255)
                               cATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGaa                                                                                                                     >  1:504318/1‑150 (MQ=255)
                                aTTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGaaa                                                                                                                    <  2:816510/150‑1 (MQ=255)
                                 ttCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAaac                                                                                                                   >  1:1250639/1‑150 (MQ=255)
                                  tCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAaaca                                                                                                                  >  1:415425/1‑150 (MQ=255)
                                     gCTGGTAGTACCTGTCGCAAATTATTTACAGTTTTTAAAATAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACg                                                                                                               <  1:1173646/150‑1 (MQ=255)
                                      cTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGa                                                                                                              >  2:471016/1‑150 (MQ=255)
                                         gtagtaCCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTTAAACAACGAAtt                                                                                                           >  2:666935/1‑150 (MQ=255)
                                         gtagtaCCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAAtt                                                                                                           >  2:218337/1‑150 (MQ=255)
                                               ccTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCtgtg                                                                                                     >  2:1066040/1‑150 (MQ=255)
                                                     gaaaaTTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGt                                                                                               <  2:1546608/148‑1 (MQ=255)
                                                     gCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACTAATTGCTGTGTGTAga                                                                                               >  1:1140393/1‑149 (MQ=255)
                                                      cAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTc                                                                                              <  2:682263/150‑1 (MQ=255)
                                                                tcattttttAAACTAATGAGAAGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCAtc                                                                                    <  1:377219/146‑1 (MQ=255)
                                                                 cAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCACCATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAAGTGCTGTGTGTAGTCTTTGCCCAtcg                                                                                   >  1:937651/1‑149 (MQ=255)
                                                                 cAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCAtct                                                                                   >  1:773292/1‑150 (MQ=255)
                                                                   gTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTcc                                                                                 >  1:1254436/1‑150 (MQ=255)
                                                                     ttttAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCa                                                                               >  2:1357676/1‑150 (MQ=255)
                                                                          aaCTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATg                                                                          <  2:706673/150‑1 (MQ=255)
                                                                           aCTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATgg                                                                         >  2:654133/1‑150 (MQ=255)
                                                                              aaTGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCt                                                                      >  2:1000828/1‑150 (MQ=255)
                                                                                    aCGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCttttttt                                                                <  1:951912/150‑1 (MQ=255)
                                                                                       aaTCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTaa                                                             <  1:1305157/150‑1 (MQ=255)
                                                                                           tGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGAAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACAtt                                                         >  2:1279620/1‑150 (MQ=255)
                                                                                            gATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACAttt                                                        >  1:528659/1‑150 (MQ=255)
                                                                                               cGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTc                                                     <  2:1675036/150‑1 (MQ=255)
                                                                                                gACGCAATAATTCCGTATGCCTACTATTAGCTCTCGGTTATAAATCATCATATTGTTTTATAAGCATGGAAATCCCCTGAGTGAATCAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTcc                                                    <  1:299890/150‑1 (MQ=255)
                                                                                                gACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTcc                                                    <  1:304561/150‑1 (MQ=255)
                                                                                                  cGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTCTAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGc                                                  <  2:1577254/150‑1 (MQ=255)
                                                                                                  cGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGc                                                  >  2:365980/1‑150 (MQ=255)
                                                                                                   gCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCt                                                 >  2:124745/1‑150 (MQ=255)
                                                                                                   gCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCt                                                 >  1:1073691/1‑150 (MQ=255)
                                                                                                     aaaaaGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATCGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCAACGATGGGCTTTTTTTTAACATTTTTCCGCtta                                               >  2:571736/1‑149 (MQ=255)
                                                                                                      aaaaGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCg                                              <  1:1332046/150‑1 (MQ=255)
                                                                                                      aaaaGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCg                                              >  2:1331958/1‑150 (MQ=255)
                                                                                                      aaaaGTCCGTATGCCTAATATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCg                                              <  1:534350/150‑1 (MQ=255)
                                                                                                       aaaGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGc                                             >  1:605705/1‑150 (MQ=255)
                                                                                                        aaGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCt                                            <  1:1000828/150‑1 (MQ=255)
                                                                                                        aaGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCt                                            <  1:1177305/150‑1 (MQ=255)
                                                                                                          gTCCGTACGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACAAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCATTATTTTAACATATTTCCTCTTCGCAAc                                          >  2:406846/1‑150 (MQ=255)
                                                                                                           tCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTAcc                                         <  1:339978/150‑1 (MQ=255)
                                                                                                              gTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCg                                      >  1:363316/1‑150 (MQ=255)
                                                                                                               tATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGc                                     >  1:1578604/1‑150 (MQ=255)
                                                                                                                aTGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGcc                                    <  2:1595778/150‑1 (MQ=255)
                                                                                                                 tGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGccc                                   >  1:770443/1‑150 (MQ=255)
                                                                                                                  gCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGcccc                                  >  2:1695053/1‑150 (MQ=255)
                                                                                                                   ccTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCt                                 >  2:1694410/1‑150 (MQ=255)
                                                                                                                   ccTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTCCGCTTCGCTACCTCGCCCCTc                                >  1:344661/1‑150 (MQ=255)
                                                                                                                    ctactaTTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTc                                <  2:1550459/150‑1 (MQ=255)
                                                                                                                    ctactaTTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTc                                >  1:212986/1‑150 (MQ=255)
                                                                                                                    ctactaTTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGATTCGCTACCTCGCCCCTc                                >  1:218360/1‑150 (MQ=255)
                                                                                                                     tactaTTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCa                               >  2:441464/1‑150 (MQ=255)
                                                                                                                      actaTTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCAc                              <  1:1546904/150‑1 (MQ=255)
                                                                                                                       ctaTTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCAct                             >  1:571807/1‑150 (MQ=255)
                                                                                                                            aGCTCACGGTTATAAATCAACATATTGATTTATAAGAATGTAAATCCCCTGAGTGAAACAACGAATTGCTGGGTGTCGCCTTGGCCCAGCTCCCACTATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTc                        <  1:218337/150‑1 (MQ=255)
                                                                                                                             gCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGATTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGAAGGGCTATTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTcc                       >  2:252575/1‑150 (MQ=255)
                                                                                                                               tCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCaa                     >  1:675328/1‑150 (MQ=255)
                                                                                                                                cACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCaat                    >  1:1697981/1‑150 (MQ=255)
                                                                                                                                cACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCaat                    <  1:1181752/150‑1 (MQ=255)
                                                                                                                                 atgcttATAAATCAATATGTTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCaata                   >  1:420844/5‑150 (MQ=255)
                                                                                                                                  cGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTACGGTTCGGTACCTCGCCCCTCACTCGTTGCaataa                  >  1:286440/1‑150 (MQ=255)
                                                                                                                                   gcttATAAATCAATATGTTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCaataat                 >  1:1304492/3‑150 (MQ=255)
                                                                                                                                   gcttATAAATCAATATGTTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGATGTGTGTAGTCTTTGACCATCTCACACGATGGGCTTTTTTTTAACATTTTTCCACTTCGCTACCTCGCCACTCACTCTTTCCaataat                 >  1:81142/3‑150 (MQ=255)
                                                                                                                                    gTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATg                <  1:1357676/150‑1 (MQ=255)
                                                                                                                                        taAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGTTGc            <  1:633545/150‑1 (MQ=255)
                                                                                                                                         aaaTCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGTTGca           <  1:522996/150‑1 (MQ=255)
                                                                                                                                         aaaTCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGTTGca           >  2:1144997/1‑150 (MQ=255)
                                                                                                                                          aaTCAACATATTGATTTATAATCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGTTGcac          >  2:515215/1‑150 (MQ=255)
                                                                                                                                          aaTCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGTTGcac          <  2:773292/150‑1 (MQ=255)
                                                                                                                                          aaTCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGTTGcac          >  1:773073/1‑150 (MQ=255)
                                                                                                                                            tCAACATATTGATTTATAAGCATGGAAATCCCATGAGTGAAACAACGAATTGCTTTGTGTAGTCTATGCCCATCTCCCACGATGGGCTTTTTATTAACATTTTTACGCTTCGCTACCTCGCACCTCACTCTTTCCAATAATGTTTAACaa        >  2:675379/1‑150 (MQ=255)
                                                                                                                                             cAATATGTTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGTTGCACaat       >  2:36142/1‑149 (MQ=255)
                                                                                                                                               aCATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGTTGCACAAGAt     <  2:1250639/150‑1 (MQ=255)
                                                                                                                                               aCATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTGTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGTTGCACAAGAt     >  2:440011/1‑150 (MQ=255)
                                                                                                                                                 atatTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGTTGCACAAGATAc   <  1:1370250/150‑1 (MQ=255)
                                                                                                                                                  tatTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGAGGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGTTGCACAAGATACt  <  2:1095432/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGTTGCACAAGATACT  >  NC_000913/1770243‑1770538

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: