Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091375,0790CT41.1% 94.6 / 45.3 73L134L (CTG→CTAthiBthiamine/thiamine pyrophosphate/thiamine monophosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (4/39);  new base T (30/0);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 1.66e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGG  >  NC_000913/74937‑75226
                                                                                                                                               |                                                                                                                                                      
cTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGttt                                                                                                                                                   >  2:536301/1‑148 (MQ=255)
   gtgGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTggtttt                                                                                                                                             >  2:385505/1‑146 (MQ=255)
    tgGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTGTtttt                                                                                                                                             >  1:317056/1‑150 (MQ=255)
     ggggTGTCATGGCCGTGGACTTTTCGCACCCATAGCAACAGACCCACCCCCGGTGTACTGGTGCGCGAATCCTGATAAATCGCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGt                                                                                                                                           <  1:69523/151‑1 (MQ=255)
      gggCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTc                                                                                                                                          <  1:838085/151‑1 (MQ=255)
          gTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATaa                                                                                                                                      <  2:73309/151‑1 (MQ=255)
          gTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTGTTTTTGtttt                                                                                                                                        >  2:139457/1‑146 (MQ=255)
            cATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAAc                                                                                                                                    <  1:737450/151‑1 (MQ=255)
            cATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAAc                                                                                                                                    <  1:1000676/151‑1 (MQ=255)
               cGCCATAGCCTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGaa                                                                                                                                 <  2:112384/151‑1 (MQ=255)
               cGCCATAGACTTTTTGCATCCATAGCAACAGCCCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGaa                                                                                                                                 <  2:275482/151‑1 (MQ=255)
               cGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTCCTGGTGCGCGAATCCTGATAAATCGCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACg                                                                                                                                   <  2:826180/149‑1 (MQ=255)
                 ccATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAgg                                                                                                                               <  1:590080/151‑1 (MQ=255)
                  cATAGACTTTTTGCATCCATAGCAACAGACACAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTCAAATTTTTCtttttgttaaaaaagg                                                                                                                               >  1:427683/1‑136 (MQ=255)
                                  ccATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCa                                                                                                             <  2:1038979/151‑1 (MQ=255)
                                    aTAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTTTTTTTTTTATAAAAGAAGGCGAAGGAGCCATAATCaa                                                                                                            >  1:902584/1‑150 (MQ=255)
                                      aGCAACAGCCCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAAc                                                                                                          <  2:1023265/150‑1 (MQ=255)
                                        cAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTTTTTTTTTTCTCAACAAAAGAGAAGCAGCCATAAACAAAAgg                                                                                                        >  2:125281/1‑150 (MQ=255)
                                          aCAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTGTTTTTGTCATAAAAAAAAGCGAAGGAGCCATAATCaaa                                                                                                           >  1:1018061/1‑145 (MQ=255)
                                          aCAGACCCACCCCCGGTGTCCTGGTGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGt                                                                                                       <  2:730035/149‑1 (MQ=255)
                                           cAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTTTTTTTTTTTTAAACAAAAGCGGAGGAGCCAAAAACAAAAgg                                                                                                        >  2:571545/1‑147 (MQ=255)
                                            aGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTTATTTTTTTTTTTCATAAAAAAAGGGGAAGGAACCACAATCaaa                                                                                                           >  2:1065760/1‑143 (MQ=255)
                                            aGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTAGTTTTTTTTTGTTTTAAACAAAGGGGAAGGAGCCCTAAACAAAAGGTACg                                                                                                    >  2:1124873/1‑150 (MQ=255)
                                               cccACCCCCGGTGTACTGTGGCGCGGAGCCTGATAAATTCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAg                                                                                                <  2:413275/151‑1 (MQ=255)
                                                  acccccGGTGTCCTGGTGCGCGGATCCTGATAAATCGCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAgt                                                                                               <  1:300398/147‑1 (MQ=255)
                                                   gCCCCGGTTTTCTGGTTCCGGTATCTTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAgtgt                                                                                             <  2:856892/150‑1 (MQ=255)
                                                    cccgGGTGTAGTGGGGCGCGATTCTTCATAAAGCTCCCGCCAGTTTTAATCGCTCTCAGCCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCAGAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCa                                                                                           <  2:953720/151‑1 (MQ=255)
                                                       cGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTATTTTTTTTTTTTTTTAAAAAAAGGCGGAGGAGGCATAATTAAAAAGGACGaaaaagttctt                                                                                         >  2:495835/1‑142 (MQ=255)
                                                        ggTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTAGTTTTTTTTTTTTATAAAAGAAAGCGAAGGAACCCTAATAAAACAGGAAGAAAGAGTGatt                                                                                         >  1:24030/1‑149 (MQ=255)
                                                        ggTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCattat                                                                                      <  1:945630/151‑1 (MQ=255)
                                                        ggTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTCTTTTTTTTTTTCATAAAAAAAAGGGAAGAAGCCAAAATCCAAAGGTACGacaaagttatt                                                                                         >  2:1102781/1‑139 (MQ=255)
                                                          tgtACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCattatt                                                                                     <  1:625515/150‑1 (MQ=255)
                                                            tACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTTTTTTTTTTTTTAAAAAAAGCGGAAGAGCCATAAACAAAAAGGAAGAAAGGGTCattatttc                                                                                   >  2:1127986/1‑148 (MQ=255)
                                                            tACTGGTGCGCGGACCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTcc                                                                                   <  1:46274/150‑1 (MQ=255)
                                                             aCTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTAGTTTTTTTTTTTCATAAACGAAGGGGGAGAAGCCCAAAACAAAAGGGGCGaaaaaggtattatt                                                                                      >  1:1106227/1‑136 (MQ=255)
                                                             aCTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTATTTTTTTCTTTTCTTAAAAGAAGGGGGAGTAGCCCAAATCAAAAGGTGAGaaaaaggtattatt                                                                                      >  1:377444/1‑136 (MQ=255)
                                                                ggggcgcGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAgc                                                                                <  2:845962/146‑1 (MQ=255)
                                                                ggTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAgccg                                                                              <  2:72708/151‑1 (MQ=255)
                                                                   gcgcGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAgccgcc                                                                            <  2:931874/150‑1 (MQ=255)
                                                                   gcgcGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAgccgcc                                                                            <  2:423990/150‑1 (MQ=255)
                                                                        gATCCTGAAAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAAc                                                                      <  2:549933/151‑1 (MQ=255)
                                                                          tCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTTTTTTTTTTTTTCATAAAAGAAGGCGAAGGAGCCATAATTAAACGGGACGAAAGTGTCATTTTTTCAACCGCCcgg                                                                         >  1:373444/1‑146 (MQ=255)
                                                                           ccTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGtt                                                                  <  2:317014/151‑1 (MQ=255)
                                                                                 aaaTCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTa                                                                 <  1:281729/146‑1 (MQ=255)
                                                                                 aaaTCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTTTTTTTGTCTTAAACGAAAGCGAAAAAACCCTAAATAAAAAGTAAGAAAAAGTCATTATTTTAGCCCCCGGGGACggt                                                                  >  2:1046374/1‑143 (MQ=255)
                                                                                    tCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTTTTTTTTTCATAAAAGAAAGGGAAGGAGCCATAATAAAAAGGGGCGAAAGAGTCATTATTTCAGCCGCCGGGGACGGTTACgg                                                             >  1:1124444/1‑147 (MQ=255)
                                                                                     cccccGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATcc                                                        <  1:747635/149‑1 (MQ=255)
                                                                                     cACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATcc                                                        <  2:448468/151‑1 (MQ=255)
                                                                                     cACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATcc                                                        <  1:345280/151‑1 (MQ=255)
                                                                                      aCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTCGTTTTTTTTTTTCTTAAACGAAGGGGAAGGAACGCCAAACAAACAGGAAGAAATTGTCAATATTTCTACCCCCCGGGAAGCTAACgg                                                             >  2:361013/1‑145 (MQ=255)
                                                                                       cccGCCATTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTGCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGc                                                      <  2:426413/151‑1 (MQ=255)
                                                                                        ccGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTTTTTTTTTTCTAAAAAAAGGGGAAGGAACCATAAACAAAAAGGACGAAAAAGTCATTATTTCACCCGCCCGGGAAGCTAAAGGGATcc                                                        >  1:883951/1‑148 (MQ=255)
                                                                                          gCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTAGTTTTTTCTTTTTATTAAAGAAGGCGAAGGAGGCCTAATTAAACAGTACGAAAGAGTGATTATTTCAGCCGCCCGGGACGGTAACGGCAAcc                                                        >  2:547492/1‑146 (MQ=255)
                                                                                          gCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTGTTTTTGTTTTAAACGAAAGCGAAGTAGCCATAATCAAACGGGACGAAAGTGTCATTATTTCAGCCGCCGGGGAAGGTTACGGCATcccccgc                                                   >  1:62986/1‑146 (MQ=255)
                                                                                           ccAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGcc                                                  <  1:505736/151‑1 (MQ=255)
                                                                                            cAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCca                                                 <  1:926438/151‑1 (MQ=255)
                                                                                                 ttGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTAGTTTTTTTTTTTATAAAAAAACGGGGAGGAGACATAAATAAAAAGGACGAAAGAGGCTTTATTTCAGCCGGCGCGGAAGGTAAAGGCCTCCgccgcccc                                                >  1:618641/1‑145 (MQ=255)
                                                                                                  tGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTTGTTTTTTTTTTTCATAAAAAAAGGCGGAGTAGCCCTAATCAAACGGGACGAAAGAGTCATTTTTTCTTCCCCCCGGGAGATTAAAGGCATCCGCCGCCCCAcccc                                           >  1:1121339/1‑149 (MQ=255)
                                                                                                   gATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGc                                           <  2:1132160/150‑1 (MQ=255)
                                                                                                    aTCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTTTTTTTTTTTTAAAAGAAGGCGGAGGAGCCATAAAAAAAAAGGGCGGAAAGGTTATTATTTCCGCCGCCGGGGACGGTTAAAGGATCCGCTGccccccccc                                           >  1:894073/1‑142 (MQ=255)
                                                                                                     tCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTTTTTTTTTTATAAAAAAAAGGGGAGGAGGCATAATTAAAACGGACGGAAGGGGCATTATTTCAACCGCCCGGGAAGTTAAAGGCCTccccccccccccccc                                           >  2:130889/1‑135 (MQ=255)
                                                                                                       gCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTTTTCTTTTTCTAAAAAAAGGAGAAGGAGGCATAATTAAAAGGTACGGAAGAGTTATTATTTCTGCCGCCGGGGAAGGTAACCGCCTCCCCCGCCCCACCCCcttt                                       >  2:63128/1‑150 (MQ=255)
                                                                                                        ctctCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTGTTTTTTTTATAAAAAAAAGGGGAGAAGGCATAATTAAACGGGACGGAAGAGGCATTATTTCAACCGCCGGGGAAGGTAAAGGCATCCGCTGCCACCCCCCcttt                                       >  2:818235/1‑149 (MQ=255)
                                                                                                         tctcAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTTTGTTTTTTTTTTTTATAAACGAAAGCGGAGGAGCCCTAATCAAAACGGTCGGAAGGGTCCTTTATTCAGCCCCCGGGGAAGGTAACGGCATCCGCCGCCACACCCCtttt                                       >  2:201495/1‑148 (MQ=255)
                                                                                                            cAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATGATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGc                                    <  1:537213/148‑1 (MQ=255)
                                                                                                                   ttCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTcc                           <  2:948300/150‑1 (MQ=255)
                                                                                                                         cAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGtttt                     <  2:993414/150‑1 (MQ=255)
                                                                                                                          aGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTAc                   <  2:502938/151‑1 (MQ=255)
                                                                                                                           ggCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACt                  <  2:645848/151‑1 (MQ=255)
                                                                                                                                       gggTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCt        >  1:634949/1‑150 (MQ=255)
                                                                                                                                           ttttCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAAc     >  1:103269/1‑149 (MQ=255)
                                                                                                                                           ttttCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACa    <  1:489604/150‑1 (MQ=255)
                                                                                                                                            tttCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCAGCGTCTAACAgg  >  1:50194/1‑151 (MQ=255)
                                                                                                                                               |                                                                                                                                                      
CTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGG  >  NC_000913/74937‑75226

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