Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,187,597 C→T 92.0% P750L (CCG→CTG)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,187,5970CT92.0% 351.3 / 14.3 130P750L (CCG→CTG) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (5/6);  new base T (53/64);  total (59/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

Only 100 of 130 total aligned reads displayed.

CAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGT  >  NC_000913/4187452‑4187742
                                                                                                                                                 |                                                                                                                                                 
cAAACGGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGa                                                                                                                                                <  1:992111/149‑1 (MQ=255)
 aaaCTGAAACCGTGATTAACCGCGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGATTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGAt                                                                                                                                               <  2:783747/149‑1 (MQ=255)
   aCTGAAGCCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGCTGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTCGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATgg                                                                                                                                             <  2:811984/149‑1 (MQ=255)
   aCTGAAACCGTGATTAACCGTGAAGGTCGGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGTTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATgg                                                                                                                                             >  1:951816/1‑149 (MQ=255)
    ctcaacccGTGATTAACCGTGACGGGCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGACTCGTCAGCGTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGc                                                                                                                                            <  1:926848/143‑1 (MQ=255)
      gAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTc                                                                                                                                          <  1:1306761/149‑1 (MQ=255)
            gTGATTACCCGGGACCGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCcatcat                                                                                                                                    <  1:1269215/149‑1 (MQ=255)
            gTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCcatcat                                                                                                                                    <  2:598698/149‑1 (MQ=255)
            gTGATTAACCGTGACGGGCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTCCATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCcatcat                                                                                                                                    <  2:761003/149‑1 (MQ=255)
             tGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCatcatc                                                                                                                                   <  2:775301/149‑1 (MQ=255)
             tGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCatcatc                                                                                                                                   <  2:290168/149‑1 (MQ=255)
             tGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCatcatc                                                                                                                                   <  2:1426369/149‑1 (MQ=255)
             tGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCatcatc                                                                                                                                   >  1:687346/1‑149 (MQ=255)
              gATTAACCGTGACGGTCAGGACGGGAATCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCg                                                                                                                                  <  1:1051841/149‑1 (MQ=255)
              gATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCg                                                                                                                                  <  2:582494/149‑1 (MQ=255)
                  acccGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAAc                                                                                                                              <  1:86112/147‑1 (MQ=255)
                  aaCCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATCGCCGAATCCTGTGCTCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAAc                                                                                                                              <  1:1427671/149‑1 (MQ=255)
                   aCCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTGGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCGTCGGAACg                                                                                                                             >  2:1093121/1‑149 (MQ=255)
                     cGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGcc                                                                                                                           >  2:290678/1‑149 (MQ=255)
                      gTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCa                                                                                                                          <  2:1343746/149‑1 (MQ=255)
                             tCAGGAAGGGAAGCAGGTTTCCTTCAACAGCATCTGCATGATGGCCGGCTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACcg                                                                                                                   <  2:246750/149‑1 (MQ=255)
                             tCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGCATGCGCGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAaccaccc                                                                                                                   >  2:559500/1‑148 (MQ=255)
                                         gCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAAGGCCAATCACCGCGAACTTCCGTg                                                                                                       >  1:409697/1‑149 (MQ=255)
                                          cAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGa                                                                                                      <  1:744886/149‑1 (MQ=255)
                                            ggTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAg                                                                                                    >  2:1342682/1‑149 (MQ=255)
                                            ggTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAg                                                                                                    <  2:158928/149‑1 (MQ=255)
                                            ggTTTCCTTCAACAGCATCTACATGATGGCCGAATCCGGTGCGCGTGGGTCTCCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAg                                                                                                    <  2:1390967/149‑1 (MQ=255)
                                            ggTTTCCTTCAACAGCATCTACAGGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAg                                                                                                    <  2:424901/149‑1 (MQ=255)
                                            ggTTTCCTTCAACAGCAGCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAg                                                                                                    <  2:343378/149‑1 (MQ=255)
                                             gTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAgg                                                                                                   >  1:3997/1‑149 (MQ=255)
                                                   ttCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGCTTCGTCCGCTTGCTGGCATGCGAGGTCTGATGGCGAAGCTGGATGGCTCCATCACCGGAACGCCAATCGCCGGGACCTTCCGTGAGGGTCGGaa                                                                                             >  2:1004763/1‑149 (MQ=255)
                                                   ttCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGGGCAGGTCTGaa                                                                                             >  2:101615/1‑149 (MQ=255)
                                                     cAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGGCAATCACCGCGCACTTCCGTGAAGGTCTGAACg                                                                                           >  1:1100560/1‑149 (MQ=255)
                                                      aaCAGCATCTACATGATGGCCGACTGCGCTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGt                                                                                          >  1:816969/1‑149 (MQ=255)
                                                      aaCAGCATCTACATGATAGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGAGTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAGCGt                                                                                          >  2:864532/1‑149 (MQ=255)
                                                       aCAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTGTGGACGTa                                                                                         >  2:12634/1‑149 (MQ=255)
                                                       aCAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTa                                                                                         <  1:724850/149‑1 (MQ=255)
                                                          gCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTc                                                                                      <  1:659178/149‑1 (MQ=255)
                                                          gCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTc                                                                                      <  2:687346/149‑1 (MQ=255)
                                                              cTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGt                                                                                  >  1:37807/1‑149 (MQ=255)
                                                              cTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCCAACTTCCGTGAAGCTCTGCACGTACTCCAGt                                                                                  >  1:737153/1‑149 (MQ=255)
                                                               tACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTa                                                                                 >  1:1294296/1‑149 (MQ=255)
                                                               tACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCACAATATTAAGCTTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTa                                                                                 <  1:1060254/149‑1 (MQ=255)
                                                               tACATGATGGCCGACTCCGGTGCGCGTCGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTcaaata                                                                                 >  2:361434/1‑144 (MQ=255)
                                                                   tgatgGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTc                                                                             >  1:149357/1‑149 (MQ=255)
                                                                   tgatgGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGGGAAGGTCTGAACGTACTACAGTACTTc                                                                             >  2:120407/1‑149 (MQ=255)
                                                                    gatgGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCa                                                                            >  2:812724/1‑149 (MQ=255)
                                                                     atgGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCAt                                                                           >  1:1359578/1‑149 (MQ=255)
                                                                      tgGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCAtc                                                                          <  1:437402/149‑1 (MQ=255)
                                                                      tgGCCGACTCCGCTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCAtc                                                                          <  1:188665/149‑1 (MQ=255)
                                                                      tgGCCGACTCAGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCAtc                                                                          <  1:188441/149‑1 (MQ=255)
                                                                          cGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTcca                                                                      >  1:307155/1‑149 (MQ=255)
                                                                           gACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACAGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATcccccc                                                                     >  1:633607/1‑147 (MQ=255)
                                                                              tCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTccaccca                                                                  >  1:675905/1‑149 (MQ=255)
                                                                                 ggTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTTCAGTACTTCGGCTCTACCCACgg                                                               >  1:836299/1‑149 (MQ=255)
                                                                                      gctgggTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTc                                                          <  2:88807/145‑1 (MQ=255)
                                                                                       cgTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCt                                                                                                                                           <  2:1206456/67‑1 (MQ=255)
                                                                                       cgTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCt                                                                                                                                           >  1:1206456/1‑67 (MQ=255)
                                                                                            ttCTGCGGCACAGATTCGTCGGCGTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAACCACCGCTAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAg                                                    <  2:203276/149‑1 (MQ=255)
                                                                                            ttCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCTAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAg                                                    <  1:290496/149‑1 (MQ=255)
                                                                                              cTGCGGCACAGATGCGTCAGCTTGGTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGt                                                  <  2:672366/149‑1 (MQ=255)
                                                                                                gCGGCACAGATTCGTCAGCTTGCTGGTAGGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGGACGTACTCCAGTACTTCATCTCCACCCACGGTGCCCGTAAAGGTCt                                                >  2:13622/1‑149 (MQ=255)
                                                                                                 cGGCACCGAGTCGGCAACTTGCGGGTATGCGGGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTg                                               <  2:658989/149‑1 (MQ=255)
                                                                                                 cGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGTTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTg                                               <  2:120396/149‑1 (MQ=255)
                                                                                                 cGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCCTCTCCACCCACGGTGCTCGTAAAGGTCTg                                               >  1:395917/1‑149 (MQ=255)
                                                                                                   gCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTGAAGGTCTGGc                                             >  2:471173/1‑149 (MQ=255)
                                                                                                   gCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGc                                             <  1:1178274/149‑1 (MQ=255)
                                                                                                   gCACAGATTCGGTAGTTTGCTGGAAGGCGTGGTCTGATGGCGAAGCTTGATGTCTCCATCAACGAAACGCCCATCACCGCCAACTTCCGTGAAGGTGTGAACGTACTCCAGTACTTTATCTCCACCCACGGTGCTTGTAAAGGTCTGGc                                             <  2:1143099/149‑1 (MQ=255)
                                                                                                      caGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAGCGCCAATAACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCCATACTTCATCTCCACCCACGCAGCTCTTGAAGGTCTGGCGGa                                          >  1:525351/1‑149 (MQ=255)
                                                                                                       aGATGCGCCAGCTTGCTTGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGAt                                         <  1:1157429/149‑1 (MQ=255)
                                                                                                         aTTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTCCCGTGAAGGTCTGAACGTCCTCCAGTACTTCATCTCCACCCCCGGTGCTCGTAAAGATCTGGCGGATAc                                       >  1:414670/1‑149 (MQ=255)
                                                                                                           tcgtcgGCTGGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAACCACCGCGAACTTCCGTGACGGTCTGAACGTACTCCAGGACTTCATCTCCCCCCACGGTGCTCGTAAAGGTCTGGCGGATACCg                                     <  2:1008252/149‑1 (MQ=255)
                                                                                                           tCGGCAGCTTGCTGGTATGCGTTGCCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCg                                     <  2:272690/149‑1 (MQ=255)
                                                                                                             gTCATCGTGCTGGGAGGCGTGGTCTGCTGGCGAATCTTGATTGCTTCATTAACTGTACGCTAATCACCGCCAACTTCCGTGAAGGTCTTAACGTACTCCAGTACTTCATCGCCATCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCa                                   <  2:489049/149‑1 (MQ=255)
                                                                                                                aGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTg                                >  1:475375/1‑149 (MQ=255)
                                                                                                                 gCTTGCTGGTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGa                               >  1:542907/1‑149 (MQ=255)
                                                                                                                 gCTTGCTGGGATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGa                               <  1:566273/149‑1 (MQ=255)
                                                                                                                        ggTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGc                        <  2:262574/149‑1 (MQ=255)
                                                                                                                        ggTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGc                        <  1:853699/149‑1 (MQ=255)
                                                                                                                         gTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCg                       >  1:30914/1‑149 (MQ=255)
                                                                                                                         gTATGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCg                       >  1:1221643/1‑149 (MQ=255)
                                                                                                                          tATGCGTGGGCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAACCTGCGa                      >  2:857146/1‑149 (MQ=255)
                                                                                                                           aTGCGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGaa                     <  2:1294296/149‑1 (MQ=255)
                                                                                                                             gagTGGTCGGACGGCGAAGCTGGATGGCTCCATCACCGACACGCCCGTAACCGTGACCTTCCGTGAAGGTCTGAACGTACTCCAGTCCTTCATCTCCTCCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACt                   <  2:417301/147‑1 (MQ=255)
                                                                                                                              cGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGCACTc                  >  2:259147/1‑149 (MQ=255)
                                                                                                                              cGTGGTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTc                  >  2:354985/1‑149 (MQ=255)
                                                                                                                                  gTCTGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCAGGCGGATACCGCACTGAAAACTGCGAACACCGGt              >  1:455076/1‑149 (MQ=255)
                                                                                                                                     tgctGGCGAGGGTGGATGGCTCCATCATCGAAACGCCAACCACCGCGAACTTCCGGGAAGGTCTTAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTAc           <  2:18601/146‑1 (MQ=255)
                                                                                                                                     tGATGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTAc           <  1:794833/149‑1 (MQ=255)
                                                                                                                                      gATGGCGAAGCTGGATGGCTCCATCATCGAAACGCTAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGTGGATACCGCTCTGAAAACTGCGAACTCCGGTTAcc          <  1:812724/149‑1 (MQ=255)
                                                                                                                                       aTGGCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCt         <  1:15360/149‑1 (MQ=255)
                                                                                                                                       aTGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAGGGTCTGGCGGATACCGCACTGAAAACTGCGGACTCCGAGTACCt         >  2:71717/1‑149 (MQ=255)
                                                                                                                                        tGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTg        >  1:872080/1‑149 (MQ=255)
                                                                                                                                         ggCGAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTCCCTGa       >  2:1167892/1‑149 (MQ=255)
                                                                                                                                         ggCCAAGCCGGATGGCTCCAGCACCGAACCGCCAATCACCGCGATCTTCCGTGCAGGTCTGAACGTACTCCAGTACTGCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGa       <  2:29859/149‑1 (MQ=255)
                                                                                                                                          gCGAAGCTGGCTGGCTTCATCATCGACACGTCAATCCCCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGAc      <  1:1090360/149‑1 (MQ=255)
                                                                                                                                            gAAGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCGGACtc    >  1:1045803/1‑149 (MQ=255)
                                                                                                                                             aaGCAGGATGACTCCAACATCGAAACCCCAATCACCGCAAACTACCCTGATGCTCCGAACGTACTCCAGTACTTCACCTCCCCCCACAGTGCCCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACtcg   <  2:836299/149‑1 (MQ=255)
                                                                                                                                              aGCTGGATGGCTCCCTCATCGACCCGCCAATCCCCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACtcgt  <  1:864532/149‑1 (MQ=255)
                                                                                                                                              aGCTGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACtcgt  >  1:1367126/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGT  >  NC_000913/4187452‑4187742

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: