Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,1450CG12.9% 116.4 / 7.1 46intergenic (+128/‑70)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (33/7);  new base G (0/6);  total (33/13)
Fisher's exact test for biased strand distribution p-value = 1.83e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
Rejected: Not supported by required number of reads on both strands.

CAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGC  >  NC_000913/3780847‑3781444
                                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                              
cAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCATCTCCCTGAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCgg                                                                                                                                                                                                                                                                                                              <  2:56085/301‑1 (MQ=255)
 aGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCggg                                                                                                                                                                                                                                                                                                             <  2:76756/301‑1 (MQ=255)
         ctttcccttcgcgccgctttttcgatctaaGTGATGAAAGTGGCGATGACGCTGACTGGCTCGAAAGCGATCAGCGAAAGTACGCAAGAGTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGAgggg                                                                                                                                                                                                                                                                                                    <  1:627056/273‑1 (MQ=255)
            gCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGCGGGCCGGGG‑GACGGGG‑aaa                                                                                                                                                                                                                                                                                                >  1:22450/1‑301 (MQ=255)
                        aaTCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCC‑CCGGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACCGc                                                                                                                                                                                                                                                                                    >  1:354869/1‑301 (MQ=255)
                                 gAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATTAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGACTGCGGCATAGTCGTTTGAGCCCCCTCACCCTACCCCTCTCCCTCAGCGAGAGGGGACC‑GTTCGGACCCGTATGTACTCCATAAACATG‑AAACGAAAACGCACTCTTATTTTCTCCCTC‑GCCACA‑CCGGGCAGAGGCCCGGGG‑TAAGGCG‑AACAGGCCGccacacccacccatt                                                                                                                                                                                                                                                                           >  1:536552/1‑287 (MQ=255)
                                 gAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACCGCACACCAAtt                                                                                                                                                                                                                                                                           >  1:358793/1‑301 (MQ=255)
                                     aaGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGCGAGGGCCGGGG‑TGAGGGG‑agatcggaagagcaaacgtctgaactcc                                                                                                                                                                                                                                                                       >  2:627933/1‑274 (MQ=255)
                                           tGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTCCTCCCTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCACACTGCTCCCCAAATTTCTGCTATc                                                                                                                                                                                                                                                                 >  1:227978/1‑301 (MQ=255)
                                                                      cgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGACAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCCCTCCCTCAGGGAGAGGGGACC‑GTTAGGCGCCGTATGTACTCCCTCACTCTG‑ACACGACACCGCCCTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGCGCGACGGACGGGA‑CGAGGGC‑AAAACACCGCACAGCTCACCAATACTCAGCTATCCTCACCACACACCAACGGGTCCGTAtc                                                                                                                                                                                                                                      >  1:28650/1‑300 (MQ=255)
                                                                                    gAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGATGTATGTACTCCCTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGCGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGATAAACATCGTACt                                                                                                                                                                                                                        >  2:414257/1‑301 (MQ=255)
                                                                                                    cGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAAtt                                                                                                                                                                                                        >  2:455348/1‑301 (MQ=255)
                                                                                                          tGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCACGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAAAAATTCCGccc                                                                                                                                                                                                  >  1:34286/1‑300 (MQ=255)
                                                                                                                  gggggttggtggtggcgtggggCACGGGCTCCCATGGGGAAAGGGAGTGGGGAATGGTTGTTTGCCCCCCCGCCCCCCACCCCTCCCCCGCGGGGGGGGGGGTCT‑GTGGGGCGCGTTGTGTCCTCCCTCACCCGG‑AAACGACACCGCCCTTTTTTTGTCCCCCCC‑GCCCCT‑CGGGGGGGGGGGCGGGGG‑GGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTgg                                                                                                                                                                                          <  1:454784/281‑1 (MQ=255)
                                                                                                                    gggTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGTGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCCCTAAGGAGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCAAAAAACAGGcc                                                                                                                                                                                        >  1:446171/1‑300 (MQ=255)
                                                                                                                                        ttttcggtcgtgtgctttccctcttgggggCATGTTCTTTTGCCCCCCCCCCCCCT‑CCCCCCCCCCCCGGGGGGGGGGGACCTGTTCGGCTGTGTTTTACCCCCCCCACTTGA‑AAACAACACCCCCCTTTTTTTTTCCCCCCC‑GCCCCC‑CGGGGGAG‑GGGCCGGGG‑GGGGGGGAAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCt                                                                                                                                                                   <  2:314454/273‑1 (MQ=255)
                                                                                                                                            ggCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGACACATCACCCGC‑CCTTTCCCTAAtt                                                                                                                                                               >  1:367161/1‑300 (MQ=255)
                                                                                                                                             gCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCTAATaa                                                                                                                                                              >  2:417852/1‑300 (MQ=255)
                                                                                                                                              cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTG‑AACCGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGAGGGG‑AAACGGCCGCACTGCCCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGCGGCAGTACGCTCAACATCGCACTTAACGAACCAGAAATTCCGCCAAATACCGCCCACATCATCCGC‑TTTTGCGCTCATAcc                                                                                                                                                             >  1:10676/1‑301 (MQ=255)
                                                                                                                                                       cGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGCGGGG‑ACAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCCCACTACGGGGGCAGTATGCTAAACATCGTCCTTTCCGCACCAGCAATTCACCCACATACTGGCAACATCATCCAT‑CTTTGCGCTACTCCCGTCTTTCgc                                                                                                                                                    >  1:323274/1‑300 (MQ=255)
                                                                                                                                                            ggtattgtgatttggtgttttGCCCCCCCTCCCCCTACCCCTCTCCCTCGGGGAGGGGGCCCC‑TTTCGGCTCTTTATTTCCCCCCACCCTCTA‑AAACGACACCGCACTTTTTTTTTCCCCCCC‑GCCCCC‑CGGGGGAGGGGGCGGGGG‑GGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCTAATACCGGCTTTCGTCTGCa                                                                                                                                               <  2:355353/283‑1 (MQ=255)
                                                                                                                                                              ctgttttgcgggtgtgtgtgctcccctcgccccTACCCCTTCCCCTCGGGGGGGGGGGTCC‑TTGGGGGGCTGTATGTCCTCCCTCACTCTG‑AAGCGACCCCGAACTCTTTTTTGCCCCCCG‑GCCCCT‑CCGGGGGGAGGGCCGGGG‑GGAGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCTAATACCGGCTTTCGTCTGCata                                                                                                                                             <  1:430186/271‑1 (MQ=255)
                                                                                                                                                                  gCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCTagatcggaagagcgtcgtacagggaa                                                                                                                                         >  1:314050/1‑276 (MQ=255)
                                                                                                                                                                            cGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTACGGGGGCAGTATGCTAAACATCGTACTTTCCGAACACGAAATTCCGCCAAATACTGGCACCATCATCCGT‑ATTTGCGCCACTACCGGCCTTCGTCTGCATttcgtagcagcctcg                                                                                                                               >  1:2979/1‑289 (MQ=255)
                                                                                                                                                                             gTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCTAATACCGGCTTTCGTCTGCataagataggaagagcaa                                                                                                                              >  2:430750/1‑286 (MQ=255)
                                                                                                                                                                                    cccccTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCCGAAATTCCGCCAAATACTGGCAACATCCACCGT‑CTTTGCGCTACTACCAACTTTCCCCTGCATCCCACCAAACCGAAGGGCTTCGc                                                                                                                       >  1:370941/1‑301 (MQ=255)
                                                                                                                                                                                            cccTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTG‑AAACGACACTGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCCAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACg                                                                                                               >  2:445705/1‑301 (MQ=255)
                                                                                                                                                                                                  cccTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAgc                                                                                                         >  2:208150/1‑301 (MQ=255)
                                                                                                                                                                                                     tctcCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATggagatcggaagagcacacgtctgaa                                                                                                      >  2:107862/1‑277 (MQ=255)
                                                                                                                                                                                                                  gagGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGAGGGG‑ACAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGAAAAGCGCCTGCGCCGCGCggg                                                                                         >  1:75616/1‑301 (MQ=255)
                                                                                                                                                                                                                   gggggACT‑GTTCGGCGTTTTTTGTCCTCCCTC‑CTTTGTAACCGACACCGCCCTTTTTTTTTCCCCCTC‑GCCCCG‑CCGGGGAGAGGGCGGGGG‑TGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCgggg                                                                                        <  2:287917/300‑1 (MQ=255)
                                                                                                                                                                                                                    ggggACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATCCTGGCAACATCATCCGT‑CTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCTCCGCGCgggac                                                                                       >  2:350142/1‑299 (MQ=255)
                                                                                                                                                                                                                               ggtctttttatgtaCTCCCTCA‑TCTGAAAACGCCCCCGCCCTTTTTTTTTCTCCCTC‑GCCCCT‑CGGGGGGGAGGCCGGGGG‑TGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACTATCa                                                                             <  2:367645/293‑1 (MQ=255)
                                                                                                                                                                                                                                             cTCCCTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCCCTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATCCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCAAATACAGGCTTTCGTATGCATATCATCGACCCGATGGGATTTGCCTGGGACAACAAGCGCATCCGCCGCGCGGGGCGGGACTATCCCGAGTTAAcagccc                                                               >  1:617341/1‑300 (MQ=255)
                                                                                                                                                                                                                                              tccctttattctaaaccgcccccgccttttttttttCTCCCTC‑GCCCTC‑CGGGGGAGGGGGCGGGGT‑GGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACTATCACGAGTTTAccgccg                                                               <  2:34337/272‑1 (MQ=255)
                                                                                                                                                                                                                                               cccTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACTATCACGAGTTTACCGCCGtt                                                             >  2:606131/1‑301 (MQ=255)
                                                                                                                                                                                                                                                ccTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCTAATACCGGCTTTCGTCTGCATATCACCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACTATCACGAGTTTACCGCCGTTa                                                            >  2:210063/1‑301 (MQ=255)
                                                                                                                                                                                                                                                 cctcgcgtgtgaagtgatcgtgatgtgttttttCTCCCTC‑CCCCCT‑CCGGGGGGGGGGGCGGGGTTGGGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACTATCACGAGTTTACCGCCGtt                                                             <  1:417289/274‑1 (MQ=255)
                                                                                                                                                                                                                                                       ctG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAAGGGGCCGGGG‑TGAGGGG‑AAAAGCCCGCCCTGCCCACCATTTTCTTCCTATCCTGCCCCCGCCCTAAGGGGGCAGAATGCTAACCATCGACCTTTCGGAACAAAAAATTCCGCCAAAT‑CCGGCCACACCATTCGTCCTTTCCCCTAATACCGCCTTCCGTCTCCATAAAATCAAACCGACGGCCTTTGCCTGGGACGACAAGCCCCCGCCCCCCCCGGGCCTGGCCTATCACCCCTTTACCACCtttccccctcc                                                     >  1:321696/1‑293 (MQ=255)
                                                                                                                                                                                                                                                        tG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCTAATACCGGCTATCGTCTGCATATCATCGAACCGATGGGATTTCCCTGGGACGCTAAGCGCCTGCGCAGCGCGGGGCTGGACTCTCACGAGTTTACCCCCCTTACGCGtcat                                                    >  1:589628/1‑301 (MQ=255)
                                                                                                                                                                                                                                                         g‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACTATCACGAGTTTACCGCCGTTAAGAGtcatc                                                   >  2:290944/1‑301 (MQ=255)
                                                                                                                                                                                                                                                          aacactccccctttctctTTTTTCCTCCCTC‑CCCCCC‑CGGGGGGGGGGGCGGGGG‑GGGGGGA‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGtcatc                                                   <  2:75706/287‑1 (MQ=255)
                                                                                                                                                                                                                                                                        cgttttttttGTCCCCCTGCCCCTGCGGGGGGGGGGGCCGGGT‑GGGGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATcgcg                                       <  1:207907/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                               cccggcccTC‑TCCCCG‑GAGGGGCGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACTATCACGAGTTTACACCAGTTACCCAACATCATAACTACCGCCCCTTCCTcc                              >  1:281139/6‑300 (MQ=255)
                                                                                                                                                                                                                                                                                      cTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCTAATACCGGCTTTCGTCTGCATATCCTCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACTACCACGAGTTTACCGCCGTTACGCGTCATCATGACTAACGCCCGTTCCCCGAAGccaa                       >  1:2152/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                                                  cGGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAAACCCCaa             >  2:179062/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                             ccGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGc  >  2:638252/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                              
CAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACC‑GTTCGGCGCTGTATGTACTCCCTCACTCTG‑AAACGACACCGCACTCTTTTTTTCTCCCTC‑GCCCCT‑CCGGGGAGAGGGCCGGGG‑TGAGGGG‑AAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGT‑CTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGCTGGACTATCACGAGTTTACCGCCGTTACGCGTCATCATGACTATCGCGCGTTCCTCGAAGCAGAAAATCCCCAGCGCCTGTTCGC  >  NC_000913/3780847‑3781444

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