Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913164,6010TG24.4% 518.5 / 49.5 258intergenic (+67/‑129)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base T (19/151);  new base G (55/0);  total (106/152)
Fisher's exact test for biased strand distribution p-value = 1.57e-36
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 285 total aligned reads displayed.

GAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGC  >  NC_000913/164453‑164751
                                                                                                                                                    |                                                                                                                                                      
gAACGGGCGTTATCCCAAACATGGCTGGCCGGACGACCCGCAAAATACTGCACCCACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGCAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctca                                                                                                                                                      <  2:783270/151‑1 (MQ=38)
  aaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcccc                                                                                                                                                    >  2:1296662/1‑151 (MQ=255)
   gggggtttACCCAAAACATGGGTGCCCGGCCGACCGGCAAAATACGCCCCCGCCGGCACGGCCCAAAATGTTTGCTAAAGGTTGTTTTTCTCCCCCGCCCTGGGGTATTGGGCGGGGGGGGGGGCATCAGCGCGCACGTtcaccctcaccc                                                                                                                                                   <  5:982305/145‑1 (MQ=255)
    gggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAGGTtcacccgcaccc                                                                                                                                                   >  1:1111375/1‑150 (MQ=255)
      gCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcgcccgcacccaa                                                                                                                                                 >  1:1248516/1‑148 (MQ=255)
       cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCACCCAAAc                                                                                                                                               >  8:1153340/1‑151 (MQ=255)
        gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGAAGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCGACCTCACCCGAAcc                                                                                                                                              >  3:1464205/1‑151 (MQ=255)
           aTCCCAACCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTGCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                                                           <  3:121184/151‑1 (MQ=255)
           aTCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGCCCCCCACCCAAAcccac                                                                                                                                           >  4:1226316/1‑149 (MQ=255)
           aTCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCGCCCCCGAAcccac                                                                                                                                           >  5:539737/1‑149 (MQ=255)
           aTCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCGCCCCCACCCGAcccccc                                                                                                                                           >  4:1184396/1‑149 (MQ=38)
           aTCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTACCTGTGGGAGAGGGCCGGGGTTAGGGCATCAGCGCGCACGTGGACCCGGACCCGAAccca                                                                                                                                            >  3:1135434/1‑149 (MQ=255)
           aTCCAAAACTTGTGTGCCCGGACGCCCGGGCAAATACTGCACCGACGCGACGCAGAAAAAAGTATGCGAACGGTTGTTTTTCACCCTCTCCGTGTGGGATTGGGGCGGGGGGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                                                           <  5:220378/151‑1 (MQ=255)
             cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCGCCCCCCAAcccaca                                                                                                                                          >  3:234169/1‑147 (MQ=37)
                   cATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGAAACCACTCCcaca                                                                                                                                    >  3:581449/1‑150 (MQ=255)
                    aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGCCCCGCCCCCGAACCCGCGCCCGCaac                                                                                                                                  >  6:1499694/1‑150 (MQ=37)
                    aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCGCACCCGAACCCACTCCCTCaac                                                                                                                                  >  1:138393/1‑150 (MQ=255)
                    aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCCCCCCCTAACCCTCTCCCGCa                                                                                                                                    >  3:595199/1‑149 (MQ=255)
                     tggcTGGCGGGACGCCCCTGAAATTACTGCACCTCCGCGAGGTACTAAAAAGTATGCGTATGGTTGTTTTTCTCCCTCTCCGTGGGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                                  <  3:1498685/147‑1 (MQ=255)
                     tGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCTGCCCCTAACCCAATCCCACaacg                                                                                                                                 >  3:312447/1‑149 (MQ=38)
                     tGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCAAACCCTCTCCCGCaacg                                                                                                                                 >  1:61557/1‑149 (MQ=255)
                     tGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCCCCCGAACCCTCTCCCGCAAgg                                                                                                                                 >  7:1544857/1‑151 (MQ=255)
                       gcTGGCCGCACGACCCGGCAAATACTGCACCGACGCGTCGGACGAAAACGTATTCGTAAGGTTGGTTCTCACCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                                                                                               <  8:875244/150‑1 (MQ=255)
                        cTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGCCCCGCCCCCGAAACCCCTCCCCCactgggg                                                                                                                              >  6:643497/1‑145 (MQ=37)
                         tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTAACCCTCCCCCTAACCCTCTCCCTCAAGGgagc                                                                                                                             >  8:102340/1‑150 (MQ=255)
                            ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGGGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACACCAACCCCCACCCCCaaccgcgccgg                                                                                                                          >  4:237272/1‑142 (MQ=25)
                            ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAACCACACCCACAACGCAGAggg                                                                                                                          >  4:724189/1‑151 (MQ=37)
                             cGGCCGACCCGGCAAATACTGCCCCGAGGCGACGGACGAAAAAGTATTCGTAAGGTTGGGTTGCTCCCGCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgggg                                                                                                                         <  6:1373285/151‑1 (MQ=38)
                             cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGCCCCGCCCCCGAACCCACTCCCGCAAGGGCGAgggg                                                                                                                         >  8:205184/1‑151 (MQ=37)
                                    cccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCCCGTTCCCCCGCCCCCAAAACCAATACCACAAGCGAGCGGGGGcccctg                                                                                                                  >  7:1190685/1‑147 (MQ=25)
                                    cccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCAAACCCTCACCCGCAAGGGAGAGGGGACCGCTc                                                                                                                  >  3:912815/1‑151 (MQ=255)
                                     cgggaaaaTACTGCACCCACGCGCCGCACCAAAAAGTATTCGTAAGGTGGGTTTGCACCCGCCCCCTGGGGGAGGGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCg                                                                                                                 <  5:1285137/146‑1 (MQ=255)
                                     ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGAACCCACGACCGCAAGGGAGCGGGGGCCgcg                                                                                                                   >  2:223865/1‑147 (MQ=37)
                                      cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGAACCGCACCCGAAACCTCTCCCCCAACGGAGAGGGGACCGCTCgc                                                                                                                >  8:313415/1‑150 (MQ=37)
                                         cAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGCCCCTCACCCGAACCCTCTCCCGGAAGGGAGCGGGGCCCGCTCgcgcc                                                                                                             >  6:116907/1‑150 (MQ=37)
                                            aTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCGTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACaa                                                                                                          <  3:662731/151‑1 (MQ=255)
                                                ggccccgccgCGCCGGCCGAAAATGTATGCGTAGGGTGGTTTTGCGCCCGCCCCGTGTGGGTGTGGGCGGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAttt                                                                                                      <  8:1473695/147‑1 (MQ=255)
                                                  cACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                                                                    <  2:133119/151‑1 (MQ=255)
                                                    ccGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                                                                  <  1:1472722/151‑1 (MQ=255)
                                                        cgcgacgCACGAAAAAGTATTCGTAAGGTGGGTTTTCTCCGTCTCCCTGTGGGCGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                                                              <  2:1473422/151‑1 (MQ=255)
                                                          cgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCGCCCCCAAACCCACTCCCGCAAGGGAGGGGGGGCCGCGCGCGCCCAAAATTTTAGAGata                                                                                             >  2:951849/1‑150 (MQ=25)
                                                          cgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCTAACCCACTCCCACAAGGGAGAGGGGCCCGAGAGAGCACAAAATTTTAGAGata                                                                                             >  4:638807/1‑150 (MQ=37)
                                                           gagggactAAAACTTATTGTTAAGGTGGGTTTCCTCCCTCTCCCTGGGGGAGTGGGCGGGGGGGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                                           <  2:438830/148‑1 (MQ=255)
                                                           gacggactAAAACGTTTTCTTAAGGTTGGTTTTCTCCCTCTCCTTGTGGTAGAGGGCCGGGTTGAGGGAATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                                           <  2:1435653/151‑1 (MQ=255)
                                                           gacggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCCCCCAAACCCACAAACACAAAGGAGAGGGGGCCCCTCGACCACAAAATTTGAGAGctat                                                                                            >  2:971860/1‑150 (MQ=25)
                                                            acggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtct                                                                                          >  8:630673/1‑151 (MQ=255)
                                                              gtacgAAAATGTTTTCTTATGGTTGGTTTTCTCCCTCTCCGTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttc                                                                                        <  4:499033/149‑1 (MQ=38)
                                                               gacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttct                                                                                       <  6:933978/151‑1 (MQ=255)
                                                                acgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCCCCCGAACCCTCTCCCGCAAGGGAGAGGGGCCCGCTCGAGCACAAATTTTTATAGATAtattaat                                                                                      >  3:1368359/1‑148 (MQ=37)
                                                                acgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCACTCCCGCAAGGGAGAGGGGGCCGCTAGAGCACAAATTTTTAGAGATAtcatctt                                                                                      >  1:1189775/1‑151 (MQ=37)
                                                                 cgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                                                                     <  6:177648/151‑1 (MQ=255)
                                                                  caaaaaGTATTCGTACGGTTGGTTTGCTCCCTCACCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                                                    <  5:677525/150‑1 (MQ=255)
                                                                    aaacgtctgcGTACGGTTGGGTTTCTCCCGCTCCCTGTGGGAGAGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                                                                  <  1:675816/142‑1 (MQ=39)
                                                                       atttttgttagggttgtttTGCGCCCGCGCCGTGTGGGAGGGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                               <  7:387691/140‑1 (MQ=255)
                                                                       aGTATTCGTAAGGTTGGTTTTCTCCCTCTCCGTGGGGTTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCGTCTTCTGTCtt                                                                               <  2:411615/151‑1 (MQ=255)
                                                                       aGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCACCCCCTAACCCTATCCCCCAACGGAGAGGGGACCGATAGAGCACAAAATTTGAGAGATATCTTCTTCTGTCtt                                                                               >  1:823693/1‑151 (MQ=37)
                                                                         ttttcttattgttgcttTTCTCCGTCTCCATTTGTTATTGGTCCGTGTTGAGGGCATCATCGCGCAAGTTCACCCAGAGCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                                              <  3:59551/140‑1 (MQ=255)
                                                                          catCGTATGGGTGGTTTTCTCCCTCTCCCTGTGGTATAGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                                             <  2:186534/148‑1 (MQ=39)
                                                                           ttCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTACCACAAGGGAGAGGGGACCGCTCGAGCACAAATTTTTAGAGATATATTGTTCTGTCATGTa                                                                            >  3:80789/1‑150 (MQ=37)
                                                                             gttAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                                         <  3:418151/149‑1 (MQ=255)
                                                                               tatggtgggttTTCTCCCGTTCCGTCTGGGTGAGGGCCGGGGTTAGGGCATCCGCGCGCACGTTCACCCTCACCCTAACGCTCTGCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                                          <  4:1459954/145‑1 (MQ=255)
                                                                                  ggtgggttTTCTCCGTATCCGTGTGGGCCTGGGCGGGGGGCCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                                                    <  5:1183687/151‑1 (MQ=255)
                                                                                  gggtggttTGCGCCCTCTCCCTGGGGGAGAGGGCCGGGGGGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                                                    <  5:639737/148‑1 (MQ=255)
                                                                                    ttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCCCCCCAACCCGCTCCCGCCAGGGCGGGGGGGCCGCGGGCGCCCACATTTTTTGCCTTATCTTCTTGTGTGGTGGTACAGAAGAAc                                                                  >  5:714230/1‑151 (MQ=12)
                                                                                     tggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCACACCCGAACCCCCACCCTCAAGGGAGAGGGGGCCGCAGGACGACCAAATTTTAGAGATATATTGTTGTGTCTTGTTAAAGAAGAAc                                                                  >  1:275512/1‑150 (MQ=18)
                                                                                     tggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCGCACCCGAACCCACTCCCCCAACGGGGCGGGGCCCGCTCGCGCCCAAATTTTTGGGGATTTCTTCTTCTGTGGTGGTTCAgaagaaa                                                                  >  8:1394075/1‑149 (MQ=12)
                                                                                              ccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCGCCCCCGAAACCCCTCCCGCAAGGGGGGGGGGGCCGCTCCCGCCCAAAATTTTTGAGATTTATTATTCGGGGGTGTTAAAGCAGAACACCAAATCgg                                                        >  6:961159/1‑151 (MQ=9)
                                                                                                ctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCAATACCACAACGGAGAGGGGCCCCATACCGCACAAAATTTTATAGAAATAATATTATGTATTGTAAAAGAAGAACAGAAAAtacggc                                                      >  1:5207/1‑146 (MQ=16)
                                                                                                  ctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGCCCCGCCCCCAAAACCACTCCCTCCACGGGGTGGGGGCCGCTGGCGCACCAATTTTTTGAGATATCTTATTGTGTGTTGTTAAAGCACAAAAGCAAATAgggggt                                                    >  6:563647/1‑148 (MQ=12)
                                                                                                  ctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCCCCCGAACCCACTCCCGGAAGGGGGAGGGGCCCGCTCGCGCCCAAATTTTTAGAGCTATATTCTTGTTGGGTGTTACACAAGAACAGAAAATAGGGcgt                                                    >  6:837923/1‑149 (MQ=16)
                                                                                                   tccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAAACCACTCCCTCCAGGGCGAGGGGGCCGATCGCGCACAAAATTTGAGAGCTATCTTAATCTTTGATGGAACAGAAGAAccccaaaaagggcgtt                                                   >  2:1022438/1‑136 (MQ=19)
                                                                                                   cccctGGGGGAGCGGGCGGGGGTCTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                                                   <  3:1081743/150‑1 (MQ=255)
                                                                                                     cctggggtcgagGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                                                 <  4:806706/142‑1 (MQ=255)
                                                                                                      gtgtgGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                                                <  2:876788/150‑1 (MQ=255)
                                                                                                       tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTAACCCACCACCGAAACCTAGACCACAAGGGAGAGGGGGCCGCACGAGCACCAATTTTTAGAGATATCTTCATATGTGTTGTATCATAATAAAAGAAAATAGGGCTTTTg                                                 >  3:1425516/1‑149 (MQ=16)
                                                                                                        gggggggggggggggggggCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcg                                               <  7:25419/137‑1 (MQ=255)
                                                                                                          gtgagaGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                                            <  6:1129077/149‑1 (MQ=255)
                                                                                                            gcttgggCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                                          <  3:1273099/147‑1 (MQ=255)
                                                                                                              gagGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAat                                        <  1:377224/151‑1 (MQ=255)
                                                                                                              gagGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCACCCCCCAAACCCCTCCCGCCAGGGGGCGGGGGGCGCGCGCGCCCCAAATTTTTTAGATATTTTATTTTGGGGTGTTGGAGAACAACAGAAAAAAGGGCGTTTggggcggcg                                         >  8:510636/1‑142 (MQ=255)
                                                                                                                gggCGGGGGGGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat                                      <  8:1269847/151‑1 (MQ=255)
                                                                                                                gggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat                                      <  5:231526/151‑1 (MQ=255)
                                                                                                                 ggCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt                                     <  7:1320241/151‑1 (MQ=255)
                                                                                                                  gCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg                                    <  3:963168/151‑1 (MQ=255)
                                                                                                                    ggggtggcggccATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                  <  7:50605/140‑1 (MQ=255)
                                                                                                                       ggTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgag                               <  8:565245/151‑1 (MQ=255)
                                                                                                                         tGAGGGCATCAGCGCGCACGTTCACCCGCACCCGAACCCAATCCCTCAAAGGAGAGGGGCCCGATCGCGCACAAAATTTTAGAGATATCTTATTGTTTGTTGTAAAAGAAGAACAAAAAATAGGGCGTTTGGGCCGGAATATTTCGgaga                              >  1:333283/1‑150 (MQ=255)
                                                                                                                         ggcggccATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaa                             <  6:219118/145‑1 (MQ=255)
                                                                                                                            gggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa                          <  1:204733/151‑1 (MQ=255)
                                                                                                                                aTCAGCGCGCAGGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCAt                      <  4:1130597/151‑1 (MQ=255)
                                                                                                                                  cAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATgg                    <  2:1190076/151‑1 (MQ=255)
                                                                                                                                  cAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATgg                    <  4:388598/151‑1 (MQ=255)
                                                                                                                                  cAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATgg                    <  5:599405/151‑1 (MQ=255)
                                                                                                                                   aGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGc                   <  4:679146/151‑1 (MQ=255)
                                                                                                                                      ccgcACGTGCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCgg                <  5:410773/150‑1 (MQ=255)
                                                                                                                                         cACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGaa             <  8:1516169/151‑1 (MQ=255)
                                                                                                                                          aCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAAt            <  5:1551329/151‑1 (MQ=255)
                                                                                                                                              tcaccctcacccTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGAc         <  5:324418/150‑1 (MQ=255)
                                                                                                                                               caccctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGAcc        <  5:869241/150‑1 (MQ=255)
                                                                                                                                               caccctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACcg       <  1:999489/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGC  >  NC_000913/164453‑164751

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑