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

Only 100 of 119 total aligned reads displayed.

TTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGTTGCACAAG  >  NC_000913/1770238‑1770533
                                                                                                                                                     |                                                                                                                                                  
tttGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAAc                                                                                                                                                    <  2:239820/150‑1 (MQ=255)
tttGCTTGCTCATCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTGGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCTCGGTTATAAATCAAc                                                                                                                                                    <  2:1322975/150‑1 (MQ=255)
 ttgcttgcTCTTCCCCTTATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACa                                                                                                                                                   <  1:1900387/150‑1 (MQ=255)
      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)
                   atTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATg                                                                                                                                 <  1:504410/150‑1 (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)
                            cTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATcccc                                                                                                                        >  2:701425/1‑150 (MQ=255)
                             ttACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTATTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAATAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCt                                                                                                                       <  1:1449134/150‑1 (MQ=255)
                              tACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTg                                                                                                                      <  2:846125/150‑1 (MQ=255)
                              tACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGACCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTg                                                                                                                      <  1:1511739/150‑1 (MQ=255)
                              tACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTATTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTg                                                                                                                      <  2:1246437/150‑1 (MQ=255)
                                ccGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACTAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAg                                                                                                                    <  2:569334/150‑1 (MQ=255)
                                  gCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTg                                                                                                                  >  1:917280/1‑150 (MQ=255)
                                    cATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGaa                                                                                                                <  1:899373/150‑1 (MQ=255)
                                          gCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACg                                                                                                          >  2:709741/1‑150 (MQ=255)
                                             ggTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAAt                                                                                                       >  1:1177553/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)
                                                         cGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAg                                                                                           <  2:1942863/150‑1 (MQ=255)
                                                           cAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTc                                                                                         >  1:1905171/1‑150 (MQ=255)
                                                                tCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCTTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGc                                                                                    <  2:51877/150‑1 (MQ=255)
                                                                tCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGc                                                                                    >  2:1924153/1‑150 (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:1636998/1‑150 (MQ=255)
                                                                                   aaTGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCt                                                                 >  1:647514/1‑150 (MQ=255)
                                                                                     tGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGGGTAGTCTTTGCCCATCTCCCACGATTGGGttt                                                               >  2:1282545/1‑150 (MQ=255)
                                                                                       agaCGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCttttt                                                             >  1:2001596/1‑150 (MQ=255)
                                                                                       agaCGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCttttt                                                             <  1:1343022/150‑1 (MQ=255)
                                                                                         aCGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCttttttt                                                           <  2:261474/150‑1 (MQ=255)
                                                                                          cGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTTAAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGTCCATCTCCCACGAAGGGCtttttttt                                                          >  1:721446/1‑150 (MQ=255)
                                                                                          cGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCtttttttt                                                          >  1:1997430/1‑150 (MQ=255)
                                                                                           gAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTa                                                         <  2:1805773/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)
                                                                                                            aaaGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGc                                        <  1:945575/150‑1 (MQ=255)
                                                                                                             aaGTCCGTATGCCTACTATTAGCTGACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCt                                       <  1:149681/150‑1 (MQ=255)
                                                                                                                 ccGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCt                                   >  2:1892044/1‑150 (MQ=255)
                                                                                                                  cGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCAGGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTc                                  <  1:1316901/150‑1 (MQ=255)
                                                                                                                    tATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGc                                <  1:664468/150‑1 (MQ=255)
                                                                                                                    tATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGc                                >  1:132769/1‑150 (MQ=255)
                                                                                                                    tATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGCAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGc                                <  1:1381825/150‑1 (MQ=255)
                                                                                                                     aTGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAAAATTTTTCCGCTTCGCTACCTCGcc                               >  2:1334363/1‑150 (MQ=255)
                                                                                                                      tGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCGTTTTTTTAACATTTTTCCGATTCGCTACCTCGccc                              >  2:1856103/1‑150 (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:1089623/1‑150 (MQ=255)
                                                                                                                          tactaTTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCa                          >  1:1495146/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)
                                                                                                                                       cGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCaataa             <  2:779054/150‑1 (MQ=255)
                                                                                                                                        ggTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCaataat            >  1:762211/1‑150 (MQ=255)
                                                                                                                                        ggTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCaataat            <  2:860403/150‑1 (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:972940/1‑150 (MQ=255)
                                                                                                                                           tataAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGtt         <  2:961390/150‑1 (MQ=255)
                                                                                                                                           tataAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGtt         <  1:1574108/150‑1 (MQ=255)
                                                                                                                                           tataAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGtt         >  2:1897158/1‑150 (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)
                                                                                                                                                     |                                                                                                                                                  
TTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGTACCTGTCGCAAATTCTTTACAGTTTTTAAACTAATGAGACGAATCTGATCGACGCAAAAAGTCCGTATGCCTACTATTAGCTCACGGTTATAAATCAACATATTGATTTATAAGCATGGAAATCCCCTGAGTGAAACAACGAATTGCTGTGTGTAGTCTTTGCCCATCTCCCACGATGGGCTTTTTTTTAACATTTTTCCGCTTCGCTACCTCGCCCCTCACTCTTTCCAATAATGTTGCACAAG  >  NC_000913/1770238‑1770533

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: