Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,019,570 A→C 14.5% T138P (ACA→CCA)  ssnA → putative chlorohydrolase/aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,019,5700AC14.5% 173.4 / 2.3 70T138P (ACA→CCA) ssnAputative chlorohydrolase/aminohydrolase
Reads supporting (aligned to +/- strand):  ref base A (4/52);  new base C (9/2);  total (15/55)
Fisher's exact test for biased strand distribution p-value = 8.08e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGGAGATGCTGCGTGAAGCCG  >  NC_000913/3019338‑3019807
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    
cTCGCACAACCATTTTTACTCGGGGCTTTTCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGACGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTCTATCGGCGGT‑TGGCT‑AT‑Ccccacttc                                                                                                                                                                                                                                                >  2:806522/1‑235 (MQ=255)
        aCCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCa                                                                                                                                                                                                                                     <  1:615674/239‑1 (MQ=255)
        aCCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCa                                                                                                                                                                                                                                     <  2:1358502/239‑1 (MQ=255)
        aCCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACAGCGCCCCTTGCCCGGATTTCATCCCAACGCTGAAAAATCTCTGGTGGCGGCTCGAGCGCGCCCTTGATGAAGAGGCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCa                                                                                                                                                                                                                                     <  2:268893/239‑1 (MQ=255)
                 aCTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCC‑CT‑CCCCACATCG‑C‑GA‑CGCC‑TTTTTtaaa                                                                                                                                                                                                                           >  1:572312/1‑239 (MQ=255)
                  cTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAg                                                                                                                                                                                                                          <  1:608539/239‑1 (MQ=255)
                    cGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCAGT‑CC‑CCACATTCCG‑C‑GAGCCCT‑TTTTTTAAAGt                                                                                                                                                                                                                         >  1:147565/1‑239 (MQ=255)
                     ggggCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TGGCT‑CT‑CCCCCTTCCC‑C‑GC‑CCCC‑TTTTTTAAAATTg                                                                                                                                                                                                                       >  1:480220/1‑239 (MQ=255)
                            tCCCGCGGAATTATGGCAAACCTCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTgcg                                                                                                                                                                                                                <  1:645639/239‑1 (MQ=255)
                            tCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CTCCCCCCTTCCC‑C‑GC‑CCCA‑TTTTTAAAAATTGGCCtc                                                                                                                                                                                                                  >  1:336441/1‑237 (MQ=255)
                                gcgGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCTGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGcgcgcg                                                                                                                                                                                                            <  2:643324/239‑1 (MQ=255)
                                     atcatcccAAGCATCGCCCCTTGCCCGTACTTCATCTCAACGCCGAAAAATCTCTGGTGGTGGCTCGTGCCCGCCTCTGATTAAGAGTCGCTCTATTACAGCGGACTGATTTGTACCCTGGCAGCGATTAAGCGCCGATGTTCATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑AT‑CCCCATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGAc                                                                                                                                                                                                       <  2:1020656/232‑1 (MQ=255)
                                     aTTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGAc                                                                                                                                                                                                       <  1:115925/239‑1 (MQ=255)
                                     aTTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGAc                                                                                                                                                                                                       <  2:286083/239‑1 (MQ=255)
                                     aTTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGAc                                                                                                                                                                                                       <  1:1266862/239‑1 (MQ=255)
                                     aTTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGAc                                                                                                                                                                                                       <  1:860337/239‑1 (MQ=255)
                                                cATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTTCGCT‑TT‑CC‑CTCTTCA‑C‑CA‑CGCA‑CTTTTTATATTTTGCCTgcgcgcccgggccgttttttt                                                                                                                                                                                              >  2:953527/1‑222 (MQ=255)
                                                cATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑GCCCT‑CT‑CCCCACTTCC‑C‑GG‑CGCA‑CTTTTTTAAATTGGCCCGCGCGCGGTGACC‑TGTGTTTGaa                                                                                                                                                                                           >  1:145530/1‑239 (MQ=255)
                                                   cGCCCCGTGCGCGGAATCCAGCGCAACGCTGAAAAAGCGCTGGTGGCGGCTCGATCGCGCACTTGAAGAAGAGTCGCTGTAATACAGCGGACTGATTTGTACTCTGGAAGCGAGTAAGAGTGGATGTACATCGGTTATAGATCGCCATGCCTCTCCGGCGTATAACGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAtcta                                                                                                                                                                                       <  2:1340629/239‑1 (MQ=255)
                                                      cccTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTg                                                                                                                                                                                    <  1:22612/239‑1 (MQ=255)
                                                           gCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TGGCC‑CT‑CCCCCCTGCT‑C‑GG‑CCCC‑TTTTTTAAAATCTGCCTGCGCGCGGTGTCCTTGt‑tttatgatgctaacag                                                                                                                                                                                 >  2:805000/1‑225 (MQ=255)
                                                                 aTTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGGGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTaaca‑a                                                                                                                                                                         <  1:969440/238‑1 (MQ=255)
                                                                 aTTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTaaca‑a                                                                                                                                                                         <  1:333550/238‑1 (MQ=255)
                                                                        ctcAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCa                                                                                                                                                                 <  2:556875/239‑1 (MQ=255)
                                                                                 gAAAAATCTCTGGTGGCGGCACGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGGAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TGTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTg                                                                                                                                                         <  1:1293771/238‑1 (MQ=255)
                                                                                     aaTCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGa                                                                                                                                                    <  2:1119850/239‑1 (MQ=255)
                                                                                     aaTCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGa                                                                                                                                                    <  1:589082/239‑1 (MQ=255)
                                                                                     aaTCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGa                                                                                                                                                    <  1:1093766/239‑1 (MQ=255)
                                                                                     aaTCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGa                                                                                                                                                    <  1:430480/239‑1 (MQ=255)
                                                                                     aaTCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGa                                                                                                                                                    <  1:1308940/239‑1 (MQ=255)
                                                                                     aaTCTCTGGGGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGa                                                                                                                                                    <  2:1107534/239‑1 (MQ=255)
                                                                                     aaTCTCTGATGGCTGCTCGATCGCGCCCTTGATGAAGAGTAGCTCTATTACAGCGGACTGATTTGTCCCCTGGAAGCGATTAAGAGCGGATGTCCATCTGTTATCGATCACCATGCCTCTGCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGa                                                                                                                                                    <  2:1068358/239‑1 (MQ=255)
                                                                                     aaTCGCTGGGGGCGGCTCGAGCGCGTCCTTGACGAAGAGTCGCGCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGa                                                                                                                                                    <  1:859148/239‑1 (MQ=255)
                                                                                              tggCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTag                                                                                                                                            <  1:977407/238‑1 (MQ=255)
                                                                                                 cggcTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TGGCT‑CT‑CCCCCCTTCT‑C‑CC‑CTCC‑TTTTTTTAAATTTGGCTGCGCGCGGGGGCC‑CGT‑TTTTTTACTATTTACCCTCAAA‑ACGGCATCACACAGATGTGGGGAGGTGGaggaga                                                                                                                                        >  2:977407/1‑239 (MQ=255)
                                                                                                  ggcTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAgaagaa                                                                                                                                       <  1:1285163/239‑1 (MQ=255)
                                                                                                  ggcTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAgaagaa                                                                                                                                       <  2:550246/239‑1 (MQ=255)
                                                                                                   gcTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TGGGT‑CTCCCCCCTTGCC‑G‑CC‑CCCG‑TTTTTAAAATTTGG‑CTGCGCGCGCAGGCC‑CTT‑TTTTTAACCACTTACCGCAACACACGCCA‑AAAAGAGCTTCAGGAGGGGGTATGAGaaa                                                                                                                                      >  2:1271450/1‑239 (MQ=255)
                                                                                                                             cGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTtgattgat                                                                                                            <  1:1231936/239‑1 (MQ=255)
                                                                                                                                tctATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTaaca‑a                                                                                                                                                                         <  1:381511/175‑1 (MQ=255)
                                                                                                                                tctATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAg                                                                                                         <  2:997020/239‑1 (MQ=255)
                                                                                                                                  tATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGa                                                                                                                                                    <  1:111187/194‑1 (MQ=255)
                                                                                                                                          cGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGGT‑CT‑CCACAATTCG‑CGGC‑CCCA‑TTTTTTAAA‑TTGTGCTTCGGGCGCGGTCC‑CGT‑TTTTTAAATCTTGTCCGTAACA‑AAAGCCTCCAAGAGGTCCAGGAAGGTGTAGAAGAAAAAACCCGTTTTGCCCCTTTGTTtgttgacggggagaca                                                                                                >  2:831626/1‑228 (MQ=255)
                                                                                                                                          cGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGA‑GCGCT‑CT‑CCACATTTCGCC‑GC‑CCCA‑TTTTTAAAA‑TTGGCCTTGGCGCGGGGGAA‑TGT‑TTTGAAACTACTTACAGTAACA‑ACCGCATCAAAGAGATGCAGGAGAGGGGATAAGAAAAAATACGATTCGCCTGttagttttttatggcggggacaa                                                                                               >  2:411750/1‑218 (MQ=255)
                                                                                                                                            gACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TGGCT‑CT‑CCTCACTTCG‑C‑GG‑CGCC‑TTTTTTTAACATGGGGTGCGGGGGGGGACC‑CGT‑TTTTTAACTACTGAACCCAAAA‑ACGGGATCAAAAAGATGCAGGGGGGGGGAGAAAAAAAAAACCTGTTCGCCGCCTTGTTTGTTGAGGGGGAGAAAgag                                                                                             >  2:1316219/1‑237 (MQ=255)
                                                                                                                                                      ttCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGc                                                                                   <  1:1111813/239‑1 (MQ=255)
                                                                                                                                                      ttCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGc                                                                                   <  1:597906/239‑1 (MQ=255)
                                                                                                                                                            ggAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑GCGGC‑TT‑CCCCCCTTCG‑C‑CG‑CGCC‑TTTTTTAAAAATGGCGTGCGGGCGACGACC‑CCT‑TTTGTAACTACTGAACCTAAAA‑AAAGCAGCAAAGAGTTGTGGGAAGGGGGAGAAGAAAAAAACCACTTTGCCCCCTTGATTGTTGGGGGGGAGAAAACGACAAACGAGCCgaa                                                                               >  2:87448/1‑235 (MQ=255)
                                                                                                                                                                  gATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCtggtgg                                                                       <  1:173111/239‑1 (MQ=255)
                                                                                                                                                                   aTTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCtggtgg                                                                       <  2:691083/238‑1 (MQ=255)
                                                                                                                                                                       agagCGGATGCAAACCGGTTATCTAGCACCATGCCTCTCCGGCGTAAATGGGCGGG‑TCGCA‑CG‑CCACATTGCG‑C‑GA‑CGAATTTTTTAAAAGTTGGCCTGGGCGCGAGGACC‑TGT‑TTTGAAACTCCTGACCGTAACA‑ACGCCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGaa                                                                     <  2:147565/237‑1 (MQ=255)
                                                                                                                                                                                 aacagccgtcacgtaTCCCCATGCGTCTCCGGCGTATATAGCCGGG‑TGGCT‑CT‑CCACATTGCG‑C‑GA‑GGCA‑TTTTTAAAAGTTGGCCAGCGCGCGAGGACC‑TGT‑GTTGAAACTACTCACCGTAACA‑ACGGCATGAATGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAGGCGACAAGCGAGCCGTAGCTGGTGGAAGCACATATCGGTg                                                        <  2:263597/225‑1 (MQ=255)
                                                                                                                                                                                      cttttctCGATTATCATGTCTCTTAGTCGTTGATCGCCGGG‑TCCCT‑CT‑ACACATTGCG‑C‑GA‑CGCA‑TTGTTAAAAGTTGGCCGGCGCGAGATGACC‑TGA‑TTTGAAGCTACCGACCTTAACA‑CCGGCAACAAGGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTGAGCCGGTTTGATTCATGGGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACCTATCGGTGCTCAc                                                   <  2:533883/233‑1 (MQ=255)
                                                                                                                                                                                      cGGTTATCGATCCCCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCAc                                                   <  1:1103431/239‑1 (MQ=255)
                                                                                                                                                                                      cGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCAc                                                   <  1:351249/239‑1 (MQ=255)
                                                                                                                                                                                       ggTTATCCATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGCCAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCAcg                                                  <  2:841411/239‑1 (MQ=255)
                                                                                                                                                                                       cgtttaccaatccccagcccccccGGGGTATATCGGGGGG‑GCGCTCCT‑CCACATGCCG‑G‑GA‑‑GGA‑ATTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCAcg                                                  <  1:1116507/217‑1 (MQ=255)
                                                                                                                                                                                           aTCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑GCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGcc                                              <  2:640143/239‑1 (MQ=255)
                                                                                                                                                                                              gATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCC‑CT‑CCCCCTTTCC‑C‑CC‑CGCC‑TTTTTTAAAAATGGGCTGCGCGCGGTGACC‑TGT‑TTTTTAACTACCGACCGCAAAA‑AAAGCATCCAAAAGGTGGAGGAAGGGGTGGAAGAAAAAAACCGCTTCTCCCCCTTGATTGATTAGGGGGAGGAAAAGAACAACAACCCGGATCTGGTGGGAGGACATATCAGGGCTCCCGCGCCgct                                           >  2:960102/1‑237 (MQ=255)
                                                                                                                                                                                                  aCCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTAcc                                       <  1:770124/239‑1 (MQ=255)
                                                                                                                                                                                                    cATGCCTCTCGGGCGTATAGCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGAA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGt                                     <  1:992369/239‑1 (MQ=255)
                                                                                                                                                                                                            tcCGGCGTAGATCGGCGGG‑TCGCT‑CT‑CCCCATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGTGCCGGATg                             <  2:406966/239‑1 (MQ=255)
                                                                                                                                                                                                            gccGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TGTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGTGCCGGATg                             <  2:1268008/238‑1 (MQ=255)
                                                                                                                                                                                                            cccGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGTGCCGGATg                             <  1:434749/238‑1 (MQ=255)
                                                                                                                                                                                                                gCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGc                                                                                       <  1:1263726/177‑1 (MQ=255)
                                                                                                                                                                                                                  gTATATCGGCGGT‑TCGCG‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTc                       <  1:953527/239‑1 (MQ=255)
                                                                                                                                                                                                                      atCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGga                   <  1:596948/239‑1 (MQ=255)
                                                                                                                                                                                                                       tCGGCGGG‑TCGCT‑CT‑CCTCATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TAT‑TTTGAAACTACTGACCTTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCTTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGgag                  <  1:570800/239‑1 (MQ=255)
                                                                                                                                                                                                                                      ct‑cCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGGAGATGCTGCGTGaa      >  1:340475/1‑238 (MQ=255)
                                                                                                                                                                                                                                      ct‑cCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGGAGATGCTGCGTGAAg     <  1:1251250/239‑1 (MQ=255)
                                                                                                                                                                                                                                      ct‑cCAAATTGCG‑C‑GT‑CGGT‑TTTTTAAAATGTGGCCTGCTCGCGATGACC‑TTG‑TTTGATACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTGCTCCCGTTTGATTGATGAGTCGAAGAAAGCGACAAGCGAGCCGTATCTTGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGGAGATGCTGCGTGAAg     <  2:499375/239‑1 (MQ=255)
                                                                                                                                                                                                                                          cacaTTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGGAGATGCTGCGTGAAGCCg  >  1:1008459/1‑239 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    
CTCGCACAACCATTTTTACTCGGGGCTTTCCCGCGGAATTATGGCAAACATCGCCCCTTGCCCGGATTTCATCTCAACGCTGAAAAATCTCTGGTGGCGGCTCGATCGCGCCCTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGG‑TCGCT‑CT‑CCACATTGCG‑C‑GA‑CGCA‑TTTTTAAAAGTTGGCCTGCGCGCGATGACC‑TGT‑TTTGAAACTACTGACCGTAACA‑ACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGGTGGAAGCACATATCGGTGCTCACGCGCCGTTTACCGTGCCGGATGCCGGTCTGGAGATGCTGCGTGAAGCCG  >  NC_000913/3019338‑3019807

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