Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091375,0790CT43.8% 74.0 / 33.9 64L134L (CTG→CTAthiBthiamine/thiamine pyrophosphate/thiamine monophosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (4/32);  new base T (28/0);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 6.43e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCT  >  NC_000913/74931‑75220
                                                                                                                                                      |                                                                                                                                               
gCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTc                                                                                                                                                 <  2:782232/149‑1 (MQ=255)
    ggCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGCCCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTg                                                                                                                                           <  1:75060/151‑1 (MQ=255)
       ttGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTAGttgt                                                                                                                                           >  1:533654/1‑147 (MQ=255)
        tgtgGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTc                                                                                                                                        <  2:1024232/150‑1 (MQ=255)
           ggggCGTCATCGCCATAGGCTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTg                                                                                                                                     <  1:763248/150‑1 (MQ=255)
           ggggCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGCCCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCTCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGt                                                                                                                                    <  1:846682/151‑1 (MQ=255)
             ggCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCa                                                                                                                                  <  2:121183/151‑1 (MQ=255)
                gTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTGt                                                                                                                                          >  1:515736/1‑140 (MQ=255)
                  cATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTGTtttttttct                                                                                                                                 >  2:421117/1‑142 (MQ=255)
                   aTCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAtttttttttttttctaaa                                                                                                                              >  1:954496/1‑141 (MQ=255)
                   aTCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTtttttttttttctaaa                                                                                                                              >  2:606469/1‑141 (MQ=255)
                     cGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGaa                                                                                                                          <  2:612331/151‑1 (MQ=255)
                     cGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTTTTTGTTTTTTTTATaaa                                                                                                                              >  2:799302/1‑147 (MQ=255)
                      gCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTAGGTTTCTTTTTTCTTTAACGaa                                                                                                                          >  1:822137/1‑150 (MQ=255)
                       ccATAGACTTTTTGCATCCATAGCAACAGCCCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGaa                                                                                                                          <  2:31790/149‑1 (MQ=255)
                       ccATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTTTTTTTTTCATAAAAGaa                                                                                                                          >  2:3783/1‑149 (MQ=255)
                        cATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTTTTTTTTTCTTTTCATTAAAGaa                                                                                                                          >  1:686883/1‑148 (MQ=255)
                        cATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTTGTtttttttttttcttaaagaaag                                                                                                                        >  2:580973/1‑136 (MQ=255)
                           aGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGCTTTTTTATTTTTTCTTTTCTTACAAGaaaggta                                                                                                                     >  2:552826/1‑145 (MQ=255)
                            gatattgtGCATCCAAAGCAACGGACCCAACCCAGGTGACGTGGTGTGCAGACCCTGAAAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCGTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGaa                                                                                                                    <  1:945735/143‑1 (MQ=255)
                              cTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTTGTTTTTTTTTTTTTTAAACCAAggggga                                                                                                                    >  2:971700/1‑144 (MQ=255)
                               tttttGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTTTTTTTTTCATAAAAAAAGGGGAAg                                                                                                                   >  2:137389/1‑148 (MQ=255)
                                     cATCCATAGCAACAGACCCAGCCCCGGTGTACTGATGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGCTTTTTCCCTTGTTCTTGTCATAAACCCAGGcccagccggc                                                                                                              >  1:108938/1‑138 (MQ=255)
                                       tCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTTGTTTTTTTTTTTTATAAAAAAAAGCGAAGGAGCCCTa                                                                                                           >  2:638834/1‑148 (MQ=255)
                                             gCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTTTTTTTTTTATAAACAAAAGCGCAGTAGGAATAATaaaa                                                                                                     >  1:294433/1‑148 (MQ=255)
                                               aaCAGACCCAGCCCCGGTGTCCTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACg                                                                                                   <  1:804705/148‑1 (MQ=255)
                                                  aGACCCAGCCCCGGTGACTGGGTGCGCGGATCCTGAGAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGTGAAGTAGCCATAATCAAACGGTa                                                                                                <  1:821031/148‑1 (MQ=255)
                                                      ccAGCCCCGGTTTACTGGTGCGCGGATCCTTATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTTTTGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAAc                                                                                                    <  2:563686/140‑1 (MQ=255)
                                                      ccAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAgt                                                                                         <  2:827209/151‑1 (MQ=255)
                                                         gCCCGGGTGTACTGGTGCGCGGTTCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAATAGCCATAATCAAACGGTACGAAAgtgt                                                                                       <  1:474603/150‑1 (MQ=255)
                                                         gCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTGTTTTTGTCTTAAACAAAGGGGGAGTAGGCCCAATCAAAAGGGACggcagag                                                                                        >  1:515420/1‑143 (MQ=255)
                                                         gCCCCGGTGTACGGGTGCGCGGATCCGGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAgtgt                                                                                       <  1:226925/150‑1 (MQ=255)
                                                          ccccGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCa                                                                                     <  2:130300/151‑1 (MQ=255)
                                                           cccGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTTTTTTTTTCATAAACGAAGGCGAAGTAGCCAAAATTAAACAGGATGAAAGGGTCat                                                                                    >  1:691142/1‑151 (MQ=255)
                                                                tgtACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTAGTTTTTTTTTGTTATAAAAGAAAGGGAAGGAGCCATAATCAAACAGGACGaaaagggcattttt                                                                                >  2:873586/1‑139 (MQ=255)
                                                                  tACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTcc                                                                             <  2:722221/150‑1 (MQ=255)
                                                                   tccgggTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTc                                                                              <  1:605522/145‑1 (MQ=255)
                                                                       ggTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAgccg                                                                        <  2:496465/151‑1 (MQ=255)
                                                                         tgcggggccctggATAGATCGCCCGCCAGTTTTGATCGCTGTCAACGATGTCTTTCAGGCTTTGTGGCGGGTTTTGCAGTTTGAGCGGGTCCTAAGGGAAGACGAAGTAGCCATATTCAAGCGGGACGAAAGTGTGATCATTCCAtcc                                                                         <  1:562951/136‑1 (MQ=255)
                                                                          gcgcGGACCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAgccgc                                                                      <  2:161164/149‑1 (MQ=255)
                                                                            gcgGATCTTGATAATTCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCggg                                                                  <  1:391096/151‑1 (MQ=255)
                                                                              ggAGCCTAATAAATCTCCCCCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGa                                                                 <  1:964197/150‑1 (MQ=255)
                                                                                    tGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTTTTTTTTTCCTAAACCAAGGGGCAGGAGGCCTAATTAAACAGGACGGAAGTGTCCTTATTTTAGCCCCCggg                                                                  >  1:780059/1‑143 (MQ=255)
                                                                                        aaaTCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTAGTTTGTTCTTTTCATAAAAGAAGGCGAAGGAGCCATAATTAAACGGGACGAAAGTGTCATTATTTCCACCCCCCGGGAAGTTaa                                                          >  2:494560/1‑147 (MQ=255)
                                                                                           tCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTTGTTTTTTTTTTTTTTAAACCAAGGCGAAGGGAGCATAATTCAACGAGGTGAAAAGGTCATTAATCCAGCCGCCGGGCACGTTAACGGCaac                                                   >  1:212704/1‑149 (MQ=255)
                                                                                           tCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTCAGTTTTTTTTTTTCATAAAACAAAGCGAAGAAGCCCAAATCAAACGGGACGAAAGGGTCATTATTCCAGCCGCCGGGGACCTTAAAGGCATc                                                   >  1:115408/1‑151 (MQ=255)
                                                                                                cGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTCGGTTTTTTTTTTCATAAAAGCAGGGGGAGGGGGCATTATCAAAAGGGACGGAAATGTCATTATTTCAGCCGCCGGGGAAATTTACGGCATcccc                                                >  1:758205/1‑147 (MQ=255)
                                                                                                 gCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTCGTTTGTTTTTTTTATAAAAGAAAGCGGAGGAGCCATAATCAAACGGGACGAAAAAGTGAATATTCCCACCGCCGGGAACGGTAACGGCAAcc                                                  >  2:482359/1‑146 (MQ=255)
                                                                                                   cAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCca                                           <  1:918560/151‑1 (MQ=255)
                                                                                                    aGTTTGGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTGCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCca                                           <  2:873655/150‑1 (MQ=255)
                                                                                                        ttGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACAcc                                       <  2:581329/150‑1 (MQ=255)
                                                                                                        ttGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGCCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACAcc                                       <  2:629440/150‑1 (MQ=255)
                                                                                                        ttGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTTTTTTTTTCTTAAACGAAGGCGAAGGAGCCATAATCAAAAGGGACGAAAGAGGCATTATTTCAACCGCCGGGGACGGTAAAGGGAACCGCTGCCACAcc                                       >  2:700021/1‑150 (MQ=255)
                                                                                                             cGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGGTTTTTCTTTTTATTAACTAAGGCGAAGGAGTCATAATCAAACAGTAGGAAAGAGTGATTTTTTCAGCCCCCGGGGACAGTAACAGGATCCCCAGcccaaccaccttt                                 >  1:850034/1‑140 (MQ=255)
                                                                                                               ctctCAACCATTTCTTTCAGGCTTTGTCGCTGTTTTTTCAGTTTATTCTTGTCATAAACAAAGGCGAAGTAGCCATAATCAAACCCTACGAAACTCTCATTATTCCAGCCCCCGAAAACGTTAACCTCATCCCCTCCCACACCCCtat                                  <  1:86097/148‑3 (MQ=255)
                                                                                                               ctctCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTCAGTTTTTTTTTTTCATAAAAGAAGGCGAAGTAGCCATAATTAAAAGGTACGAAAGGGTCATTATTTCAGCCGCCCGGGACGGTTACGGCATCCCCTGCCCCACCGCctttt                                >  2:656544/1‑149 (MQ=255)
                                                                                                                     aCCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACa                         <  2:996369/151‑1 (MQ=255)
                                                                                                                        aGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTTTTTTTGTCATAAACCAAGGCGGAAGAGGCACAATCAAAAAGGACGAAAGGGTTATTATTTCAGCCGCCCGGGAAGGTTAAGGGATCCGCAGCCCCACCCCCTTTTGCaaaaaat                       >  2:701223/1‑146 (MQ=255)
                                                                                                                         gTTCTTTCAGGCTTTGTGGCGGGTTTTTTAGTTTGTTTTTTTCATAAACGAAGGGGGAGGAGGCCTAATCAAAAAGGACGAAAAAGTCATTATTTCAGCCGCCGGGAAAGTTAACGGCATCCGCCGCCACACCCCCTTTTGCAAACAGAcc                     >  2:1397/1‑151 (MQ=255)
                                                                                                                              ttCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGt                  <  2:60978/149‑1 (MQ=255)
                                                                                                                                      ttGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTggc         <  2:790869/150‑1 (MQ=255)
                                                                                                                                        gtgGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTggcggc      <  2:324247/151‑1 (MQ=255)
                                                                                                                                            cGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCt  <  2:1004805/151‑1 (MQ=255)
                                                                                                                                            cGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCt  >  1:116146/1‑151 (MQ=255)
                                                                                                                                            cGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGCCAAACAGTCCGGTTTTACTGgcggcg     >  1:121574/1‑148 (MQ=255)
                                                                                                                                                      |                                                                                                                                               
GCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCT  >  NC_000913/74931‑75220

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