Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,188,710 T→C 5.3% L1121P (CTG→CCG)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,188,7100TC 5.3% 189.1 / 2.0 77L1121P (CTG→CCG) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base T (27/43);  new base C (2/2);  total (29/48)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

AAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAA  >  NC_000913/4188563‑4188853
                                                                                                                                                   |                                                                                                                                               
aaGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTg                                                                                                                                                <  1:269560/149‑1 (MQ=255)
aaGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTg                                                                                                                                                <  1:91280/149‑1 (MQ=255)
aaGATCTGCGGCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTg                                                                                                                                                <  1:556474/149‑1 (MQ=255)
  gATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGgc                                                                                                                                              >  1:53570/1‑149 (MQ=255)
    ccTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGgcgc                                                                                                                                            <  1:412133/148‑1 (MQ=255)
      tGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTACCGACGATCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAACGCGATTATTCAGCTGGAAGATGGCGTGCAGATCAGCGCTGGTGACACCCTGGCGCGt                                                                                                                                          >  2:519203/1‑149 (MQ=255)
       gcgcCCGGCACCTAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTa                                                                                                                                         <  1:169712/149‑1 (MQ=255)
       gCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTa                                                                                                                                         <  2:480713/149‑1 (MQ=255)
       gCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTa                                                                                                                                         <  2:678964/149‑1 (MQ=255)
        cGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTAt                                                                                                                                        >  2:111320/1‑149 (MQ=255)
            cGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCTGGTACCGATATTCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCGCTGGTGACACCCTGGCGCGTATTCCg                                                                                                                                    <  1:61185/149‑1 (MQ=255)
              gCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCa                                                                                                                                  <  2:309157/149‑1 (MQ=255)
               taCTGAAAACGGTTGAGGCCTAGGGGCACGACGTTCTGACACCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGAGTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAg                                                                                                                                 <  2:554066/148‑1 (MQ=255)
                  tGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGaa                                                                                                                              >  2:369135/1‑149 (MQ=255)
                  tGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGAGGACACCCTGGCGCGTATTCCGCAGGaa                                                                                                                              >  1:497338/1‑149 (MQ=255)
                           ttGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCGGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGt                                                                                                                     >  2:348476/1‑149 (MQ=255)
                                 tccAGGGTAACGACGTGCTGGTCGCCGGTACCGATATGCCAGCGCAGTACGTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAg                                                                                                               <  2:746844/147‑1 (MQ=255)
                                 cgcAGGTTACCGACGCTCTGATCCCAGGTCGCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTCTTCAGCTGGAAGATGGCGTACAGATCAGCCCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAg                                                                                                               <  2:620221/147‑1 (MQ=255)
                                 cTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCCGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAg                                                                                                               >  2:696630/1‑149 (MQ=255)
                                 cTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCCGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAg                                                                                                               >  2:50664/1‑149 (MQ=255)
                                 cTCAGGGTAACGACGGTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAg                                                                                                               <  1:777482/149‑1 (MQ=255)
                                  tCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAgg                                                                                                              >  1:463707/1‑149 (MQ=255)
                                   cAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGa                                                                                                             >  1:11733/1‑149 (MQ=255)
                                   cAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGa                                                                                                             >  1:393773/1‑149 (MQ=255)
                                            gacgTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCgg                                                                                                    <  2:296987/149‑1 (MQ=255)
                                              cgTGCTGTTCCCACGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCCGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCggtg                                                                                                  <  1:5449/149‑1 (MQ=255)
                                                ttCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCggtggt                                                                                                <  2:498160/149‑1 (MQ=255)
                                                      tCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGCCATCACCGGTGGGCTGCcg                                                                                          >  1:335622/1‑149 (MQ=255)
                                                              aCCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGCGGTGGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGc                                                                                  <  1:607643/149‑1 (MQ=255)
                                                               ccGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCGGCGCGTTGCg                                                                                 >  1:272747/1‑149 (MQ=255)
                                                                cGATATGCCAGCGCAGTACTTACTGCCGGGTAAAGCGATTGTTAAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACGAAGGACATCACCGGTGGTCTGCGCGGCGGTGCgg                                                                                >  2:202486/1‑149 (MQ=255)
                                                                 ggtatGCGAGCGCAGTCCTTCCTGCCGGGTAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGTCACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGAc                                                                              <  2:43507/147‑1 (MQ=255)
                                                                    atGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCt                                                                            <  1:495286/149‑1 (MQ=255)
                                                                         aGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCg                                                                       >  2:474518/1‑149 (MQ=255)
                                                                         aGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCg                                                                       >  2:65919/1‑149 (MQ=255)
                                                                          gcgcAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGa                                                                      >  2:506827/1‑149 (MQ=255)
                                                                           cgcAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGACGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGaa                                                                     <  2:683177/149‑1 (MQ=255)
                                                                             cAGTACTTCCGGCCGGGTAAAGCGATTGTTCACCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCGGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACAGCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGc                                                                   <  2:750014/149‑1 (MQ=255)
                                                                               gTACTCCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCAc                                                                 <  1:81726/149‑1 (MQ=255)
                                                                                  cTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGAACTGTTCGAAGCAcgtc                                                              >  1:750102/1‑149 (MQ=255)
                                                                                  cTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGAGCGTTGCGGACCGGTTCGAGGCAagtc                                                              >  2:588072/1‑149 (MQ=255)
                                                                                  cTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGGGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGGCTGCCGCGCGTTGCGGACCTGTTCGAAGCAcgtc                                                              >  1:750107/1‑149 (MQ=255)
                                                                                      cTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTcc                                                          >  2:527678/1‑149 (MQ=255)
                                                                                         ccGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGaa                                                       <  1:783067/149‑1 (MQ=255)
                                                                                         ccGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGaa                                                       <  1:588692/149‑1 (MQ=255)
                                                                                         ccGGGTAAAGCGAGTGTTCGGCTGGAAGACGGCTGACGGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGaa                                                       <  2:624723/149‑1 (MQ=255)
                                                                                                aaGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCg                                                <  1:690117/149‑1 (MQ=255)
                                                                                                  gCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGc                                              >  1:501450/1‑149 (MQ=255)
                                                                                                  gCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGc                                              >  1:192596/1‑149 (MQ=255)
                                                                                                        gTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACCGCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCt                                        <  2:484682/149‑1 (MQ=255)
                                                                                                        gTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCt                                        <  1:195479/149‑1 (MQ=255)
                                                                                                            ggCTGGAAGACGGTGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCt                                    <  1:489241/148‑1 (MQ=255)
                                                                                                             gcgggaggaGGGCGTACCGATCAGCTCGGGTGCCCCCATGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTg                                   <  2:468504/146‑1 (MQ=255)
                                                                                                             gCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTg                                   <  1:432316/149‑1 (MQ=255)
                                                                                                              cTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGa                                  <  1:369135/149‑1 (MQ=255)
                                                                                                                 gAAGATGGCGTACAGATCAGCTCTGGTGACACCCCGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAAt                               <  1:50664/149‑1 (MQ=255)
                                                                                                                 gAAGATGGCGTACAGATCAGCTCTGGTGACACCCCGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTGCGATTCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAAt                               <  1:696630/149‑1 (MQ=255)
                                                                                                                  aaGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATc                              <  1:637996/149‑1 (MQ=255)
                                                                                                                      tGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCCATCCTGGCTGAAATCAGca                          >  2:91082/1‑148 (MQ=255)
                                                                                                                      tGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCg                          >  2:467684/1‑149 (MQ=255)
                                                                                                                           tACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATc                     >  1:703955/1‑149 (MQ=255)
                                                                                                                                aacAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTc                >  2:559733/3‑149 (MQ=255)
                                                                                                                                aTCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTc                <  1:155643/149‑1 (MQ=255)
                                                                                                                                aTCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTc                <  1:533188/149‑1 (MQ=255)
                                                                                                                                   cgCTCTGGGGACACCCTGGCGCGTGTTCCGCAGGAATCCGGCGGTACCGAGGACCTCCCCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGTAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCtt             <  2:63208/148‑1 (MQ=255)
                                                                                                                                   agcccgggTGACACCCGGGCCCGTATTCCGCAGGAATCCGGCGGTACCAACGACATCACCGGTGGGCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCtt             <  2:185755/143‑1 (MQ=255)
                                                                                                                                   aGCTGTGGTGGCACCCTGGCGCGTATTCCGGAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCtt             <  1:242996/149‑1 (MQ=255)
                                                                                                                                   aGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCtt             <  1:293884/149‑1 (MQ=255)
                                                                                                                                     gtctGGCGGCGCCCTGGCGGCTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGGGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCg           <  2:497338/148‑1 (MQ=255)
                                                                                                                                     ctctGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCTAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCg           <  1:196970/149‑1 (MQ=255)
                                                                                                                                     ctctGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGCAATCAGCGGTATCGTTTCCTTCg           >  1:659880/1‑149 (MQ=255)
                                                                                                                                     ctctGGTGACACCCGGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACAGCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCg           <  1:81/149‑1 (MQ=255)
                                                                                                                                     ctcgGGTGCCACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCg           <  1:668209/149‑1 (MQ=255)
                                                                                                                                      tctGGTCGCACCCTGGCGCGTATTCCGGAGGAATCCCGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCgg          <  2:709674/149‑1 (MQ=255)
                                                                                                                                           tGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAg     >  1:43451/1‑149 (MQ=255)
                                                                                                                                              cccccTGGCGGGTATCCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGaaa  <  2:750102/147‑1 (MQ=255)
                                                                                                                                              cacCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGGCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGaaa  <  2:750107/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
AAGATCTGCGTCCGGCACTGAAAATCGTTGATGCTCAGGGTAACGACGTTCTGATCCCAGGTACCGATATGCCAGCGCAGTACTTCCTGCCGGGTAAAGCGATTGTTCAGCTGGAAGATGGCGTACAGATCAGCTCTGGTGACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAA  >  NC_000913/4188563‑4188853

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: