Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,770,387 C→T 6.1% noncoding (16/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,3870CT 6.1% 366.4 / 7.8 115noncoding (16/106 nt)rprAsRNA antisense activator for RpoS translation, Hfq‑dependent
Reads supporting (aligned to +/- strand):  ref base C (63/45);  new base T (4/3);  total (67/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 115 total aligned reads displayed.

TTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGTTGCACAAGAT  >  NC_000913/1770239‑1770535
                                                                                                                                                    |                                                                                                                                                    
ttgcttgcTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACa                                                                                                                                                     <  1:1258277/150‑1 (MQ=255)
    ttgctCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATAtt                                                                                                                                                 <  2:207954/150‑1 (MQ=255)
      gctCTTCCACATATTTTCCTGCTTAACGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAATATGTTGa                                                                                                                                               <  1:361464/150‑1 (MQ=255)
       ctctTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGAt                                                                                                                                              >  1:682263/1‑150 (MQ=255)
       ctctTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGAt                                                                                                                                              >  2:1055761/1‑150 (MQ=255)
       ctctTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGAt                                                                                                                                              >  1:1675036/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)
                       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                                                                                                                      >  2:1474715/1‑150 (MQ=255)
                                 gCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTg                                                                                                                    <  1:932459/150‑1 (MQ=255)
                                 gCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTg                                                                                                                    >  1:980267/1‑150 (MQ=255)
                                   cATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGaa                                                                                                                  >  1:504318/1‑150 (MQ=255)
                                    aTTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGaaa                                                                                                                 <  2:816510/150‑1 (MQ=255)
                                     ttCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAaac                                                                                                                >  1:1250639/1‑150 (MQ=255)
                                     ttCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAaac                                                                                                                <  1:313906/150‑1 (MQ=255)
                                      tCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAaaca                                                                                                               >  1:415425/1‑150 (MQ=255)
                                         gCTTGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACg                                                                                                            <  1:757894/150‑1 (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)
                                                                     cAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCAtct                                                                                >  1:600057/1‑150 (MQ=255)
                                                                       gTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTcc                                                                              >  1:1254436/1‑150 (MQ=255)
                                                                         ttttAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCa                                                                            >  2:1357676/1‑150 (MQ=255)
                                                                         ttttAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCa                                                                            >  1:425274/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)
                                                                                         cGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCtttttttt                                                            >  1:1468096/1‑150 (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)
                                                                                                      cGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTCTAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGc                                               <  2:1577254/150‑1 (MQ=255)
                                                                                                      cGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGc                                               >  2:365980/1‑150 (MQ=255)
                                                                                                       gCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCt                                              >  1:1073691/1‑150 (MQ=255)
                                                                                                       gCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCt                                              >  2:124745/1‑150 (MQ=255)
                                                                                                         aaaaaGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATCGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCAACGATGGGCTTTTTTTTAACATTTTTCCGCtta                                            >  2:571736/1‑149 (MQ=255)
                                                                                                          aaaaGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCg                                           >  2:1331958/1‑150 (MQ=255)
                                                                                                          aaaaGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCg                                           <  1:1332046/150‑1 (MQ=255)
                                                                                                          aaaaGTCCGTATGCCTAATATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCg                                           <  1:534350/150‑1 (MQ=255)
                                                                                                           aaaGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGc                                          >  1:605705/1‑150 (MQ=255)
                                                                                                            aaGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCt                                         <  1:1177305/150‑1 (MQ=255)
                                                                                                            aaGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCt                                         <  1:1000828/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                             >  1:212986/1‑150 (MQ=255)
                                                                                                                        ctactaTTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTc                             <  2:1550459/150‑1 (MQ=255)
                                                                                                                        ctactaTTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGATTCGCTACCTCGCCCCTc                             >  1:218360/1‑150 (MQ=255)
                                                                                                                         tactaTTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCa                            >  2:441464/1‑150 (MQ=255)
                                                                                                                          actaTTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCAc                           <  1:1546904/150‑1 (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)
                                                                                                                                      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       >  1:773073/1‑150 (MQ=255)
                                                                                                                                              aaTCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGTTGcac       <  2:773292/150‑1 (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)
                                                                                                                                                    |                                                                                                                                                    
TTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGTTGCACAAGAT  >  NC_000913/1770239‑1770535

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: