Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,027,2330TA29.5% 73.1 / 11.6 44intergenic (+5/+376)fre/fadANAD(P)H‑flavin reductase/3‑ketoacyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base T (20/11);  new base A (1/12);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 7.03e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAGGCGCTTTCGTTGAAGCATCCTGGTCTGCAAGTGGTGCCGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGG‑CGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCG‑AAAA‑TTGCCCGCGATCTGTTTTGCAGTG‑AGCGTAATGCGCGGGAAGATCGCCTGTTT‑GGCGA‑TG‑CGTTTGCATTTATC‑TGAGATATAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGAT  >  NC_000913/4026983‑4027480
                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        
gAGGCGCTTTCGTTGAAGCATCCTGGTCTGCAAGTGGTGCCGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGG‑CGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCG‑AAAA‑TTGCCCGCGATCTGTTTTGCAGTG‑AGCGTAATGCGCGGGAAGATCGCCTGTTT‑GGCGA‑TG‑CGTTTGCATTTATC‑TGAGatat                                                                                                                                                                                                                                                          <  2:70722/251‑1 (MQ=255)
                         gtgggaaaGTGGTGCGGGTGGTGAAACACCCGGAAGCGGTCTGG‑GTGGGGCGTACTGGCCCGGTTTTACGGGGGGTATTGCAGGATCGCGTAACGCTGGCAGAGAATGTTATCTATATTGCCGGACGTTTTGAGATGGCG‑AAAA‑TTGCCCGTTTTCTTTTTGGCAGTT‑AGCGTAAGGCGCGGAAAGACCGCCTTTTT‑GGGGA‑TG‑CGTGTGCTTTTAAC‑TGAAATAAAAAAAACCCCGCCCCTGACAGGCgg                                                                                                                                                                                                                                 <  2:180343/245‑1 (MQ=255)
                                     ggcgggggggtaaacaccgggaggCGGGGTGGCGGGGG‑CGTACGGGCAGCGTGTTAAGGGCGGTATTGAAGGATCACGGTACGCTGGCAGAGCATAATATCTATATGGCGGGACGTTTTGAGAGGGCA‑AAAA‑TGCCCCGCGATTTTTTTGGAAGTG‑AGCGTAATGCGGGGAAAGATCGCCTGTTT‑GGCGA‑TG‑CGTTTGAATTTATC‑TGAGATAAAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCaa                                                                                                                                                                                                                    <  2:266805/228‑1 (MQ=255)
                                        cGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGG‑CGTACTGGCACCGTGTTAACGGCGGTAATTTAGGATAACGGGAAGCTGGCAGAGAAAGATATATTTATTTGCGGGCGTTTTGTGATGGGGGAAAA‑TGGCCC‑CGGTCTGTTTTTCAGGG‑AGGGGAAAGCGCGGGGAGAACGCCTGGTT‑GGGGC‑TGTGGTTTG‑ATTTTTT‑TGAGATAAAAAAAAAACCCCCCCCGACAGGCggggaaaaaaggaaac                                                                                                                                                                                                                   >  2:31317/1‑237 (MQ=255)
                                                       cgaaagcgggtggggggggg‑CGTCCTGCCACCTTGTTACCGGGGGTATTGCAGGACCACGGTAGGCTGGCAGAGAATTATATCTATTTGGCCGGAGTTTTTGAGAGGGGA‑AAAA‑TTGCCCGCGTTCTGTTTGGCAGTG‑AGCGTAATGCGCGGGAAGATCGCTTGTTG‑GGCGG‑TG‑CGTTGGCTTTTATC‑TGAGATAAAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTg                                                                                                                                                                                                  <  1:247319/234‑1 (MQ=255)
                                                          aaGCGGGCTGGCGTGGG‑CGTACTGGCACCGTGTTAACGGCGGTATTGCAGGGTCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCG‑GAAA‑TTGCCCGCGATCTGTTTTGCAGTG‑AGCGGAATGCGCGGGAAGATCGCCTGTTT‑TGCGA‑TG‑CGTTTGCATTTATC‑TGAGATATAAAAAAACCCGCCCCTGACAGGCGGG‑GAGAAAGGCAACCAAACAGTTTTTCAGTggg                                                                                                                                                                                                >  1:284700/1‑249 (MQ=255)
                                                             ggggttgggtgggg‑CTTACGGGCCCCGTGTTACCGGCGGTATGGCAGGATCACGGTACGCTGGCAGAGCATGATATTTATATTGCGGGACGTTTTGAGATGGCG‑AAAA‑TTGCCCGCTATTTTTTTTGCAGTG‑AGCTTAATGGGTGGGAAGATCGCTGGTTG‑GGCTA‑TG‑TGTTTGCTTTTATC‑TGAGATATAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCAttt                                                                                                                                                                                            <  1:360369/240‑1 (MQ=255)
                                                               ggCTGGCGTGGG‑CGTACTGGCACCGTGTTAACGGCGGTATTGCAGGGTCACGGTACGCTGGCAGAGCATGATATCTATATTTCCGGACGTTTTGAGATGGCG‑GAAA‑TTGCCCCCGATCTGTTTTGCAGTG‑AGCGTAATGCGCGGGAAGATCGCCTGTTT‑TGCGA‑TG‑CGTTTGCATTTTTC‑TGAGATATAAAAAAAACCCCCCCTGACAGGCGGG‑GAGAACGGCAAATaaaaatttttttcgggggatttt                                                                                                                                                                                           >  1:268116/1‑230 (MQ=255)
                                                                    ggcgggggcgagcTGGCACCGTGTTAACGGCGGTATTGAGGGATCACGGTACGCTGGCAGAGCATGATATCTTTATTGCCGGACTTTTTGAGATGGCG‑AAAA‑TTGCCCGCGATTTGTTTGGCAGTG‑GGCGAAATGCGGGGGAAGACCGCCTGTTT‑GGCGA‑TG‑GGTTTGAATTGATC‑GGAGATATAAAAAACCCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATGCAGTGGCATTTAGATCt                                                                                                                                                                                      <  1:374880/239‑1 (MQ=255)
                                                                         gggcgtccgggcccggtgtTAACGGTGGTAGTGCAGGATCACGGTCCGCTGGCAGAGCATGATATCTATATTGCCGGACTTTTTGAGATGGGG‑AAAA‑TTCCCCGCGTTCTGTTTTGCAGTG‑AGCGTAATGCGCGGAAAGATCGCCTGTTG‑GGCGA‑TG‑CGTTTGCTTTTATT‑TGGGAGATAAAAAACCCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGAc                                                                                                                                                                                 <  1:55604/236‑1 (MQ=255)
                                                                          g‑CGTACTGGAACCGTGTTAACGGCGGTATTTCCGGATCACGGTAAGCTCGCAGAGCATGAAATCTATATTGAGCGACGTTTTGAAATAGCG‑AAAA‑TTGACCGCGATCTGTTTTGTAGGG‑GGAGTAATGCGCGGGGAAATCGCCTGTTT‑GGGGG‑TG‑GGTTTGCATTTATC‑TGAGATATAAAAAAACCCGCCCCAGAACGGCGGGGAAAAACGGCAA‑TAAAATATTTTTCAGGGGCAtttttatctttaacgg                                                                                                                                                                               >  1:362568/1‑238 (MQ=255)
                                                                                  gtctccgtgttacgggggatATTGAAGGATCCCGGTCCGCGGGAAGGGTATGATATCTATTTTGCCGGGCTTTTTGAGTTGG‑G‑AAAAGTTGCCCGCGATCTGTTGTGCGGTG‑GGCGTAATGCGCGGGAAGTTCGCTTTTTT‑GGCGA‑TG‑CGTTTGCATTTTTCGTG‑GATAAAAAAACTCCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGGATCTGGc                                                                                                                                                                        <  1:57483/232‑1 (MQ=255)
                                                                                  ggCACCGTGTAAGCGGCGTTATTGCGGGTGCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGAGGTTTGGAGATGGCA‑AAAA‑TTGCCCGCGATCTGTTTTGCAGTG‑GGCGTAATGCGCGGGAAGATCGCTTGTTT‑GGCGA‑TG‑CGTTTGAATTTATC‑TGAGATATAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGc                                                                                                                                                                        <  1:11206/251‑1 (MQ=255)
                                                                                       cGTGTTAACGGGGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATAGTGCCGGACGTTTTGAGATGGCG‑AAAA‑TTGCCCGCGATCTGTTTTGCAGTG‑AGCGTAATGCGGGGGAAGACCGCTTGTTT‑GGCGA‑GG‑CGTTGGCATTTATC‑GGAGATAAAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAAcc                                                                                                                                                                   <  1:201472/251‑1 (MQ=255)
                                                                                                        gggaggatacgggTAGGCTGGAAGAGCATGATATTTATATTGCCGGACGTTTTGAGATGGCG‑AAAA‑TTGCCCGAGATTTTTTTTGAGGTG‑AGGTTAATGCGCGGGAAGACCGCCTGTTT‑GGCAA‑TG‑GGTTGGCATTTATC‑TGAAATATAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCt                                                                                                                                                  <  2:417440/240‑1 (MQ=255)
                                                                                                                   ggTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCG‑AAAA‑TTGCCCGCGATCTGTTTTGCAGTG‑AGCGTAATGCGCGGGAAGATCGCCTGTTT‑GGCGA‑TG‑CGTTTGCATTTATC‑TGAGATATAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCtt                                                                                                                                        >  1:263248/1‑250 (MQ=255)
                                                                                                                    gTACGCTGGAAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCG‑AAAA‑TTGCCCGCGATCTTTTTTGCAGTG‑AGCGTAATGCGCGGGAAGATCGCCTGTTT‑GGCGA‑TG‑CGTTTGCATTTATC‑TGAGATATAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTAt                                                                                                                                      <  1:92474/251‑1 (MQ=255)
                                                                                                                               gagCATGATATCTATATTGCCGGACGTTTTGAGATGGCG‑AAAA‑TTGCCCGCGATCTGTTTTGCAGTG‑AGCGTAATGCGCGGGAAGATCGCCTGTTT‑GGCGA‑TG‑CGTTTGCATTTATC‑TGAGATATAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACa                                                                                                                            >  1:113499/1‑250 (MQ=255)
                                                                                                                                  cATGATATCTATATTGCCGGACGTTTTGAGATGGCG‑AAAA‑TTGCCCGCGATCTGTTTTGCAGTG‑AGCGTAATGCGCGGGAAGATCGCCTGTTT‑GGCGA‑TG‑CGTTTGCATTTATC‑TGAGATATAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAAt                                                                                                                         >  1:68709/1‑250 (MQ=255)
                                                                                                                                          cTATATTGCCGGACGTTTTGAGATGGCG‑AAAA‑TTGCCCGCGATCTGTTTTGCAGTG‑AGCGTAATGCGCGGGAAGATCGCCTGTTT‑GGCGA‑TG‑CGTTTGCATTTATC‑TGAGATATAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGCATCCGGCAAAACATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAAc                                                                                                                >  1:59635/1‑251 (MQ=255)
                                                                                                                                            atatTGCGGGACGTTTTGAGATGGGG‑AAAA‑TTGCCCCGGATCTGTTTGGCAGTGGAG‑GTAATGCGCGGGAAGATCGC‑TGTTTTGGCGA‑TG‑CGTTTGCATTTATC‑TGAGATATAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAAcc                                                                                                               <  1:135834/250‑1 (MQ=255)
                                                                                                                                            atatTGCCGGACGTTTGGAGAGGGCA‑AAAA‑TGGCCCGCGTTTTGTTTTGCAGTG‑AGCTTATTGCGGGGGAAGATCGCCTGTTT‑GGCGA‑TG‑CGTTTGCTTTTATC‑TGAGATAAAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAAcc                                                                                                               <  1:236531/250‑1 (MQ=255)
                                                                                                                                              atggccggtgtttttGAGAGGGCG‑AAAA‑TTGCCCGAGTTGTGTTTTGCCGTG‑AGCGTATGGCGCGGGAAAATCGTCTGTTT‑GGCGA‑GG‑CGTTTGCTTTAATG‑TGAGATATAAAAAACCCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTa                                                                                                            <  1:103388/240‑1 (MQ=255)
                                                                                                                                                 gCCGGACGTTTTGAGATGGCG‑AAAA‑TTGCCCGCGATCTGTTTTGCAGTG‑AGCGTAATGCGCGGGAAGATCGCCTGTTT‑GGCGA‑TG‑CGTTTGCATTTATC‑TGAGATATAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAAATAAACTGTTATTCCGAGGCATTTAGATCTATGACGTAACTGGCAAACCATGCCCGATGCGAACCTGTCGCGTCTTATTGTTCCTCCAAATACTCCAAACCCTAGGc                                                                                                         >  1:412691/1‑251 (MQ=255)
                                                                                                                                                   gggACGTTTTGGGATGGCG‑AAAA‑TGGCCCGCAATCTGTTTTGCAGGG‑GGTGTAAGGCGGGGGAAGATCCCCTGTTG‑GGCGA‑TG‑CGTTTGCATTTATC‑TAAGATAAAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCg                                                                                                       <  2:268116/250‑1 (MQ=38)
                                                                                                                                                     gACGTTTTGAGATGGCG‑AAAA‑TTGCCCGCGATCTGTTTTGCAGTG‑AGCGTAATGCGCGGGAAGATCGCCTGTTT‑GGCGA‑TG‑CGTTTGCATTTATC‑TGAGATATAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCCACTTAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACCAATAGTCCGAACCGGAGGcccgg                                                                                                     >  2:59127/1‑248 (MQ=255)
                                                                                                                                                      aCGTTTTGAGATGGCG‑AAAA‑TTGCCCGCGATCTGTTTTGCAGTG‑AGCGTAATGCGCGGGAAGATCGCCTGTTT‑GGCGA‑TG‑CGTTTGCATTTATC‑TGAGATATAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTAAAAATAGTCCGAACCGTAGGCCggg                                                                                                     >  2:67860/1‑249 (MQ=255)
                                                                                                                                                                    ggtgaaaaaggCCCGAGATCTGTTTTGCTGTG‑AGGGTAAGGCGGGGGAAGAGCGACTTTTG‑GGCGA‑TG‑CTTTTGATTTTATC‑TGAGAAATAAAAAAGCCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCGACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTAcg                                                                                        <  1:189945/241‑1 (MQ=25)
                                                                                                                                                                         aa‑TTGCCCGCGATCTGTTTTGCAGTG‑AGCGTAATGCGCGGGAAGATCGCCTGTTT‑GGCGA‑TG‑CGTTTGCATTTATC‑TGAGATATAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAAAAGTCCGAACCGGAGGCCGGGAAAGGCGTTTACGCCGCAt                                                                                  >  2:241630/1‑251 (MQ=255)
                                                                                                                                                                              gCCCGCGATCTTTTTTGCAGTG‑AGCGTAATGCGCGGGAAGATCGCTTGTTT‑GGCGA‑TG‑TGTTTGCATTTATC‑TGAGATATAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                                                                               <  1:387813/250‑1 (MQ=255)
                                                                                                                                                                                  gcgATCTGTTTTGCAGTG‑AGCGTAATGCGCGGGAAGATCGCCTGTTT‑GGCGA‑TG‑CGTTTGCATTTATC‑TGAGATATAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAAAGGCAACTAAACTGTTAATCAGTGGCATTTAGATCTATGAAGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGCGCCTACAAAAAGTCCGAACCGTAGGCCGGATAAGGCGGTTACGCCGCCTCCGGCaac                                                                          >  1:400072/1‑250 (MQ=255)
                                                                                                                                                                                    gATCTGTTTTGCAGTG‑AGCGTAATGCGCGGGAAGATCGCCTGTTT‑GGCGA‑TG‑CGTTTGCATTTATC‑TGAGATATAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATTTATGACGTATCTGGCAAACCATGCCCGCTGCGACGCTGTCGCGTCTTATCGTGGCTACAAATAGACCGAACCGGAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTg                                                                        >  1:317877/1‑251 (MQ=255)
                                                                                                                                                                                                   g‑AGCGTAATGCGCGGGAAGATCCCCTGTTG‑GGCGAGTG‑‑GTTTGTATTTATC‑TGAGATATAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGa                                                         <  2:284700/251‑1 (MQ=255)
                                                                                                                                                                                                          aaTGCGCGGAAAGATCGCTTTTTT‑GCCTA‑TG‑CTTTTGCATTTACC‑TGAAAAAAAAAAAACCCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGAc                                                   <  2:68709/251‑1 (MQ=25)
                                                                                                                                                                                                           aggcgcGGGAAGATCGCTTGTTT‑GGCGT‑TG‑CTTTTGCTTTTATC‑GGAGATAAAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACg                                                  <  1:157689/249‑1 (MQ=25)
                                                                                                                                                                                                                   gAAGATCGCCTGTTT‑GGCGA‑TG‑CGTTTGCATTTATC‑TGAGATATAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTcgcg                                          >  1:288807/1‑251 (MQ=255)
                                                                                                                                                                                                                       aTCGCCTGTTT‑GGCGA‑TG‑CGTTTGCATTTATC‑TGCGATATAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCACTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGACCCGTAGGCCGGGTAAGGCGTTTACTCCCCATCCGGCAATTGGTGCATGCTGCCTGATGCGACGCTGTCGCGTCtt                                      >  2:337418/1‑251 (MQ=37)
                                                                                                                                                                                                                            cTGTTT‑GGCGA‑TG‑CGTTTGCATTTATC‑TGAGATATAAAAAAACCCGCCCCTGACAGGCGGG‑AAGGACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGGTTACGCCGCATCCCGCCATTGGTGCATGATGCCTGATGCGACGCCGTCGCGTCTTTTCGt                                 >  2:131473/1‑251 (MQ=255)
                                                                                                                                                                                                                                tt‑GGCGA‑TG‑GGTTTGCTTTTATC‑TGAGATAAAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGcc                              <  1:303689/250‑1 (MQ=255)
                                                                                                                                                                                                                                   ggCGA‑TG‑GTTTTGCTTTTATC‑TGAAATAAAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTAc                           <  1:311235/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gA‑TG‑CGCTTGCATTTATC‑TGAGATCTAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGGACCCTAACTGTTATTCAGGGGAATTTAGATCTATGACGTATCTGGCAAACCCTGACCCATGCGGCTCCGCCGCGTCTTATCCTCCCCTCAAATAATCCGAACCGCAAGCCGGCCAAGGCGTTTTCGCCGCAACCGGCATTTGGGGCCTCATTCCTGATGCGGCGCTCCCGCGCCCTGTCGTGCCTACaa                         >  2:338060/1‑250 (MQ=255)
                                                                                                                                                                                                                                          g‑CGTGTGCTTTTCTC‑TGAGATAAAAAAAACCCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTGACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATa                      <  1:366698/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 gCATTTATC‑TGAGATATAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGa                >  1:64105/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     ttATC‑TGAGATATAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGGGCATGATGCCTGATGCGACGCCGTCGCGTATTATCGTGCCTACAAATAGTCCGAACCGt           >  1:42907/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       agtcgggagaataAAAAACCCCGCCCCAGACAGGCGAG‑AAGATCGGCAACTAAACTGTTATTAAGTGGAAGTGAGTTCTATGAAGATTCGGGCAAACCAGGCCCGATTCGACGATGTCGCGTCTTATCGTGCCGACAAATATTCCGAACAGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTTCAAATAGTCCGAACCGt           <  2:251573/240‑1 (MQ=18)
                                                                                                                                                                                                                                                            gagaTATAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTc                                             >  1:429657/1‑211 (MQ=255)
                                                                                                                                                                                                                                                               atatAAAAAACCCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGAt  <  1:91014/251‑1 (MQ=255)
                                                                                                                                                                                                                                                                tataAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAAACGTAGGCCGGATAAGGCGTTTTCGCCGCATCCGGCAATTGGGGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGGCCGAACCGtcgggcgg    >  2:398829/1‑241 (MQ=255)
                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        
GAGGCGCTTTCGTTGAAGCATCCTGGTCTGCAAGTGGTGCCGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGG‑CGTACTGGCACCGTGTTAACGGCGGTATTGCAGGATCACGGTACGCTGGCAGAGCATGATATCTATATTGCCGGACGTTTTGAGATGGCG‑AAAA‑TTGCCCGCGATCTGTTTTGCAGTG‑AGCGTAATGCGCGGGAAGATCGCCTGTTT‑GGCGA‑TG‑CGTTTGCATTTATC‑TGAGATATAAAAAAACCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAGGCCGGAT  >  NC_000913/4026983‑4027480

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