Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,520,7510GT33.3% 47.8 / 14.4 49T90K (ACA→AAA) hofPDNA utilization protein HofP
Reads supporting (aligned to +/- strand):  ref base G (8/24);  new base T (16/0);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 6.52e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGT‑AAAATTGTCCAGCC‑GTTTTCCAG‑CACAT‑C‑GTTTTGCTGCA‑CCCGTCGCC‑ATTTCTTTTGCCCGTCTTTTATTACACCGATGATGCGCTCGCCTCGCCCTACCATCCCCTGATAGCGCCACTGGCTAAGTTCGCTAATCCGGCATAGAT‑CTTCCGGCGG‑TTTAAAAGGGTCACGCATACCGGTTAAAAGGCACAATGCAATACC  >  NC_000913/3520509‑3520992
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                          
aaCATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTAAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCCG‑TTTTTccccta                                                                                                                                                                                                                                                     >  1:283078/1‑247 (MQ=255)
       ggtAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTGCGCAG‑TTTGTCCCGGTACCCAg                                                                                                                                                                                                                                               <  2:39418/247‑1 (MQ=255)
          aGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTGTCCCGGTACCCAGCGtt                                                                                                                                                                                                                                           <  1:110508/250‑1 (MQ=255)
                 aCGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCGTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTGTCCCGGTACCCAGCGTTAGTACGTc                                                                                                                                                                                                                                   <  1:1808/251‑1 (MQ=255)
                      acCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTTTCCCCGTACCCCGCGGTAGGACGTCTGGCg                                                                                                                                                                                                                              >  1:68993/1‑251 (MQ=255)
                       ccGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTGTCCCGTTACCCAGCGTTAGTACGTCTGGCGt                                                                                                                                                                                                                             <  2:324078/251‑1 (MQ=255)
                        cGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTGTCCCGGTACCCAGCGTTAGTACGTCt                                                                                                                                                                                                                                  <  1:105851/245‑1 (MQ=255)
                        cGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTc                                                                                                                                                                                                                            <  1:248966/251‑1 (MQ=255)
                                           gACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTGTCCCGGTAACCAGCGTTAAGACGTCTGGCGTCAGCTGT‑AAAATTGTCCAGc                                                                                                                                                                                                        >  1:234379/1‑251 (MQ=255)
                                             cacGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTTTCCCCGTACCCAGCGTTAGTACTTCTTGCGTCCGCTGTTAAAAtatcacgccg                                                                                                                                                                                                       >  1:47708/1‑242 (MQ=255)
                                                         cAACTTCTCCTGTTCCGCCAGCGCCTGCAACCCCTGAGCTACCGGAACGTCACCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCAGTGTGGCGGTTCGCAG‑TTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGT‑AAAATTGTCCAGCC‑GTTTTCCAG‑caca                                                                                                                                                                                        <  2:68993/251‑1 (MQ=255)
                                                                tcCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCCTCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGCGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAC‑TTTTTCCCCGTACCCAGCGGTCGTACGTCTGGCGTACGCTGT‑AAAATTGTCCCGCC‑GTTTTCCAG‑CACAT‑C‑Gtttt                                                                                                                                                                               >  1:210878/1‑251 (MQ=255)
                                                                            agcGCCTGCACCACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGT‑AAAATTGTCCAGCC‑GTTTTCCAG‑CACAT‑C‑GTTTTGCTGCA‑CCCGTc                                                                                                                                                                  <  1:53004/251‑1 (MQ=255)
                                                                               gcCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGT‑AAAATTGTCCAGCC‑GTTTTCCAG‑CACAT‑C‑GTTTTGCTGCA‑CCCGTCGcc                                                                                                                                                               <  1:293069/251‑1 (MQ=255)
                                                                                 cTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTTTCCCGGGTCCCCGCGTTAGTACGTCTGGGGGCAGCTGT‑AAAAATTTCCCGCC‑GTTTTTCAG‑CACAtaaatcgggagagcgtactgtggg                                                                                                                                                                >  2:85112/1‑228 (MQ=255)
                                                                                 cTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGT‑AAAATTGTCCAGCC‑GTTTTCCAG‑CACAt                                                                                                                                                                                       <  1:85112/228‑1 (MQ=255)
                                                                                     aaCACCTGAGCTACCGGACCGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGACGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCAGTTGCCACAACCATTGTGGCGGTTCGCCG‑TTTGTCCCGGGACCCCGCGTTCGTACGTCTGGGGTCAGATGT‑AAACATGTACCGCC‑GTTTTTCAG‑CACCT‑C‑GTTTTGCTGCA‑CCCgcttccatttat                                                                                                                                                         >  1:265816/1‑240 (MQ=255)
                                                                                        acCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCCGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTCTCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTTGCAG‑TTTTTCCCGGTACCCAGCGCTAGTCCGACTGGGGTCAGCTGT‑GTAAATGTCCACCC‑GGTTTCCCA‑CACCT‑CCGTTTT‑CTGCA‑CCCGTCCCC‑CTTTCtttt                                                                                                                                                     >  1:141873/1‑251 (MQ=255)
                                                                                                tACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTTTCCCGGTACCCAGCGATAGTACGTCTGGCCTCAACTGT‑AAAATTGTCCAGCC‑GTTTTTCAG‑CATAA‑C‑GTTTTTCTTCA‑CACGTCGCC‑AGTTATTTTTTCCGTCt                                                                                                                                             >  1:208098/1‑251 (MQ=255)
                                                                                                                          tCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑GTTTTCCCGGGACCCAGCGGTAGTACGTCTTGCGTCAGCTGT‑AAAATTGGCCCACC‑GTTTTCCCA‑CACAT‑C‑GTTTTTCTGCA‑CCCCTCTCC‑ATTTCTTTTTTCCGGCCTTTATTTCACCCATGATGCGcacgcc                                                                                                                   >  1:65890/1‑251 (MQ=255)
                                                                                                                            ccTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTGTCCCGGTACCCAGCGTTAGTACGGCTGGCGTCAGCTGT‑AAAATTGTCCAGCC‑GTTTTCCAG‑CACAT‑C‑GTTTTGCTGCA‑CCCGTCGCC‑ATTTCTTTTGCCCGTCTTTTATTACACCGATGATGCGctcgcctc                                                                                                                 <  2:123004/250‑1 (MQ=255)
                                                                                                                               ttttGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGT‑AAAATTGTCCAGCC‑GTTTTCCAG‑CACAT‑C‑GTTTTGCTGCA‑CCCGTCGCC‑ATTTCTTTTGCCCGTCTTTTATTACACCGATGATGCGctcgcctcg                                                                                                                <  2:287800/249‑1 (MQ=255)
                                                                                                                                  tGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTGTCCCGGGACCCAGCGTTAGTACGTCTGGCGTCAGCTGT‑AAAATTGTCCAGCC‑GTTTTTCAG‑CACAT‑C‑GTTTTTCTGCA‑CCCGTCGCC‑CTTTTTTTTGTCCGTCTTTTTTTACACCGCTGTTGCGCTCGCCTCGCCCtc                                                                                                           >  1:234110/1‑250 (MQ=255)
                                                                                                                                     ggCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCCG‑TTTTTCCCCGGACCCCGCGTTAGTACGTCTGGCGTCCGCTGT‑TAAATTTTCCAGCC‑GTTTTCCCG‑CCCCT‑C‑GTTTTTCTGCA‑CCCGCCGCC‑CTTTCTTTTTCCCCTCTTTTTTTACCCCGCTGCTGCGCCCGCCTCGCCCTAccc                                                                                                        >  1:156684/1‑250 (MQ=255)
                                                                                                                                         ttGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTTTCCCGGTACCCCGCGTTAGTACGTCTGGCGTCAGCTGT‑AAAATTGTCCAGCC‑GGTTTCCAG‑CACAT‑C‑GTTTTGCTGCA‑CCCGTCTCC‑ATTTTTTTTGCCCGTCTTTTATTTTACCGATAGTGCGCTCCCCTCGCCCTCCCATccc                                                                                                    >  2:94189/1‑251 (MQ=255)
                                                                                                                                                  gTACGCCGGGTACCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGT‑AAAATTGTCCAGCC‑GTTTTCCAG‑CACAT‑C‑GTTTTGCTGCA‑CCCGTCGCC‑ATTTCTTTTGCCCGTCTTTTATTACACCGATGATGCGCTCGCCTCGCCCTACCATCCCCTGATAgc                                                                                            <  1:207773/250‑1 (MQ=255)
                                                                                                                                                     cGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCCCAG‑GTTCTCCCGGGACCCAGCGTTTGTACGTCTTGGGTCCGCTTT‑AAAATTGTCCAGCC‑GTTTTTCAG‑CACCT‑C‑GTTTTGCTGCAACCCG‑CGCCCATTT‑TTTTTCCCGCCTTTTTTTTCCACGCTGCTGCGCTCGCCCCGCCCCAccctacccctatcaggccc                                                                                        >  1:146539/1‑234 (MQ=255)
                                                                                                                                                       ccGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAC‑TTTGTCCCGGTACCCCCCGTTAGTCCGTCTGGCGTCCGCTGG‑AAAATTGTCCCGCC‑GTTTTCCAT‑CCCAT‑C‑GTTTTTCTGCA‑CCCGCCCCC‑CTTTCTTTTTCCCCTCTTTTTTTTCACCCATTCTGCGCTCCCCTCcccccccccaacccttatagcgccccc                                                                                      >  1:45230/1‑228 (MQ=255)
                                                                                                                                                         gggTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGATCGCAG‑TTTTTCCCGGTCCCCAGAGTTAGTACGTCTGGCGCCAGCTGT‑AAAATTTTCCAGCC‑GTTTTCCCG‑CACAT‑C‑GTTTTGCTGCA‑CCCCTCGGC‑CTTTATTTTTCCCGTCTTTTATTACACCGATGATGCGCTCGCCCCGCCCTACCATCCCCTGATTACGCCAccgg                                                                                    >  1:172453/1‑248 (MQ=255)
                                                                                                                                                               cAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGT‑AAAATTGTCCAGCC‑GTTTTCCAG‑CACAT‑C‑GTTTTGCTGCA‑CCCGTCGCC‑ATTTCTTTTGCCCGTCTTTTATTACACCGATGATGCGCTCGCCGCGCCCTACCATCCCCTGATAGCGCCACTGGCTAAGt                                                                              <  2:234379/251‑1 (MQ=255)
                                                                                                                                                                  cATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTTTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGT‑AAAATTGTCCAGCC‑Gtt                                                                                                                                                                                                   >  2:167756/1‑136 (MQ=255)
                                                                                                                                                                  cATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGT‑AAAATTGTCCAGCC‑Gtt                                                                                                                                                                                                   <  1:167756/136‑1 (MQ=255)
                                                                                                                                                                            agtGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTTTCCCGGTACCCAACGTTAGTACGTCTGGCGTCAGCTGT‑AAAATTTTCCAGCC‑GTTTTCCAG‑CACAT‑C‑GCTTTGCTGCA‑CCCGTCGCC‑ATTTCTTTTTCCCGTCTTTTTTTACACCGATGATGCGCTCGCCCCGCCCCCCCATCCCCTGGTAGCGCCACTGGCTAAGtttccctatcccgc                                                                 >  2:299581/1‑239 (MQ=255)
                                                                                                                                                                            agtGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGT‑AAAATGGTCCAGCC‑GTTTTCCAG‑CACAT‑C‑GTTTTGCTGCA‑CCCGTCGCC‑ATTTCTTTTGCCCGTCTTTTATTACACCGATGATGCTCTCGCCTCGCCCTACCATCCCCTGATAGCGCCACTGGCTAAGTTCGCTAAGCCGGc                                                                 <  2:65890/251‑1 (MQ=255)
                                                                                                                                                                                       ccATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGT‑AAAATTGTCCAGCC‑GTTTTCCAG‑CACAT‑C‑GTTTTGCTGCA‑CCCGTCGCC‑ATTTCTTTTGCCCGTCTTTTATTACACCGATGATGCGCTCGCCTCGCCCTACCATCCCCTGATAGCGCCACTGGCTAAGTCCGCTAATCCGGCATAGAT‑CTTcc                                                     <  1:90222/251‑1 (MQ=255)
                                                                                                                                                                                         aTTGCTTCATTTGTATCTCCTTGACGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTGTTCCGGCACCCAGGGTTTGGAAGTCTGGCCTCAGCCGG‑AAAATTTTTCAGCC‑CTTTTTCCG‑CACAT‑C‑GTGTTTCTGCA‑CCCGCCTCC‑CTTTCTTTTTCCCGTCTTTTTTTTCACCCCTGGTACGCCCGCCTCGCCCCTCCCCCCCCTGCTAGCGCCCCCGGCCAAATTCGCTCATcccgcccaaaccttccgg                                                    >  2:40487/1‑235 (MQ=255)
                                                                                                                                                                                           tGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCCGATTTTTCCCG‑TCCCCAACGTTTGTAAGTCTGGCGTCAGCTGT‑AAAATTTTCCCGCC‑GTTTTCCAGACCCA‑‑C‑GTTTTTCTGCA‑CCCGGCGCC‑ATTTTTTTTATCCCTCTTTTATTTAACCCCAGATGCGCTCCCCTCGCCCCACCAACCCCTCATACTGCCGCTCGCTCAGTTTGCTGATCCGGccgaatacttcccgcg                                                  >  2:123105/1‑236 (MQ=255)
                                                                                                                                                                                                cATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑GTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGT‑AAAATTGTCCAGCC‑GTTTTCCAG‑CACAT‑C‑GTTTTGCTGCA‑CCCGTCGCC‑ATTTCTTTTGCCCGTCTTTTATTACACCGATGATGCGCTCGCCTCGCCCTACCATCCCCTGATAGCGCCACTGGCTAAGTTCGCTAATCCGGCATAGAT‑CTTCCGGCGG‑TTTa                                           <  1:41470/251‑1 (MQ=255)
                                                                                                                                                                                                 aTTTGTATCTCCTTGCCGTTGCCACAACCATTGTGTCGGTTAGCAG‑GTTGTCCCGGTACCCAGGGTTAGTAATTTTTGCGTACGCTGT‑AAAAATGTCCAGCC‑GTTTTACAA‑CACAT‑C‑GTTTTGCTGCA‑CCCGTCGCC‑ATTTCTTTTGCCCGTATTTTATTAAACCGATGAGAGGCACGCCTCGCCCTACCCTCCCCCGGTAGCGCCACTGGCCAAGTTCTCTAATCCCGCacatactttccagcggtttt                                             >  1:210825/1‑231 (MQ=255)
                                                                                                                                                                                                    tGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGT‑AAAATTGTCCAGCC‑GTTTTCCAG‑CACAT‑C‑GTTTTGCTGCA‑CCCGTCGCC‑ATTTCTTTTGCCCGTCTTTTATTACACCGATGATGCGCTCGCCTCGCCCTACCATCCCCTGATAGCTCCACTGGCTAAGTTCGCTAATCCGGCATAGAT‑CTTCCGGCGG‑TTTaaaa                                        <  1:198854/250‑1 (MQ=255)
                                                                                                                                                                                                           ccTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTGTCCCGGTGCCCAGCGTTAGTACGTCTGGCGTCAGCTGT‑AAAATTGTCCAGCC‑GTTTTCCAG‑CACAT‑C‑GTTTGGCTGCA‑CCCGTCGCC‑ATTGCTTTTGCCCGTCTTTTATTACACCGATGATGCGCTCGCCTCGCCCTACCATCCCCTGATAGCGCCACGGGCTAAGGTCGCTAATCCGGCATAGAT‑CTGCCGGCGG‑GTTAAAAGGGTCACg                                <  2:210878/251‑1 (MQ=255)
                                                                                                                                                                                                             ttGCCGTTGCCACAACCATTGTGGCGGTTCGCCG‑TTTTTCCCCGGACCCAGCGTTAGTACGTCTTGCGTCAGGTGT‑AAAATTTTCCAGCCCGTTTTC‑AG‑CACATAT‑GTTTT‑CTGCA‑CCCCTCGCC‑ATTTCTTTTTCCCGTCTTTTTTTTCACCGCTGAAGCGCTCGCCTCGCCCTCCCATCCCCTGATAGCGCCACTGGGTTAGTTCTCTCATCCCGCATA‑ATACTTCCGGCGGGTTTATAAGGG‑CACGcc                              >  2:270979/1‑250 (MQ=255)
                                                                                                                                                                                                             ttGCCGTTGCCACAACCATTGTGGCGGGTCGCCG‑TTTGTCCCGGGCACCAGCGGTAGTACGTCCGGCGTCAGCTTT‑AAAAATGTCCAGCC‑GTTTTCCCG‑CCCAT‑C‑GTTTTGCTGCA‑CCCGTCGCC‑ATTTCTTTTTCCCGTCTTTTATTACCCCGATGATGCGCTCCCCCCGCCCTACCCTCCCCTGATAGCGCCCCCGGCTAAGTTCGCTAATCCGGCATAGAT‑CTTCCGGCGGGTTTAAAAggtcccgcc                               >  1:68369/1‑244 (MQ=255)
                                                                                                                                                                                                                    tGCCACAACCATTGTGGCGGTTCGCAG‑TTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGT‑AAAATTGTCCAGCC‑GTTTTCCAG‑CACAT‑C‑GTTTTGCTGCA‑CCCGTCGCC‑ATTTCTTTTGCCCGTCTTTTATTACACCGATGATGCGCTCGCCTCGCCCTACCATCCCCTGATAGCGCCACTGGCTAAGTTCGCTAATCCGGCATAGAG‑TTTCCGGCGG‑TTTAAAAGGGTCACGCATACCgg                        <  1:277529/250‑1 (MQ=255)
                                                                                                                                                                                                                     gCCACAACCATTGTGGCGGTTCGCAG‑TTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGT‑AAAATTGTCCAGCC‑GTTTTCCAG‑CACAT‑C‑GTTTTGCTGCA‑CCCGTCGCC‑ATTTCTTTTGCCCGTCTTTTATTACACCGATGATGCGCTCGCCTCGCCCTACCATCCCCTGATAGCGCCACTGGCTAAGTTCGCTAATCCGGCATAGAT‑CTTCCGGCGG‑TTTAAAAGGGTCACGCATACCGGtt                      <  1:94189/251‑1 (MQ=255)
                                                                                                                                                                                                                     gCCACAACCATTGTGGCGGTTCGCAG‑TTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGT‑AAAATTGTCCAGCC‑GTTTTCCAG‑CACAT‑C‑GTTTTGCTGCA‑CCCGTCGCC‑ATTTCTTTTGCCCGTCTTTTATTACACCGATGATGCGCTCGCCTCGCCCTACCATCCCCTGATAGCGCCACTGGCTAAGTTCGCTAATCCGGCATAGAT‑CTTCCGGCGG‑TTTAAAAGGGTCACGCATACCGGt                       <  1:299581/250‑1 (MQ=255)
                                                                                                                                                                                                                       cacaACCATTGTGGCGGATTGGAG‑TTTTTCCAGGTACCCAGCGATAGAACGTCTGGCGGCCGCAGT‑TAACTTTTCCAGCC‑GTTTTCCAA‑CACAT‑C‑CTTTTGGTGCA‑CCCGTCGCC‑CTTTCTTTTAGCCGGCTTTTATTACAACGATGATGCGCTCTCCTCGCCCTACCATCCCCTGATAGCGCCACTGGGTAAGTTCGCTAATCCGGGCTAGAT‑CATCCGGGGG‑TTTAAAAAGGGCAAGCACACCgggta                     >  2:334191/1‑247 (MQ=255)
                                                                                                                                                                                                                                      ggTTCGCAG‑TTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGT‑AAAATTGTCCAGCC‑GTTTTCCAG‑CACAT‑C‑GTTTTGCTGCA‑CCCGTCGCC‑ATTTCTTTTGCCCGTCTTTTATTACACCGATGATGCGCTCGCCTCGCCCTACCATCCCCTGATAGCGCCACTGGCTAAGTTCGCTAATCCGGCATAGAT‑CTTCCGGCGG‑TTTAAAAGGGTCACGCATACCGGTTAAAAGGCACAATGCAAt     >  1:105261/1‑251 (MQ=255)
                                                                                                                                                                                                                                          cGCAG‑TTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGT‑AAAATTGTCCAGCC‑GTTTTCCAG‑CACAT‑C‑GTTTTGCTGCA‑CCCGTCGCC‑ATTTCTTTTGCCCGTCTTTTATTACACCGATGATGCGCTCGCCTCGCCCTACCATCCCCTGATAGCGCCACTGGCTAAGTTCGCTAATCCGGCATAGAT‑CTTCCGGCGG‑TTTAAAAGGGTCACGCATACCGGTTAAAAGGCACAATGCAATAcc  <  1:40487/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                          
AACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCAGCGTCACTTTTTGCGGCTTTGCCGCCTGTACGCCGGGTATCAGCATCAACAGTAGTGCGGCTATCCATTGCTTCATTTGTATCTCCTTGCCGTTGCCACAACCATTGTGGCGGTTCGCAG‑TTTGTCCCGGTACCCAGCGTTAGTACGTCTGGCGTCAGCTGT‑AAAATTGTCCAGCC‑GTTTTCCAG‑CACAT‑C‑GTTTTGCTGCA‑CCCGTCGCC‑ATTTCTTTTGCCCGTCTTTTATTACACCGATGATGCGCTCGCCTCGCCCTACCATCCCCTGATAGCGCCACTGGCTAAGTTCGCTAATCCGGCATAGAT‑CTTCCGGCGG‑TTTAAAAGGGTCACGCATACCGGTTAAAAGGCACAATGCAATACC  >  NC_000913/3520509‑3520992

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