Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,776,305 A→C 11.5% H130P (CAC→CCC)  yfjX → CP4‑57 prophage, putative antirestriction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,776,3050AC11.5% 194.7 / 2.1 80H130P (CAC→CCC) yfjXCP4‑57 prophage, putative antirestriction protein
Reads supporting (aligned to +/- strand):  ref base A (40/29);  new base C (2/9);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 2.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGG  >  NC_000913/2776067‑2776543
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              
tttgCCTGGGTGGACCGTTTTTGGGATGGGTACTGGGGGGGTGTCTGGTCCTTTTACCCGCTCAGCAATGGCGGAGGTTTTATGGCTCCGGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCa                                                                                                                                                                                                                                                <  2:539488/236‑1 (MQ=255)
ttCGCCTGGATGGACCGCTTCTGGGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCa                                                                                                                                                                                                                                                <  1:661691/239‑1 (MQ=255)
    ccTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCCCAAAGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGGAGCCGCAGGAAGCGCTGTCTGCCTGATTGAATATAGCCAGCCAGCCTGCCTCACCGAAAGCGATGCCATGACGGAGCACTAt                                                                                                                                                                                                                                            >  1:1080696/1‑239 (MQ=255)
             aCCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATAGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATTGTGCGGGAATGAGCGCGGAAGCCGCAGGGATCGCGTTCTGCCTGGTTGAAAATAGCCATCACGCCGGCCCCACCCAATGCGATTCCCTTACGGGAAACCATTAccccgt                                                                                                                                                                                                                                    >  2:419749/1‑234 (MQ=255)
              agttttgtataagggccggcgggggggTATTTTGGCCCTTTTCCCCGCCCAAGAAGGGCGGGGCGTGTGTGGCCCCTGGTGTTGAGGGTGACGATAAAGGGCGTTGGCTCAGCGGCATGAATGGCACTGGTGCGGAAAGGAGCGGGGAAGCCGCAGGAATCGTTGTCTGCCTGATGGAATATAGCCATCCCGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTg                                                                                                                                                                                                                                   <  2:572364/217‑1 (MQ=255)
                       ttgtgtgttcggggggggTTTTTGGTCCTTTTACCCGCTCAGAAATGGGGGGGGTTTTAGGGCTCCGGATGCTGGCGGGGACGAAAAATGGCTTTTGCTCAACGGCATGAATGGCAATGGTGCGAAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACg                                                                                                                                                                                                                         <  2:664269/226‑1 (MQ=255)
                          ggatgccggggggtgTCTCTGGTCCTTTTACCCGCCAAGGAAGGGCGGGGGGTTTAGGGTCCCTGATGCTGACGGTGGCGATAAGTGGCTTTGGCTCAACGGCATGGAGGGAAATGGGTCGGAAATGGGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCACCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGct                                                                                                                                                                                                                       <  2:519643/228‑1 (MQ=255)
                           gaGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTAAGc                                                                                                                                                                                                                        >  1:852022/1‑236 (MQ=255)
                            aGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCCTGAATGGGAATGGTGCGGAAATGGGCCCGGAAGCCGCAGGTATCGCTGTCTGCCTGCTTGAATATAGCCATACAGCCCGCCTCACCGAATTCGAGGCCATGACGGGGACCTATTACGCCCTCCGGGGTTACGccctc                                                                                                                                                                                                                    >  2:230645/1‑238 (MQ=255)
                                 tGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCCATGGTGCGGAAATGAGGGCGGAAGCCGCAGGGATCGCTGTCTGCCTGATTGAATAAAGCCATCACGCCTCCCCCACCGAAGGCGATGCCCTGCGCGAGCACACTTACccccccggggttaaccccctgcaca                                                                                                                                                                                                                >  2:891437/1‑217 (MQ=255)
                                   cGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGGGCGGGAATGAGCGCGGAAGCCGCAGGGATCGCTGTCTGCCTGATTGAAAATAGCCATACAGCCTGCCCCACCGCATGCGATGCCATTACGGGGCACCATTACCCCCCGCGGGAATACGCCCCGCAacac                                                                                                                                                                                                              >  2:901875/1‑238 (MQ=255)
                                       gtgtttttGGTTCTTTTACCCGCTCAGAAAGGGGGGGGGGTTTAGGGCTCCTGGTGCGGGGGGTGGAGAAAAAGGGCATTTGCCAAACGGAATGAATGGCAATGGGGGGGAAATGGGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCt                                                                                                                                                                                                           <  2:352663/231‑1 (MQ=255)
                                        gTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGGATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACCATTACCCCCTGCGGGATTACGCTCTGCACAACCCCGaa                                                                                                                                                                                                        >  2:980815/1‑239 (MQ=255)
                                            tCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGGACTATTTCCCCCTGCGGGATTACg                                                                                                                                                                                                                         >  1:572364/1‑219 (MQ=255)
                                            tCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGGTTACGCTCTGCAGCACCCTGAATc                                                                                                                                                                                                      >  1:253788/1‑238 (MQ=255)
                                            tCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGGAACGCTGTCTGCCTGATTGAATATAGCCCTCACGCCTGCCTCACCGAATGCCATGCCATGACGGGGCCCCTTTACCCCCTCCGGGATTCCGCTCTGCACCACCCCGAAAcc                                                                                                                                                                                                     >  2:27541/1‑239 (MQ=255)
                                              ttgTCCTTTATCCCGCTCAGCATGGGGGGGGGTTTTAGGGCTCATGGTGTGAGGGGGGCCGAAAAAGGGCTTTGCCTAAACGGAAGGAATGGAAAGGGGGGGGAAAGGACCGGGGAAGCCGGAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAg                                                                                                                                                                                                   <  2:991986/237‑1 (MQ=255)
                                                gTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAgcg                                                                                                                                                                                                 <  1:25143/239‑1 (MQ=255)
                                                 tCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGGGAGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCCACGGGACTACGCCTTGCAGCACCCTGGATCCAgcg                                                                                                                                                                                                 >  1:202743/1‑238 (MQ=255)
                                                     tttACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGGGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGCTTACGCTCTGCAGCACCCTGAATCCAACGCCCtt                                                                                                                                                                                            >  1:38061/1‑239 (MQ=255)
                                                      ttACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGGGCGGAAATGAGCGCGGAAGCCGCAGGAATCGCTGTCTGCCTGATTGAATATAGCCATCACCCCTGCCTCACCGAATGCGATGCCATGACGGAGAACTATTAcc                                                                                                                                                                                                                                        >  1:913817/1‑194 (MQ=255)
                                                              tCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGGTGCCATGACGGAGAACTATTTCCGCCTGCGGGCTTAACCCCTCCACCCCCCCGAAACCGACGCCCTTATggcccta                                                                                                                                                                                   >  2:857273/1‑233 (MQ=255)
                                                               cAGCAATGGCGGGGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGAATGAATGGCAATGGTGCGGAAAGGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCatta                                                                                                                                                                                  <  1:118004/239‑1 (MQ=255)
                                                               aaggaaTGGCGGGGCTTTTAGGGCTCCTGATGCTGACGGGGACGATAAAGGGCATTTGCTCAACGGCATGAATGGAAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCatta                                                                                                                                                                                  <  1:729691/235‑1 (MQ=255)
                                                                     tGGGGGGGCGTTTATGGCCCCTGATGCGGACGGTAACGATAAATGGCATTTGCTCAAGGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCCCCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGAtt                                                                                                                                                                            <  1:857273/239‑1 (MQ=255)
                                                                          gAGCGTTTATGGCTCCTGATGCTGACGGGGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGACTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGAACTATAACCCCCTGCGGGATTACGCCCTGCAGCCCCCCGAATCCAGCCCCATTATGGCCattatttgataaggg                                                                                                                                                                       >  2:758399/1‑230 (MQ=255)
                                                                               tttATGGCTCCGTGTGCTGCCGGTGACGATAAAGGGCATTTGCTCAACGGCATGAATGGAAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTTTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACt                                                                                                                                                                   <  1:672129/238‑1 (MQ=255)
                                                                               tttATGGCTCCGGATGCTGACGGGTACGATAAAGGGCTTTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCCCCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACt                                                                                                                                                                   <  1:891437/238‑1 (MQ=255)
                                                                                   tGGCTCTTGAGGGTTAGGGTGACGATAAATGGCTTTTGCTAAACGGCATGAATGGCAATGGTGCGGAAATGGGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTTATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCCCTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCatga                                                                                                                                                              <  1:336831/239‑1 (MQ=255)
                                                                                        gcTGATGCTGAGGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAAGGGGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCCCGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGaa                                                                                                                                                         <  1:351001/238‑1 (MQ=255)
                                                                                         cTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTCCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTCTGCGCATTATTTATTAAGGGAACCCAAGATGGaaa                                                                                                                                                        >  2:3627/1‑238 (MQ=255)
                                                                                         cTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTCCGGGATTACGCCTTGCAGCCCCCTGAATCCAGCCCCATTTTGCGCCTTTTTGATTAAGGaaaacactaatgaaaa                                                                                                                                                        >  2:68548/1‑224 (MQ=255)
                                                                                              gCTGACGGTGACGAAAAATGGCATTGGCCCAACGGCATGAAGGGCAATGGTGCGGAAATGGGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCCCGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACt                                                                                                                                                                   <  1:601224/223‑1 (MQ=255)
                                                                                               cTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCGGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAAGGAGCACTATTACCGCCTGCGGGATTACGCTCTGAAGCCCCCTGAATCCAGCGCCCTTATCCGCCTTATTGATTAAGGAAACTCATGATGGAACAGtaa                                                                                                                                                   >  2:24632/1‑236 (MQ=255)
                                                                                                  gagggTGACGATAAATGGCATTTGCTAAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTc                                                                                                                                                <  1:27541/236‑1 (MQ=255)
                                                                                                     ggTGACGATAATTGGTATGTGCCCAACGGGTGAAATGGCAAGGGTGGGAAAATGAGCCCGGACGCCGGCGGTAACGCTGTCGGCCTGATTGAATATAGCCAACCCGCCTGCCTCCCCGAATGCGATGCCATGGCGGAGCCCTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCAt                                                                                                                                              <  2:38061/238‑1 (MQ=255)
                                                                                                         aCGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCAATATGCGCATTTTTGATTAAGGAAACTCATGATGGAAccggtctcttataaac                                                                                                                                         >  2:601224/1‑224 (MQ=255)
                                                                                                             aaaaTGGCATTTGCCCAAGGGCATGAATGGCAATGGGGCGGAAATGAGCGGGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCCCCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGa                                                                                                                                      <  1:758399/237‑1 (MQ=255)
                                                                                                                aTGGCATTGGCCCAACGGCTGGACGGGCAAGGGTGCTGAAATGAGCGCGGAAGCCGAAGGTATGCCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCCCCGAATGCGATGCCAGGACGGAGCCCTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACAcc                                                                                                                                  <  1:1002310/239‑1 (MQ=255)
                                                                                                                     ttttGCTCACGGGTATGAATGCAATGGGGCCGAAAAGTGGCGGGGAAGCCCGGGTGATCGCTGTCTGCCGGTTTGATTATGGCCATCCCGCCAGCCTCACCGAATGCGGTGCCCGGGCGGAGCACTTTTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTAATGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGCCACCGGTa                                                                                                                              <  2:182166/237‑1 (MQ=255)
                                                                                                                     aTTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGAt                                                                                                                                                                             <  1:407891/191‑1 (MQ=255)
                                                                                                                     aTTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGt                                                                                                                               >  1:292103/1‑237 (MQ=255)
                                                                                                                     aTTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTAATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTTTTGAt                                                                                                                                                                             >  2:407891/1‑191 (MQ=255)
                                                                                                                      tttGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACGGAACCTCATCCCACAGACACCGGAAc                                                                                                                             >  1:1031109/1‑238 (MQ=255)
                                                                                                                              aCGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATgc                                                                                                                                                                                        >  2:415947/1‑171 (MQ=255)
                                                                                                                              aCGGCATGAATGGAAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATgc                                                                                                                                                                                        <  1:415947/171‑1 (MQ=255)
                                                                                                                              aCGGAATGAATGGCAAGGGGGCGGAAATGAGCGCGGAAGCCGAAGGTATCGCTGTCTGCCTGTTTGAATAGAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCTCTGa                                                                                                                    <  1:121851/239‑1 (MQ=255)
                                                                                                                               cGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCCCTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGa                                                                                                                                      <  2:1041596/220‑1 (MQ=255)
                                                                                                                               cGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGa                                                                                                                                      >  1:1041596/1‑220 (MQ=255)
                                                                                                                                      aagggcaggggTGCGAAAATGGCCCGGGAGGCCGCAGGTTCGCGTGTCGGCCTGATTGAATAAAGCCATCACCCCTGCCTCACCGAATGCGATGCCATGACGGAGCCCTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGcaca                                                                                                             <  1:368903/229‑1 (MQ=255)
                                                                                                                                             aTGGTGCGGAAAGGAGCGCGGAAGCCCCAGGTATCGGTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACgg                                                                                                     <  1:901875/239‑1 (MQ=255)
                                                                                                                                                gTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTAACTCATCCCACAGCAACCGTTACTTCCACTGACCCCACAAAGCACGGGaa                                                                                                  >  2:708043/1‑239 (MQ=255)
                                                                                                                                                    ggAAATGAGCGCGGAAGCCGCAGGTACCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGTCACCGGTACTTCCACTGACCGCACAACGCACGGTAAAAcg                                                                                              <  1:418554/239‑1 (MQ=255)
                                                                                                                                                     gAAATGGGCGCGGAAGCCGCAGGTATCGGTGTCTGCCTGATTGAATATAGCCATCCCGCCTGCCTCCCCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAAcg                                                                                              <  1:68548/238‑1 (MQ=255)
                                                                                                                                                     gAAATGGGCGCGGAAGCCGCAGGTACGCTTGTCGTCCTGAATGAATATAGCCAGCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCCCTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAAcg                                                                                              <  1:153985/238‑1 (MQ=255)
                                                                                                                                                      gaatccgcccgccaaGCCGCGGGTACCCCTGTCTGCCTGATTGAATAGAGCCATCCCGCCTGCCCTCCCGAATGCGATGCCAGGCCGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAAcg                                                                                              <  2:1080696/225‑1 (MQ=255)
                                                                                                                                                         tGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCCTGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGAACAACGCACGGTAAAACGCGCCt                                                                                         >  1:410728/1‑239 (MQ=255)
                                                                                                                                                                   aaGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCCCTCATCCCACAGACACCGGTa                                                                                                                              >  2:1028689/1‑193 (MQ=255)
                                                                                                                                                                   aaGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTa                                                                                                                              <  1:1028689/193‑1 (MQ=255)
                                                                                                                                                                    acccGCAGTTTGCGCTGTCTGCTTGATTAAATATCGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACgctgc                                                                                <  1:1058759/236‑1 (MQ=255)
                                                                                                                                                                       cccgcgggttttcctttttGCCTGTTTGAATATGCCCACCCCCCCTGCCCCCCTGAATGGCTTCACATGACGGAGCCCTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACgctgct                                                                               <  2:1031109/220‑1 (MQ=255)
                                                                                                                                                                                         tGCCTGATTGAATATAGCCATCACGCCTGCCTCACGGAATGCGATGCCATGACGGAGCACTATTACAGCATGAGGGATTACGCTCTGCAGAACCCTCAATCCAGCGCCATTATGAGCATTATTGATTAAGGATACTCATGATGGAACAGTCAATCACCACACAGACACAGGTACTTCAACCGACCGCACAAGGCACGGTAAAACGAGCCTTAACACTGCTTGACGTACACCTGCCCGaa                                                            >  1:78435/1‑239 (MQ=255)
                                                                                                                                                                                         tGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGaa                                                            >  1:964827/1‑239 (MQ=255)
                                                                                                                                                                                                     tatgGCCATCCCCCCTGCCCCACCGAATGGGTTGCCAGACGGGAGCCCTTTTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCa                                                <  2:292103/239‑1 (MQ=255)
                                                                                                                                                                                                       tggCCATCACGCCTGCCTCCCCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCAtt                                              <  1:980815/237‑1 (MQ=255)
                                                                                                                                                                                                       aaGCCACCCCGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCAtt                                              <  1:230645/238‑1 (MQ=255)
                                                                                                                                                                                                           cctcccgtctgcctCACCGAAGGCGGTGCCATGACGGAGCCCTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCAcc                                          <  1:21875/231‑1 (MQ=255)
                                                                                                                                                                                                             tCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCCCTCATCCCACAGACACCGGTACTTCCACTGACCGCCCAAGCCACGGTACAACCCGCCTTAACGCTCCTTACCCGACACCGTCGCGAAAAAGGCGGGGCATTCACCTc                                        >  2:443715/1‑239 (MQ=255)
                                                                                                                                                                                                             tCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCAttaccct                                         >  1:871763/1‑233 (MQ=255)
                                                                                                                                                                                                             tCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCt                                         >  1:711513/1‑238 (MQ=255)
                                                                                                                                                                                                             tCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACCCTCCTTGACCCACACCTCCCCAAAAAAGCGGGGGATTCACcccc                                       >  2:459908/1‑239 (MQ=255)
                                                                                                                                                                                                             tCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTAATTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTCCGCGAAACAGGCGTGGCATTCACCt                                         >  1:1043283/1‑238 (MQ=255)
                                                                                                                                                                                                             tCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAACGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACACACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCCGCGCGAAAAAGGAGGGGCATGCAccc                                         >  2:463724/1‑237 (MQ=255)
                                                                                                                                                                                                              cACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCCCGGGAAAACCAGCCTTACCCCGGGTTACCCGAAACCTGCGGCAAACAGACCGGGCTTTCCCCTc                                        >  2:606303/1‑238 (MQ=255)
                                                                                                                                                                                                                      cctCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCCCGAAACAGGCGTGGCATTCACCCCCACCCAg                                 >  2:428596/1‑237 (MQ=255)
                                                                                                                                                                                                                      cctCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACGGACCGCCCAACGCACGGTCAAACCCGCCTTAACCCTGCTTGACCGAAACCTCCGCGACAAAGACTGGCCATTAACCTCCACCCAggt                               >  1:1055030/1‑238 (MQ=255)
                                                                                                                                                                                                                                tGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAAAGCTGCTTGACCGACCCCTGCGCGAAACAGGCGGGGCATTCACCTCCACTCAGGCCGCCCGGGAc                     >  2:506109/1‑239 (MQ=255)
                                                                                                                                                                                                                                       cccttgcgggAGCACTTTTTCCGCCTGCGGGATTTCGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACCCCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGgctg               <  2:253788/231‑1 (MQ=255)
                                                                                                                                                                                                                                                  gCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAAcgcg                                                                                            <  1:784525/151‑1 (MQ=255)
                                                                                                                                                                                                                                                  gCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAAcgcg                                                                                            >  2:784525/1‑151 (MQ=255)
                                                                                                                                                                                                                                                  gCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGGGGCATTCCCCCCCCCCCAGCCTGCCCGGGACTGgctagaacggaaaaa    >  2:294101/1‑227 (MQ=255)
                                                                                                                                                                                                                                                   cACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACgctgct                                                                               >  1:215885/1‑163 (MQ=255)
                                                                                                                                                                                                                                                    aCTTTTACCGGCTGCGGGATTACCCTCTGCAGCACCCTGAATCCCGCCCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATgg  <  1:599883/239‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              
TTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGACGGAGCACTATTACCGCCTGCGGGATTACGCTCTGCAGCACCCTGAATCCAGCGCCATTATGCGCATTATTGATTAAGGATACTCATGATGGAACAGTCACTCATCCCACAGACACCGGTACTTCCACTGACCGCACAACGCACGGTAAAACGCGCCTTAACGCTGCTTGACCGACACCTGCGCGAAACAGGCGTGGCATTCACCTCCACTCAGGCTGCTCGTGACTGGCTGAAGCTGAAAATGG  >  NC_000913/2776067‑2776543

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