Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,770,390 A→G 5.4% 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 5.4% 340.7 / 7.2 114noncoding (19/106 nt)rprAsRNA antisense activator for RpoS translation, Hfq‑dependent
Reads supporting (aligned to +/- strand):  ref base A (49/57);  new base G (2/4);  total (51/63)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01

Only 100 of 116 total aligned reads displayed.

TGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGTTGCACAAGATACTG  >  NC_000913/1770244‑1770539
                                                                                                                                                  |                                                                                                                                                     
tgctCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCTCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAATATGTTg                                                                                                                                                    <  1:1221112/150‑1 (MQ=255)
tgctCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATAttt                                                                                                                                                    >  1:258944/1‑149 (MQ=255)
     ttCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAATATGTTGATTta                                                                                                                                               <  1:730222/150‑1 (MQ=255)
        cccATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTATTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATaa                                                                                                                                            >  1:1545702/1‑150 (MQ=255)
         ccATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAg                                                                                                                                           <  1:1900246/150‑1 (MQ=255)
          cATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGc                                                                                                                                          >  2:539660/1‑150 (MQ=255)
           atatTTTCCTGCGTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGCTTATAAATCAATATGTTGATTTATaaccg                                                                                                                                         <  1:1015769/150‑4 (MQ=255)
               tttCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGAATGCCTACTATTAGCTCACGGTTATAAATCAATATGTTGATTTATaaccgtgag                                                                                                                                     <  1:895165/150‑8 (MQ=255)
               tttCCTGCTTACCGACATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGa                                                                                                                                     <  1:100708/150‑1 (MQ=255)
                ttCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGaa                                                                                                                                    <  1:1613848/150‑1 (MQ=255)
                  ccTGCTTACCGCCATTCAGCTGGTAGTACATGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTACGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAAt                                                                                                                                  <  2:928749/150‑1 (MQ=255)
                   cTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATc                                                                                                                                 >  1:519695/1‑150 (MQ=255)
                   cTGATTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCGTACTATTACCTCACGGTTATAAATCAAAATATTCATTTATAAGCATGGAAATc                                                                                                                                 >  2:320987/1‑150 (MQ=255)
                       ttACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTATTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAATAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCt                                                                                                                             <  1:1449134/150‑1 (MQ=255)
                        tACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTg                                                                                                                            <  2:846125/150‑1 (MQ=255)
                        tACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTATTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTg                                                                                                                            <  2:1246437/150‑1 (MQ=255)
                          ccGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACTAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAg                                                                                                                          <  2:569334/150‑1 (MQ=255)
                          ccGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAg                                                                                                                          <  1:1331980/150‑1 (MQ=255)
                           cGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGt                                                                                                                         <  2:869354/150‑1 (MQ=255)
                            gCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTg                                                                                                                        >  1:917280/1‑150 (MQ=255)
                                    gCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACg                                                                                                                >  2:709741/1‑150 (MQ=255)
                                    gCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACg                                                                                                                <  2:1798297/150‑1 (MQ=255)
                                     cTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGCGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGa                                                                                                               >  1:668739/1‑150 (MQ=255)
                                     cTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGa                                                                                                               >  1:793935/1‑150 (MQ=255)
                                     cTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGATCACGGTTATAAATCAACATATTTATTTATAAGCATGGAAATCCCCAGAGTGAAACAACGa                                                                                                               >  1:359998/1‑150 (MQ=255)
                                         tagtaCCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTg                                                                                                           <  2:827502/150‑1 (MQ=255)
                                            taCCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCtg                                                                                                        <  2:474221/150‑1 (MQ=255)
                                              ccTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCtgtg                                                                                                      <  2:504239/150‑1 (MQ=255)
                                                 gTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTgtgtgt                                                                                                   >  1:1219517/1‑150 (MQ=255)
                                                     cAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTc                                                                                               >  1:1905171/1‑150 (MQ=255)
                                                          tCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGc                                                                                          >  2:1924153/1‑150 (MQ=255)
                                                           cTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGcc                                                                                         <  1:1134219/150‑1 (MQ=255)
                                                            tttACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTTccc                                                                                        >  2:1350804/1‑150 (MQ=255)
                                                             ttACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCa                                                                                       >  1:1048615/1‑150 (MQ=255)
                                                              tACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCAt                                                                                      >  2:891715/1‑150 (MQ=255)
                                                                 aGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCAtctc                                                                                   >  1:1962941/1‑150 (MQ=255)
                                                                 aGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCCCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCAtctc                                                                                   <  1:1666366/150‑1 (MQ=255)
                                                                  gTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTcc                                                                                  >  1:1063196/1‑150 (MQ=255)
                                                                    ttttAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATTAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCa                                                                                >  1:373790/1‑150 (MQ=255)
                                                                          aCTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATgg                                                                          <  1:1455723/150‑1 (MQ=255)
                                                                             aaTGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCt                                                                       >  1:647514/1‑150 (MQ=255)
                                                                               tGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGGGTAGTCTTTGCCCATCTCCCACGATTGGGttt                                                                     >  2:1282545/1‑150 (MQ=255)
                                                                                 agaCGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCttttt                                                                   <  1:1343022/150‑1 (MQ=255)
                                                                                   aCGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCttttttt                                                                 <  2:261474/150‑1 (MQ=255)
                                                                                    cGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCtttttttt                                                                >  1:1997430/1‑150 (MQ=255)
                                                                                     gAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTa                                                               <  1:1582241/150‑1 (MQ=255)
                                                                                      aaTCTGATGGACGCAAAAAGTCCGTATTCCTACAATTATCTCACGGTTATAAATCAACATATTTATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTaa                                                              <  2:217477/150‑1 (MQ=255)
                                                                                      aaTCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTaa                                                              <  1:819868/150‑1 (MQ=255)
                                                                                      aaTCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTaa                                                              >  1:869900/1‑150 (MQ=255)
                                                                                          tGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACAtt                                                          <  1:690682/150‑1 (MQ=255)
                                                                                          tGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACAtt                                                          >  2:1722912/1‑150 (MQ=255)
                                                                                           gATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGACCATCTCCCACGATGGGCTTTTTTTTAACAttt                                                         >  2:91134/1‑150 (MQ=255)
                                                                                           gATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCATGAGTGAAACAACGAATTGCTGTGTGTATTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACAttt                                                         >  1:996472/1‑150 (MQ=255)
                                                                                           gATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCACTGAGTGAAACAACGAATTGCTGTGTTTAGTCTTTGACCATCTCCCCCGATTGGCTTTTTTTTAACAttt                                                         >  1:1443431/1‑150 (MQ=255)
                                                                                                 cGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGc                                                   >  2:648261/1‑150 (MQ=255)
                                                                                                  gCAAAAAGTCCGTATGCCTACTATTAGCTCAAGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTAGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCt                                                  <  2:1041496/150‑1 (MQ=255)
                                                                                                  gCAAAAAGGCCGTATGCCTACTAATAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCt                                                  <  2:848033/150‑1 (MQ=255)
                                                                                                   cAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCtt                                                 >  2:388633/1‑150 (MQ=255)
                                                                                                    aaaaaGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTc                                                <  1:1514930/150‑1 (MQ=255)
                                                                                                     aaaaGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCg                                               <  2:1545702/150‑1 (MQ=255)
                                                                                                     aaaaGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCg                                               <  2:555635/150‑1 (MQ=255)
                                                                                                      aaaGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGc                                              <  1:945575/150‑1 (MQ=255)
                                                                                                       aaGTCCGTATGCCTACTATTAGCTGACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCt                                             <  1:149681/150‑1 (MQ=255)
                                                                                                           ccGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCt                                         >  2:1892044/1‑150 (MQ=255)
                                                                                                            cGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTc                                        <  1:920781/150‑1 (MQ=255)
                                                                                                            cGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCAGGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTc                                        <  1:1316901/150‑1 (MQ=255)
                                                                                                              tATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGc                                      >  1:132769/1‑150 (MQ=255)
                                                                                                              tATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGc                                      <  2:752611/150‑1 (MQ=255)
                                                                                                              tATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGc                                      <  1:664468/150‑1 (MQ=255)
                                                                                                               aTGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAAAATTTTTCCGCTTCGCTACCTCGcc                                     >  2:1334363/1‑150 (MQ=255)
                                                                                                                tGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCGTTTTTTTAACATTTTTCCGATTCGCTACCTCGccc                                    >  2:1856103/1‑150 (MQ=255)
                                                                                                                ggactactaTTAGCTCACGGGTATAAATCAACATATTGATTTATAAGCATGGACATCCCCTGAGTGATACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGccc                                    <  2:286025/147‑1 (MQ=255)
                                                                                                                 gCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGcccc                                   <  1:1776849/150‑1 (MQ=255)
                                                                                                                  ccTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCt                                  >  1:1715033/1‑150 (MQ=255)
                                                                                                                  ccTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCt                                  <  1:997264/150‑1 (MQ=255)
                                                                                                                    tactaTTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCa                                >  1:1495146/1‑150 (MQ=255)
                                                                                                                    tactaTTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCa                                >  1:1089623/1‑150 (MQ=255)
                                                                                                                     actaTTAGCTGAAGGTTATAAATCAACATTTTGATTTATAATCATGGAAATCCCCTGAGTGTAACAAAGAATTGCTGTGTGTTGTCTTTGCCCATCTCCCACGATTGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCAc                               <  1:278483/150‑1 (MQ=255)
                                                                                                                      ctaTTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCAct                              <  2:631859/150‑1 (MQ=255)
                                                                                                                       taTTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCActc                             <  2:1063196/150‑1 (MQ=255)
                                                                                                                        aTTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTTCCAATATCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACATCGCCCATAActct                            >  1:213488/1‑150 (MQ=255)
                                                                                                                           aGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGCCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCATAAGTCTTTc                         >  2:1385609/1‑150 (MQ=255)
                                                                                                                            gCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTcc                        <  1:701425/150‑1 (MQ=255)
                                                                                                                            gCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTcc                        <  1:1412645/150‑1 (MQ=255)
                                                                                                                               cACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCaat                     <  2:1905171/150‑1 (MQ=255)
                                                                                                                                 cGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCaataa                   <  2:779054/150‑1 (MQ=255)
                                                                                                                                  ggTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCaataat                  >  1:762211/1‑150 (MQ=255)
                                                                                                                                   tttATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGGTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATATTg                 <  2:1589405/149‑1 (MQ=255)
                                                                                                                                   cttATAAATCAATATGTTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGATGTGTGTAGTCTTTGCCAATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCTCTACCTCGCCCCTCACTCTTTCCaataatt                 >  1:40848/2‑149 (MQ=255)
                                                                                                                                     tataAATCAATATGTTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGtt               >  2:1668111/1‑150 (MQ=255)
                                                                                                                                     tataAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGtt               <  2:961390/150‑1 (MQ=255)
                                                                                                                                     tataAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGtt               >  2:972940/1‑150 (MQ=255)
                                                                                                                                     tataAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGtt               >  2:1897158/1‑150 (MQ=255)
                                                                                                                                     tataAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGtt               <  1:1574108/150‑1 (MQ=255)
                                                                                                                                       taAATCAAGATATTGATTTATAAGCATGGAAATCCCCGGAGTGAAACAACTAATTGCTGTGTGTAGTCTTTTCCCATCTCCCACGATGTGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGTTGc             <  1:100221/150‑1 (MQ=255)
                                                                                                                                          aTCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGTTGcaca          >  2:1875261/1‑150 (MQ=255)
                                                                                                                                           tCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGTTGCACaa         <  1:1635349/150‑1 (MQ=255)
                                                                                                                                            cAACATTTTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATGTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGTTGCACAAg        <  2:1744384/150‑1 (MQ=255)
                                                                                                                                                  atTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGTTGCACAAGATACtg  >  2:680896/1‑150 (MQ=255)
                                                                                                                                                  atTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACACTTTCCAATAATGTTGCACAAGATACtg  >  2:1154917/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGTTGCACAAGATACTG  >  NC_000913/1770244‑1770539

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: