Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,186,706 G→A V453I (GTT→ATT)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,186,7060GA98.6% 279.8 / ‑4.3 73V453I (GTT→ATT) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (44/28);  total (45/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

CGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCG  >  NC_000913/4186558‑4186850
                                                                                                                                                    |                                                                                                                                                
cgAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGAtt                                                                                                                                                <  2:920873/151‑1 (MQ=255)
cgAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGAtt                                                                                                                                                >  1:440466/1‑151 (MQ=255)
 gaagaaGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGAttt                                                                                                                                               <  2:642008/151‑1 (MQ=255)
 gaagaaGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGAttt                                                                                                                                               >  2:1134275/1‑151 (MQ=255)
   agaagCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTtgt                                                                                                                                             <  2:1046044/151‑1 (MQ=255)
      agCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCg                                                                                                                                          <  2:413072/151‑1 (MQ=255)
        cTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGc                                                                                                                                        >  1:17305/1‑151 (MQ=255)
        cTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGc                                                                                                                                        >  1:194937/1‑151 (MQ=255)
        cTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGc                                                                                                                                        >  2:464239/1‑151 (MQ=255)
          gTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGc                                                                                                                                        <  2:624943/149‑1 (MQ=255)
          gTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCat                                                                                                                                      >  1:232255/1‑151 (MQ=255)
           tCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCat                                                                                                                                      >  2:82803/1‑150 (MQ=255)
            cGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCata                                                                                                                                     <  2:1120429/150‑1 (MQ=255)
             gTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCatat                                                                                                                                    >  2:403756/1‑150 (MQ=255)
             gTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCatat                                                                                                                                    >  2:928046/1‑150 (MQ=255)
               ttGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAAc                                                                                                                                 >  2:865486/1‑151 (MQ=255)
                  ggATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGcc                                                                                                                              <  1:510981/151‑1 (MQ=255)
                       tCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACtt                                                                                                                         >  2:157336/1‑151 (MQ=255)
                           ggACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGAt                                                                                                                     >  2:463471/1‑151 (MQ=255)
                           ggACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGAt                                                                                                                     <  2:866507/151‑1 (MQ=255)
                                  gTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGAcc                                                                                                              >  1:326070/1‑151 (MQ=255)
                                     aTCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGAcc                                                                                                              >  1:850468/1‑148 (MQ=255)
                                     aTCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAg                                                                                                            <  1:430707/150‑1 (MQ=255)
                                     aTCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGa                                                                                                           <  2:786198/151‑1 (MQ=255)
                                       ccGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGAt                                                                                                          <  1:554143/150‑1 (MQ=255)
                                       ccGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATg                                                                                                         >  1:473775/1‑151 (MQ=255)
                                         gcgAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGc                                                                                                       >  1:680260/1‑151 (MQ=255)
                                         gcgAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGc                                                                                                       >  2:322343/1‑151 (MQ=255)
                                          cgAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCt                                                                                                      <  2:1162760/151‑1 (MQ=255)
                                               acCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCCCCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGtt                                                                                                   <  1:133503/149‑1 (MQ=255)
                                               acCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTc                                                                                                  >  1:446064/1‑150 (MQ=255)
                                                      aCTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCg                                                                                          >  2:586026/1‑151 (MQ=255)
                                                      aCTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGCGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTAc                                                                                            >  2:118507/1‑149 (MQ=255)
                                                            gAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACg                                                                                    <  1:84219/151‑1 (MQ=255)
                                                              aCCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCt                                                                                  >  2:228333/1‑151 (MQ=255)
                                                                           tctGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTg                                                                     >  2:575566/1‑151 (MQ=255)
                                                                            ctGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTgg                                                                    >  1:805347/1‑151 (MQ=255)
                                                                              gCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGa                                                                   >  2:969697/1‑150 (MQ=255)
                                                                                 ccGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAgc                                                                <  1:575999/150‑1 (MQ=255)
                                                                                  cGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAgcgc                                                              <  1:340154/151‑1 (MQ=255)
                                                                                   gTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGcgcg                                                             >  2:467764/1‑151 (MQ=255)
                                                                                         gTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCt                                                       <  2:1161511/151‑1 (MQ=255)
                                                                                         gTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCt                                                       <  2:798123/151‑1 (MQ=255)
                                                                                                ggCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGatgatg                                                >  1:1042494/1‑151 (MQ=255)
                                                                                                 gCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGt                                               >  1:38966/1‑151 (MQ=255)
                                                                                                     ttGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGGGCGTGCGCTGATGATGTCTa                                            <  2:1079459/150‑1 (MQ=255)
                                                                                                     ttGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTa                                            <  2:390627/150‑1 (MQ=255)
                                                                                                       gAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACca                                         <  2:159902/151‑1 (MQ=255)
                                                                                                        aaCCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACca                                         >  2:682434/1‑150 (MQ=255)
                                                                                                         aCCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACcaac                                       >  2:508900/1‑151 (MQ=255)
                                                                                                          ccGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACcaaca                                      <  2:1050597/151‑1 (MQ=255)
                                                                                                             gTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAacaaca                                   <  2:763954/151‑1 (MQ=255)
                                                                                                                cTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATcc                                >  1:163363/1‑151 (MQ=255)
                                                                                                                   aTCGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGt                             >  2:603062/1‑151 (MQ=255)
                                                                                                                     cGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTcc                           <  1:463454/151‑1 (MQ=255)
                                                                                                                     cGAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTcc                           >  1:924082/1‑151 (MQ=255)
                                                                                                                      gAAGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTccc                          >  1:1024512/1‑151 (MQ=255)
                                                                                                                        aGGTAAAGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCcg                        >  2:24485/1‑151 (MQ=255)
                                                                                                                            aaaGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCcggcg                     >  2:434220/1‑150 (MQ=255)
                                                                                                                              aGCTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCcggcgaac                  >  1:910970/1‑151 (MQ=255)
                                                                                                                                cTATCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCcggcgaacgg                >  2:64982/1‑151 (MQ=255)
                                                                                                                                  aTCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCcggcgaacggcg              >  2:817800/1‑151 (MQ=255)
                                                                                                                                  aTCCAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCcggcgaacggcg              <  1:392864/151‑1 (MQ=255)
                                                                                                                                     cAGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCggcgaacggcgaac           >  2:699102/1‑151 (MQ=255)
                                                                                                                                      aGCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAAcc          <  1:499398/151‑1 (MQ=255)
                                                                                                                                       gCTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCa         <  2:1164605/151‑1 (MQ=255)
                                                                                                                                        cTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCa         >  2:551110/1‑150 (MQ=255)
                                                                                                                                        cTGCACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCa         >  2:841716/1‑150 (MQ=255)
                                                                                                                                           cACCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAatc      >  2:907347/1‑150 (MQ=255)
                                                                                                                                            aCCCGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAatcat    >  2:1152081/1‑151 (MQ=255)
                                                                                                                                             cccGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAatcatc   >  2:126465/1‑151 (MQ=255)
                                                                                                                                              ccGCTGATTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCg  <  1:603679/151‑1 (MQ=255)
                                                                                                                                               cGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCg  >  1:610997/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGGCTGTTCACGTACCGCTGACGCTGGAAGCCCAGCTGGAAGCGCGTGCGCTGATGATGTCTACCAACAACATCCTGTCCCCGGCGAACGGCGAACCAATCATCG  >  NC_000913/4186558‑4186850

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: