Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,187,214 A→T D622V (GAC→GTC)  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,2140AT98.4% 243.5 / ‑3.3 64D622V (GAC→GTC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (32/31);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.19e-01

ACTGTTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGG  >  NC_000913/4187066‑4187361
                                                                                                                                                    |                                                                                                                                                   
aCTGTTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTcc                                                                                                                                                   >  2:659150/1‑151 (MQ=255)
aCTGTTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTcc                                                                                                                                                   >  2:677273/1‑151 (MQ=255)
  tGTTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTcc                                                                                                                                                   >  2:187490/1‑149 (MQ=255)
   gTTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTcc                                                                                                                                                   >  1:392811/1‑148 (MQ=255)
   gTTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAg                                                                                                                                                 <  2:511916/150‑1 (MQ=255)
   gTTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAg                                                                                                                                                 <  2:802626/150‑1 (MQ=255)
          gTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGt                                                                                                                                          >  1:1109653/1‑150 (MQ=255)
          gTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTa                                                                                                                                         <  2:756342/151‑1 (MQ=255)
            gCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTaca                                                                                                                                       <  2:44981/151‑1 (MQ=255)
            gCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTac                                                                                                                                        <  1:37618/150‑1 (MQ=255)
             ccATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTacac                                                                                                                                      >  1:17126/1‑151 (MQ=255)
                   tgtgGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCtt                                                                                                                                >  2:715660/1‑151 (MQ=255)
                    gtgGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCtt                                                                                                                                <  1:1108899/150‑1 (MQ=255)
                    gtgGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCt                                                                                                                                 <  1:628351/149‑1 (MQ=255)
                    gtgGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTc                                                                                                                               <  1:801878/151‑1 (MQ=255)
                        atgatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCt                                                                                                                           <  2:560030/151‑1 (MQ=255)
                        atgatTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCt                                                                                                                           >  1:790036/1‑151 (MQ=255)
                                  cGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGcgcg                                                                                                                 >  2:635139/1‑151 (MQ=255)
                                              cTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATc                                                                                                     >  2:39796/1‑151 (MQ=255)
                                                    ccATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGtt                                                                                                 <  2:648435/149‑1 (MQ=255)
                                                    ccATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTgg                                                                                               <  2:818991/151‑1 (MQ=255)
                                                     cATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCACCAAAATGCTGAACGCCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGt                                                                                              <  2:55284/151‑1 (MQ=255)
                                                      aTCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGt                                                                                              >  1:1113033/1‑150 (MQ=255)
                                                       tCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTAt                                                                                            >  2:418593/1‑151 (MQ=255)
                                                        cGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATc                                                                                           <  2:124224/151‑1 (MQ=255)
                                                        cGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATc                                                                                           <  1:549955/151‑1 (MQ=255)
                                                            aaCCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATg                                                                                       <  2:598029/151‑1 (MQ=255)
                                                                aGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACAt                                                                                   >  2:296786/1‑151 (MQ=255)
                                                                aGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACAt                                                                                   >  1:491418/1‑151 (MQ=255)
                                                                 ggCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATg                                                                                  >  1:527060/1‑151 (MQ=255)
                                                                     cTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTc                                                                               <  2:178119/150‑1 (MQ=255)
                                                                           aaaaaaGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCg                                                                         <  2:182788/150‑1 (MQ=255)
                                                                            aaaaaGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCgg                                                                        >  2:1031767/1‑150 (MQ=255)
                                                                             aaaaGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCgg                                                                        >  1:79069/1‑149 (MQ=255)
                                                                                  cAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAgaagaa                                                                 <  1:1005635/151‑1 (MQ=255)
                                                                                   aaTCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAgaagaa                                                                 >  1:526031/1‑150 (MQ=255)
                                                                                      ctcCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAacac                                                             >  2:676/1‑151 (MQ=255)
                                                                                                    acacCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCg                                                >  1:43038/1‑150 (MQ=255)
                                                                                                      acCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAgg                                             >  1:110840/1‑151 (MQ=255)
                                                                                                       ccTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgc                                            >  1:802029/1‑151 (MQ=255)
                                                                                                       ccTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgc                                            <  1:261755/151‑1 (MQ=255)
                                                                                                          gCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcaga                                         <  1:105466/151‑1 (MQ=255)
                                                                                                          gCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcag                                          <  1:377255/150‑1 (MQ=255)
                                                                                                          gCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgca                                           <  1:64968/149‑1 (MQ=255)
                                                                                                          gCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgc                                            >  2:197420/1‑148 (MQ=255)
                                                                                                             aCCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaa                                        <  1:734816/149‑1 (MQ=255)
                                                                                                              ccGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaagca                                     <  1:119982/151‑1 (MQ=255)
                                                                                                              ccGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaagca                                     <  2:641090/151‑1 (MQ=255)
                                                                                                                 cATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaagcagaa                                  <  2:463071/151‑1 (MQ=255)
                                                                                                                  aTTCTCGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGgcagaagcaga                                   >  2:852939/1‑149 (MQ=255)
                                                                                                                    tctcGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGtt                               >  2:609301/1‑151 (MQ=255)
                                                                                                                       cGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGtt                               <  1:528465/148‑1 (MQ=255)
                                                                                                                       cGGTCTGAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCt                            <  1:423900/151‑1 (MQ=255)
                                                                                                                             gAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGaaa                        <  2:430056/149‑1 (MQ=255)
                                                                                                                             gAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAAtt                      >  2:812774/1‑151 (MQ=255)
                                                                                                                             gAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAAt                       <  2:993864/150‑1 (MQ=255)
                                                                                                                             gAAACCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAAt                       <  2:1031673/150‑1 (MQ=255)
                                                                                                                               aaCCGACCGTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTc                     >  2:961685/1‑150 (MQ=255)
                                                                                                                                   gaccgTTATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCa                    >  2:1015821/1‑147 (MQ=255)
                                                                                                                                        ttATTTTTGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGtt           >  2:426783/1‑151 (MQ=255)
                                                                                                                                            ttttGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTcc         >  1:861629/1‑149 (MQ=255)
                                                                                                                                              ttGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAgtc     >  1:898832/1‑151 (MQ=255)
                                                                                                                                               tGCGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAgtct    >  1:949948/1‑151 (MQ=255)
                                                                                                                                                 cGGTCCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAgtctgg  <  1:604715/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACTGTTGGCCGTGCCATTCTGTGGATGATTGTACCGAAAGGTCTGCCTTACTCCATCGTCAACCAGGCGCTGGGTAAAAAAGCAATCTCCAAAATGCTGAACACCTGCTACCGCATTCTCGGTCTGAAACCGACCGTTATTTTTGCGGACCAGATCATGTACACCGGCTTCGCCTATGCAGCGCGTTCTGGTGCATCTGTTGGTATCGATGACATGGTCATCCCGGAGAAGAAACACGAAATCATCTCCGAGGCAGAAGCAGAAGTTGCTGAAATTCAGGAGCAGTTCCAGTCTGG  >  NC_000913/4187066‑4187361

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: