Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913802,1180AC43.0% 140.8 / 18.6 86T78P (ACC→CCC) ybhIputative DASS family tricarboxylate or dicarboxylate transporter
Reads supporting (aligned to +/- strand):  ref base A (38/11);  new base C (0/37);  total (38/48)
Fisher's exact test for biased strand distribution p-value = 8.25e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGGTCTGG  >  NC_000913/801969‑802268
                                                                                                                                                      |                                                                                                                                                       
tAAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTaaaa                                                                                                                                                         <  2:1002558/150‑1 (MQ=255)
  aGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAacca                                                                                                                                                     <  1:835534/151‑1 (MQ=255)
   gCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAaccac                                                                                                                                                    >  4:441547/1‑151 (MQ=255)
    cGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAaccacc                                                                                                                                                   <  4:95179/151‑1 (MQ=255)
     gAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCAccg                                                                                                                                                  <  2:221576/151‑1 (MQ=255)
       aCTGGCGTGGGTGCTTTTTGGGATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAAACAccgcc                                                                                                                                                >  2:302133/1‑151 (MQ=255)
            cGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTAtt                                                                                                                                           >  2:578174/1‑151 (MQ=255)
             gTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTAtt                                                                                                                                           >  3:363602/1‑150 (MQ=255)
               gggTGTTTTTTGGTATTTCCCTGGCGGCCATTGTGGGGCTGGTTATCAACCTTTTCCCGGAACCTGTCTTCTTTTTATTTCCCGTTGCTGCCTCAATGGTGGTGGTGGTTAACTTACCCGACGGTGCTTTAAAACCCACCGCCGTATTa                                                                                                                                          <  2:961771/149‑1 (MQ=255)
                   gCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCgg                                                                                                                                     >  2:583249/1‑150 (MQ=255)
                     tttttGGTATTTCCCTGGCGGCCATTGTGGGGCTGTTTATCAACCTTTCCCCGGACCCTGTCGTATTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTACCCGACGGTGCTTTAAAACCCACCGCCGTATTAAGCgg                                                                                                                                     <  2:984636/148‑1 (MQ=255)
                     ttttgGGTATTTCCCTGGCGGCCTTTGTGGGGCTGTTTATAAAGCCTTCCCCGGACCCTGTCGTACTGTAAATTGCCGTTGCTGCCTCAAGGGTGGTGGTCGTTAATTTACCCGACGGTGCGTTAAAACCCACCGCCGTATTAAGCGGtt                                                                                                                                   <  2:585294/150‑1 (MQ=255)
                      ttttGGTTTTTACCGGGGGGCTTTTGTGGGGCTGTTTATCAAGCCTTTCCCGGACCCTGTCGTACTGTTATTTGCCTTTGCTGCCTCAATGGTGGTGGTGGGTAACTTATCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTa                                                                                                                                  <  3:260416/150‑1 (MQ=255)
                         tGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTActct                                                                                                                              >  1:891570/1‑151 (MQ=255)
                              tttgcgggggggcCGTTGAGGGGCTGTTTATAAGGCATTCCCGGGATCCTGTCTTACTTTTTATTCCCGTTGCTGCCCCAGGGGTGGTGGTCGTTAACTTACCCGACGGTGCGTTTACAACCGCCGCCGTATTAAGCGGTTACTCTTc                                                                                                                            <  1:407472/138‑1 (MQ=255)
                              tttccgtGGCGGCTATTGTGGGGCTGGTTATCAAGCCTTCCCGGGACCCTTTCGTACTTTTAATTGCCTTTGCTGCCTCAATGGTGGTGTTCGTTACCTTATCCGCCGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCa                                                                                                                           <  3:63661/143‑1 (MQ=255)
                              tttCCCTGGCGGCCATTGGGGGGCTGGTTATCAAGCCTTTCCCGGACCCTGTCTTACTTTTATTTGCCGTTGCTGCCTCAATGGTGGTGGTCGTTAATTTACCCGACGGTGCGTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCa                                                                                                                           <  2:523536/149‑1 (MQ=255)
                               ttcccgggcggcCTTGGTGGGGTTGGTTATCACGCCTTTCCCGGACCCTGTCGTCCTGTTAATTGCCTTTGCTGCCTCAAGGGGGGTGGTCGTTACTTTACCCGAGGGTGGTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAgg                                                                                                                         <  3:21338/144‑1 (MQ=255)
                               ttACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGt                                                                                                                        >  3:282318/1‑151 (MQ=255)
                                    tGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTaccac                                                                                                                   >  2:1199762/1‑151 (MQ=255)
                                    tGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTaccac                                                                                                                   >  1:845485/1‑151 (MQ=255)
                                    tGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTaccac                                                                                                                   >  2:64407/1‑151 (MQ=255)
                                    tGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTaccac                                                                                                                   >  1:591571/1‑151 (MQ=255)
                                        ggcCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctgg                                                                                                               >  2:1047564/1‑151 (MQ=255)
                                          ccATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctggct                                                                                                             >  2:417142/1‑151 (MQ=255)
                                           cATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACctggctg                                                                                                            >  3:71872/1‑151 (MQ=255)
                                             ttGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGgt                                                                                                          >  3:161553/1‑151 (MQ=255)
                                              tttgGGGCTGTTTATCATCCCTTTCCCGGACCCTTTCGTCTTGTTATTTGCGGTGCCTCCCTCATTGGTGTTGGTCGGTATTTTCTCCGCCGTTCCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGgt                                                                                                          <  1:119800/148‑1 (MQ=255)
                                                tgGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGtt                                                                                                       >  1:880358/1‑151 (MQ=255)
                                                 ggggCTGTTTATAAAGCCTTTCCCGGAACCTGTCGTACTTTTATTTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTATTTTACCCGAGGGTGCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGtt                                                                                                       <  2:722381/150‑1 (MQ=255)
                                                  gggTTGTTTATCAAGCCTTTCCCGGACCCTGTCGTACTGTTAATTGCGGTTGCTGCCTAAATGGTGGTGGTCGTTACTTTATCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGtttc                                                                                                     <  1:877937/151‑3 (MQ=255)
                                                  gggTTGTTTATCAACCTTTCCCGGGACCCTGTCTTACTTTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTGGTTAATTTACCCGGCGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTtct                                                                                                     <  2:175803/151‑1 (MQ=255)
                                                  gggCTGTTTATCAACCTTTTCCCGGACCCTGTCTTACTGTTAATTGCCGTTGCGGCCTAAATGGTGGTGGGCGGTACTTTGCCCGAGGGGGCTTTAAAAACCTCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTtct                                                                                                     <  1:906694/151‑1 (MQ=255)
                                                   ggCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAATTTATCCGACGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTtctc                                                                                                    <  4:260516/151‑1 (MQ=255)
                                                   ggCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTtctc                                                                                                    >  2:1165139/1‑151 (MQ=255)
                                                     cTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCgg                                                                                                  >  1:1251014/1‑151 (MQ=255)
                                                     cTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCgg                                                                                                  >  2:944793/1‑151 (MQ=255)
                                                      tGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGc                                                                                                 >  2:1217702/1‑151 (MQ=255)
                                                         ttATCAAGCCTTCCCCGGACCCTGTCGTTCTGTTAATTGCGGTTGCTGCCTCAATGGTGGTGGTCGGTATCTTATCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCg                                                                                                   <  2:503883/146‑1 (MQ=255)
                                                           aTCAAGCTTTTCCCGGACCCTGTCTTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGGGGTCGGTATTTTACCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGttt                                                                                             <  2:245118/150‑1 (MQ=255)
                                                           aTCAAGCCTTCCCCGGACCCTGTCGTAGTGTTAATGGCCTTTGCTGCCTAAAGGTGGGGGGTGGTTACCTTACCCGACGGTGCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTa                                                                                            <  2:1269357/151‑1 (MQ=255)
                                                                 ccTTTCCCGGACCCTGTCGTACTTTTAATTCCCTTTGCTGCCTCAATGGTGGTGGTCGGTAATTTACCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTa                                                                                       <  1:167348/150‑1 (MQ=255)
                                                                 ccTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTa                                                                                       >  4:157532/1‑150 (MQ=255)
                                                                   ttccccGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTAAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTa                                                                                       <  2:391147/145‑1 (MQ=255)
                                                                    ttCCCGGAACCTGTCGTACTGTTAATTCCCGTTTCTCCCAAATTGGGGGTGTTCGATAACTTATCCGGCGGTGCTTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTa                                                                                       <  2:977429/147‑1 (MQ=255)
                                                                     tCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgcgc                                                                                  >  2:990042/1‑151 (MQ=255)
                                                                      cccGGAACCTGTCGTACTTTTAATTGCCGTTGCTGCCTAAATGGTGGTGGTCGGTAATTTATCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgcg                                                                                   <  4:158060/149‑1 (MQ=255)
                                                                        ccgaacccTGTCTTACTGTTATTTGCCGTTGCTCCCTAAAGGGTGGTGGTCGTAATTTTACCCGCCGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGcgccg                                                                                <  1:1237292/145‑1 (MQ=255)
                                                                        cGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGcgccgc                                                                               >  4:267980/1‑151 (MQ=255)
                                                                            aaCCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAtt                                                                            >  2:1107035/1‑151 (MQ=255)
                                                                            aaCCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGGGTTCTCGGCGTTTACCTTAAGCGCCGCAt                                                                             >  2:828961/1‑150 (MQ=255)
                                                                              ctcgtggTATTGTTAATGGCCGTTGTTCCCTATGTGGTGGTGGTGGATATCTTTCCCGAGGGTGAGTTAAAACCCCCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAttt                                                                           <  2:1058843/144‑1 (MQ=255)
                                                                              ccTGTCGTACTTTTAATGGCCGTTGCTGCCTCAATGGTGGTGGTGGGTAACTTATCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAttt                                                                           <  4:133383/150‑1 (MQ=255)
                                                                               cTGTCGTTCTGTTAATGCCCTTTGCTGCCTCAATGGTGGTGTTGGGTAACTTACCCGGGGGTGCTTTTAACACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCAttt                                                                           <  1:405807/149‑1 (MQ=255)
                                                                                tGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtgtg                                                                        >  2:115748/1‑151 (MQ=255)
                                                                                 gTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtgtg                                                                        >  1:1071355/1‑150 (MQ=255)
                                                                                  tCGTTCTTTTAATTGCCGTTGCTGCCTAAATGGTGGTGGTCGGTAACTTACCCGACGGTGCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtgtg                                                                        <  1:1202073/149‑1 (MQ=255)
                                                                                  tCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGac                                                                      >  1:1117150/1‑151 (MQ=255)
                                                                                    gTACTGTTAGTTGCCTTTGCTCCCTCAAGGGTGGTGGTCGTTAATTTATCCGACGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGacca                                                                    <  2:1270625/151‑1 (MQ=255)
                                                                                     tACTGTTAATTGCCTTTCCTGCCTCATTGGTGGTGGTCGTTACCTTATCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGa                                                                       <  3:27213/147‑1 (MQ=255)
                                                                                        tGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCgg                                                                >  3:239695/1‑151 (MQ=255)
                                                                                          ttAATTGCCGTTGCTCCCTCAAGGGTGGTGGTCTGTAACTTACCCGACGGGGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCg                                                                 <  2:245672/148‑1 (MQ=255)
                                                                                            aaTTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGttt                                                             >  4:156110/1‑150 (MQ=255)
                                                                                            aaTGGCCGTTGCTGCCACAATGGTGGTGGTCGGTAACTTAGCCGTGGGTGCTTTTAAACCCCCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGtt                                                              <  1:583363/149‑1 (MQ=255)
                                                                                              ttGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAgg                                                          >  2:355482/1‑151 (MQ=255)
                                                                                               tGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGt                                                         >  2:786769/1‑151 (MQ=255)
                                                                                                gCCGTTGCTGCCTCAGTGGTGGTGGTGGGTAACTTATCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGt                                                         <  3:442367/150‑1 (MQ=255)
                                                                                                 cctttGCTGCCTCAATGGTGGTGTTGGGTAATTTACCCGACGGTGCTTTAAAACCCGCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAgg                                                          <  1:885712/145‑1 (MQ=255)
                                                                                                    ttgcgccctgatgggtggtggtCGTTACCTTATCCGACGGTGTTTTGAAAACCTCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAAc                                                     <  1:302180/137‑1 (MQ=255)
                                                                                                         gcCTCAATGGGGGTGGTGGTTAACTTATCCGACGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAtt                                                <  4:161555/150‑1 (MQ=255)
                                                                                                          ccTCAATGGTGGTGGTCGGTAATTTATCCGACGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAtt                                                <  1:381117/149‑1 (MQ=255)
                                                                                                               aTGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATc                                         >  1:18888/1‑151 (MQ=255)
                                                                                                                 ggtggtggtCGTTAATTTACCCGACGGGGCTTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCCCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctg                                       <  3:259382/151‑1 (MQ=255)
                                                                                                                 ggtggtggtCGGTAATTTACCCGAGGGTGCTTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGct                                             <  3:307636/145‑2 (MQ=255)
                                                                                                                 ggtggtggtCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATct                                        <  3:139506/150‑1 (MQ=255)
                                                                                                                    ggtggtGGGTACCTTATCCGACGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctgct                                     <  1:528370/150‑1 (MQ=255)
                                                                                                                    ggtggtCGGTAACTTATCCGACGGTGTGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctgct                                     <  4:426171/150‑1 (MQ=255)
                                                                                                                       ggtCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGAtt                                 >  1:583946/1‑151 (MQ=255)
                                                                                                                        gtCGTTAACTTCCCCGCGGTTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCGTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGAtt                                 <  2:554482/150‑1 (MQ=255)
                                                                                                                           ggTAATTTACCCGACGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGt                              <  2:1256829/150‑1 (MQ=255)
                                                                                                                              aaCTTACCGGGCGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTaaa                           <  1:930094/150‑1 (MQ=255)
                                                                                                                              aaCTTACCCGCCGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTgg                               <  1:448015/146‑1 (MQ=255)
                                                                                                                               aCTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAAt                         >  1:468913/1‑151 (MQ=255)
                                                                                                                                 ttctCCGACGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATc                        <  2:1226282/147‑1 (MQ=255)
                                                                                                                                 ttACCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTaaaa                          <  1:308863/148‑1 (MQ=255)
                                                                                                                                 ttACCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTaaaa                          <  1:614808/148‑1 (MQ=255)
                                                                                                                                    tCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGt                     >  3:281782/1‑150 (MQ=255)
                                                                                                                                        aCGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAacac                >  1:588686/1‑151 (MQ=255)
                                                                                                                                                ttaaaacCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctg         <  2:470263/147‑1 (MQ=255)
                                                                                                                                                 ttAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctggg       >  3:331003/1‑151 (MQ=255)
                                                                                                                                                  tAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggt      >  4:195/1‑151 (MQ=255)
                                                                                                                                                   aaacCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtc     <  2:933529/151‑1 (MQ=255)
                                                                                                                                                      accaccGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGctgggtctgg  >  2:799547/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
TAAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAACCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCACGCTGGGTCTGG  >  NC_000913/801969‑802268

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