Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,609,768 T→G 24.0% G118G (GGT→GGG prfC → peptide chain release factor RF‑3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,609,7680TG24.0% 102.4 / 8.0 67G118G (GGT→GGGprfCpeptide chain release factor RF‑3
Reads supporting (aligned to +/- strand):  ref base T (18/27);  new base G (12/6);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 9.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAAAGGCCGTGGTTCCAACCAGCACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGTGATATCCGCGACCCGATGGAGCTGCTCGATGAAGTTGAGAACGAGCTGAAAATCGGCTGTGCGCCGATCACCTGGCCGATTGGCTGCGGCAAGCTGTTTAAAGGCGTTTACCACCTTTATAAAGACGAAACCTAT  >  NC_000913/4609551‑4609987
                                                                                                                                                                                                                          |                                                                                                                                                                                                                            
tataagggacaggTTTCAACCCGCCAGCTAAGGTGGACTGGGTGGAGATGGAAAAGCAGCGTGGGAGCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACCCCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCg                                                                                                                                                                                                                   <  2:596292/218‑1 (MQ=255)
tAAAAGGCCGTGGTTCCAACCAGCACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATTACTCCGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCg                                                                                                                                                                                                                   <  1:136413/229‑1 (MQ=255)
  aaaGGCCGTTTTTCCGACCAGCCCGCTAAAGTGGACTTGGTGGAGATGGAAAGGCAGCGTGGGATTTCCCTTACGACGTCTGTGGTGCAGTTTCCGTATCACGATTGCCTGGTTACCCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTa                                                                                                                                                                                                                <  2:283485/229‑1 (MQ=255)
       ccGTGGTTCCAACCAGCACGTTAAGTCGTACTGGATGGACATGGAAAAGCAGCGTGGGATCGCCATTACTACGTCTGTGATGCAGGTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTTCCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGt                                                                                                                                                                                                           <  1:591501/229‑1 (MQ=255)
          tGGTTCCAACCAGCACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATTGCTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGTCACCCCGGGGCACGAAGACTTCTCGGAAGATACCAATCGTACCCTGACGGCGGCGGGATGCGGCCTGAGGGGTCTCGCAGAACCAAAAGGAGTGGAAggcttgtctcttata                                                                                                                                                                                                         >  2:72481/1‑215 (MQ=255)
          tGGTTCCAACCAGCACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACCGCTGCCTGGTGGTTTTCGACCCCCCAAACGGGGTTggaacacttacccctaaa                                                                                                                                                                                                         >  2:719952/1‑212 (MQ=255)
                 aaccgtgaCGCTAAGTCGGACTGGATGGATATGGAAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATg                                                                                                                                                                                                  <  1:550607/221‑1 (MQ=255)
                       cACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGCCTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGCCGCGCCACAAGGAGTTTGAGCTCGTACCCGACACTTGATGGAGGtt                                                                                                                                                                                            >  1:735742/1‑228 (MQ=255)
                       cACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATGGACGCGCCACAAGGGGTTGAAGTCGTTACCCGTAAGCTGATGGAAGtt                                                                                                                                                                                            >  1:544450/1‑228 (MQ=255)
                       cACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGGGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGGGTTGAAGATCGTACCCGGAAGcccgcctcttata                                                                                                                                                                                           >  1:596293/1‑217 (MQ=255)
                       cACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGGGGACGTCTGCCCGAGGGTTATCGACGCCCACAAGAGGGGTGAAAACCGTACCCGTAAggttgtggagagaa                                                                                                                                                                                           >  2:490329/1‑216 (MQ=255)
                                                gggaAAAGCGGCGTGGGGTCTCCATTCTTCCTTGTGTGCTGCATTTTCCGATTCCCGATGGCCTGGTTAACCTGCCCGACCCCCCGGGGGACGAAGACTTCTGGGAAGATACCTACCGTCCCCGGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTCCCCGTCTGCGCGACACGCCGATCCt                                                                                                                                                                  <  2:285674/228‑1 (MQ=255)
                                                  gaaaaacgacttGGGTTCTCCTTTCCTCCGTTGTTGATGCATTTCTCGTATCCGGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTTCCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCa                                                                                                                                                                <  2:446809/218‑1 (MQ=255)
                                                  gaAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATTGAAGTTACCCGTCTGCGTGACACGCCGATCCTCa                                                                                                                                                                >  1:72006/1‑229 (MQ=255)
                                                  gaAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGGGTTGAAGATCGTACCCGTAAGCTGCTGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCa                                                                                                                                                                >  1:559326/1‑229 (MQ=255)
                                                  gaAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGCGGACTGCTGCCCGATGGTTATCGACGCCGCAAAAGGTGTTGAAACTCGTACCCGTAAGCTGATGGAAGTTACCCGTCCGCGCGACAGCCCGctcctc                                                                                                                                                                 >  1:219010/1‑228 (MQ=255)
                                                        catgtggggATCTCCATTACTAGGTCGGTGATGCGGTTTCCGTATCCCGATTGCCTGGTTAACCTGCTCAGCACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCACCttt                                                                                                                                                           <  2:153030/222‑1 (MQ=255)
                                                        cagcGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGGGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCACCTTTa                                                                                                                                                          >  1:17218/1‑229 (MQ=255)
                                                            gtgggattctcttttcttcctCTGTGAAGGTGTTTCAGTAACCCTGTTGCCCGGTTACCCTGCTCGACTCCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGGGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCCCCTTTATGaa                                                                                                                                                      <  2:391648/209‑1 (MQ=255)
                                                             tGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGCCTGCGCGACACCCCGATCCTCACCTTTATGAAc                                                                                                                                                     >  1:10550/1‑229 (MQ=255)
                                                                     ccATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGCCTTCTCGGCAGATACCTATCGTACCCTGACGGCGGCGGACTGCTGCCTGATGGTTATCGCAGCCGCCAAAGGGGGTGAAGCTCGTACCCGTAAGCTGATGGAAGTTACCCTGCTGCACGCCACGCCCGTCCTCACCCTTATGAAAAAAATTg                                                                                                                                              >  2:728503/1‑228 (MQ=255)
                                                                      ctgtactgcGTCTGGGATGCCGTTTCCGTATCACGATTGCCTGGTTTACCTGCTCGACGCCCCGGGGCACGAAGTCTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCGGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTg                                                                                                                                              <  1:222458/224‑1 (MQ=255)
                                                                      cATTACTCCGTCTGGGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTCACCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACt                                                                                                                                                <  1:299396/225‑1 (MQ=255)
                                                                       cttctcacctgtgtgAGGCAGTTTCCGTTACTCAATTGCCCGGTTAACCGGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTGTCGTACCCGGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGAc                                                                                                                                            <  2:207017/217‑1 (MQ=255)
                                                                              cGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAACAGGGGTTGAAGATCGTAACCGTAAGCTGATGGAAGTTAACCGTCTGCTCGAGACGCCGATCCTAACCTTTATGAACAAACTTGCCCGTGgta                                                                                                                                     >  1:731343/1‑228 (MQ=255)
                                                                              cGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGGGatat                                                                                                                                    >  1:119019/1‑229 (MQ=255)
                                                                                tCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGGAAGCTGATGGCAGTTACCCGCCTGCGCGACACGCCGGTCCTCACCTTTATTAACAACTTTAACCGTGatattc                                                                                                                                  >  1:600467/1‑227 (MQ=255)
                                                                                tCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGAGTTGAAAATCGTACCCGTAAGCTGGTGGAGGTTACCCGTCTGCGGGACAGGCCGATCCCAACCTTTATTAACAAACTTGACCGGGTTATcc                                                                                                                                  >  1:91849/1‑229 (MQ=255)
                                                                                tCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCAGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATTGACCACCCAAAAGGGGTTGAGGATCGTACCCGTGAGCTGATGGAAGTTAACCGCGTGCGGGCAATGCCGATCCTCACCCTTATGCAAAAAATTGACCggagtatcc                                                                                                                                  >  2:37261/1‑221 (MQ=255)
                                                                                      aggaatgttCCGTATCGCGATGGGCTGGTTAACCTGCCGGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGGGGTGGACTGCTGCCTGAGGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTCCCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGTGATATCCGCGAcc                                                                                                                            <  2:219010/222‑1 (MQ=255)
                                                                                                cGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGGGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGTGATATCCGCGACCCGATGGAgct                                                                                                                  >  1:377606/1‑229 (MQ=255)
                                                                                                 gTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTg                                                                                                                                                                                                                          <  1:158055/125‑1 (MQ=255)
                                                                                                    tCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGGGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGGGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGGCAAGCCGATCCTCACCCTTATGACCAAACTTTACCGCGTAATCCCCGAACCCATgaaccgtacg                                                                                                              >  2:709121/1‑220 (MQ=255)
                                                                                                                  tgaaccgtcTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATGCCTATCGTACCCTTACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGTGATATCCGCGACCCGATGGAGCTGCTCGATGAAGTTGAGaa                                                                                                <  1:734428/221‑1 (MQ=255)
                                                                                                                           tCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGTGCTATCCGCGACCCGATGGAGCTGCTCGATGAAGTTGAGAACGAGCTgca                                                                                       >  1:515499/1‑227 (MQ=255)
                                                                                                                           tCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGTGATATCCGCTACCCGATGGAGCTGATCGATGAAGTTGAGAACCAGCTgca                                                                                       >  1:614689/1‑227 (MQ=255)
                                                                                                                           tCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGTGATATCCGCGACCCGATGGAGCTGCTCGATGAAGTTGAGAACGAGCTGaa                                                                                       >  1:534711/1‑229 (MQ=255)
                                                                                                                              accccccGGGGCGCGAAGGCTTCTGGGAAGATACCCAAGGTATCCGGGCGGCGGTGGCCTGTGGTCTGATGGTTACCGAGGCCGCAAAAGTGGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCCCCTTTATGAACAAACTTGACCGTGATAGCCGCTACCCGATGGAGCTGCTCGATGAAGTTGAGCACGAGCTGAAAAt                                                                                    <  2:665280/226‑1 (MQ=255)
                                                                                                                                 ccgcGGGGCACGAAGACTTCTCGGAAGATACCGTTCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTGATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGTGATATCCGCGACCCGATGGAGCTGCTCGATGAAGTTGAGAACGAGCTGAAAATCgg                                                                                 <  1:360074/226‑1 (MQ=255)
                                                                                                                                 ccccGGGGCCAGGAGACTTCTCGGAAGATACCTAGCGTACCCTGACGGCGGTGGACTGATGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGTGATATCCGCGACCCGATGGAGCTGCTCGATGAAGTTGAGAACGAGCTGAAAATCgg                                                                                 <  1:623549/229‑1 (MQ=255)
                                                                                                                                 ccccGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGGGTTGAAGATCGTACCCGTAAGCTGATGGAAGGTACCCGTCTGCGCGACACGCCGATCCTCACCTTTTTGAACAAACTTTACCGCGATATACCCGCCCCGATGGAGCTGCCCGGTGAGGTTggaaaaggacgtaaaaaca                                                                                  >  2:377815/1‑210 (MQ=255)
                                                                                                                                   ctggggCACGTGAACTTCCCGGAAGTTCCCTATCGTGCCCGGACGGCGGGGGCCTGCCGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTCCCCGTCTGCGCGACCCGCCGATCCTCGCCTTTATGAACAAACTTGACCGTGATATCCGCGACCCGATGGAGCTGCTCGATGAAGTTGAGAACGAGCTGAAAATCGGCt                                                                               <  2:158391/227‑1 (MQ=255)
                                                                                                                                    cGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGCCTTCGCGACACGCCGATCCTCACCTTTTTGAAGAAAATTAGCCCGGCTATCCCCCACCCGATGAGGCGGCTCGTGGACGGTCAGAACGAAATGAAAATCGTCtg                                                                              >  2:307303/1‑229 (MQ=255)
                                                                                                                                       ggCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGGGATATCCGCGACCCGATGGAGCTGCTCGATGAAGTTGGAACCGAGCTGAAAATCGGCtgt                                                                             >  1:468361/1‑227 (MQ=255)
                                                                                                                                            gAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGAAGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGACGTTACCCGTCTGCGCGACACGCCGACCCTCACCTTTATGAACAAACTTGACCGGGATATCCCCGACCCGCTGGAGCTGCCCGAGTAATTTGGAAACGAGCTCAAAACCGGcgtggcgcgga                                                                      >  2:623548/1‑219 (MQ=255)
                                                                                                                                                   tctcGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGAGTTGAAGATCGTACCCGTAAGCTGATGGAAGGTACCCGTCTGCGCGACACGCCGTTCCTCACCTTTATGAACAAACTTGCCCGGGCTAGCCGCGACCCGATGGAGCTGCTCGATGAAGTTGAGACCGACCTGGAAATCGGCACGGCGCCGATAACCTg                                                               >  2:198121/1‑229 (MQ=255)
                                                                                                                                                           gATACCTATCGTCCCCTGACGGCGGGGGACTGCCGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGTGATATCCGCGACCCGATGGAGCTGCTCGATGAAGTTGAGAACGAGCTGAAAATCGGCTGTGCGCCGATCACCTGGCCGAtt                                                        <  1:167866/228‑1 (MQ=255)
                                                                                                                                                                   tCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGTGATATCCGCGACCCGATGGAGCTGCTCGATGAAGTTGAGAACGAGCTGAAAATCGGCTGTGCGCCGATCCCCTGGCCGATTGCGTGCGGc                                               >  1:560602/1‑229 (MQ=255)
                                                                                                                                                                        cccGGACGGGGGTGAGCTGCGGCTTGATGGTTAGCGACCCCGCAAAAGTGGTTGAAGATGGTACCGTTAAGCGTGTGGAAGTTCCCCGTCTGCGCGACCCGCCGATCCTCACCGTTATGAACAAACTTGACCGTGATATCCGCGACCCGATGGAGCTGCTCGATGAAGTTGAGAACGAGCTGAAAATCGGCTGTGCGCCGATCACCTGGCCGATTGGCTGCGGCAAg                                            <  2:618836/227‑1 (MQ=255)
                                                                                                                                                                           tGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGTGATATCCCCGACCCGATGGAGCTGCTCGATGAAGTTGAGAACGAGCTGAAAATCGGCTGTGCGCCGATCACCTGGCCGATTGGCTGCGGCACGCTGtt                                       >  1:391763/1‑229 (MQ=255)
                                                                                                                                                                             gcggcggTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGTGATATCCGCGACCCGATGGAGCTGCTCGATGAAGTTGAGAACGAGCTGAAAATCGGCTGTGCGCCGATCACCTGGCCGATTGGCTGCGGCAAGCTGTTTa                                     <  1:490329/228‑1 (MQ=255)
                                                                                                                                                                             acgggggggtgatgcgtccTGATGGTTATCCACGTCGCAAAAGGTTGTTAAGATCGTCCCCGTAAGCTGATGGAAGTTTCCCGGCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGTGATATCCGCGACCCGATGGAGCTGCTCGATGAAGTTGAGAACGAGCTGAAAATCGGCTGTGCGCCGATCACCTGGCCGATTGGCTGCGGCAAGCTGTTTa                                     <  2:468361/212‑1 (MQ=255)
                                                                                                                                                                             aCGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGTGATATCCGCGACCCGATGGAGCTGCTCGATGAAGTTGAGAACGAGCTGAAAATCGGCTGTGCGCCGATCACCTGGCCGATTGGCTGCGGCAAg                                            >  1:618837/1‑222 (MQ=255)
                                                                                                                                                                             aCGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCCCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACGAACTTGACCGTGATATCCGCGACCCGCTGGAGCTGCTCGATGAAGTTGAGAACGGGCTGAAAAACGGCTGTGGACCGGTCCCCTGGCCGGTTGGCTCGGGCAACCTGGTTa                                     >  2:172566/1‑229 (MQ=255)
                                                                                                                                                                                gcggTGGACTGCTGCCTGATGGTTATCGCCGCCGCACAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGTGACATCCGCGACCCGATCGAGCTCCTCGATGAAGTTGAGAACCGGCGCAAAACCGGCTGTGCGCCGACCCCCCTGCCGATTGGCCGGGGCAACCTTTTTaaa                                   >  2:273199/1‑228 (MQ=255)
                                                                                                                                                                                        cgtctgcCTAATGGTTAGCGGCGCCGCAAAAGGGGTTGAAGATCCTCCCCGTAAGCTGATGGAAGTTACCCGGCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGTGATATCCGCGACCCGATGGAGCTGCTCGATGAAGTTGAGAACGAGCTGAAAATCGGCTGTGCGCCGATCACCTGGCCGATTGGCTGCGGCAAGCTGTTTAAAGGCGTTTac                          <  1:37261/226‑1 (MQ=255)
                                                                                                                                                                                         tgctgcCTGATGGTTATCGGCCCCGCAAAAGTTTTTGAACATCGTACCCTTAAGCTGATGGAAGTTCCCCGTCTGCGCGACCCGCCGATCCTCACCTTTATGAACAAACTTGACCGTGATATCCGCGACCCGATGGAGCTGCTCGATGAAGTTGAGAACGAGCTGAAAATCGGCTGTGCGCCGATCACCTGGCCGATTGGCTGCGGCAAGCTTTTTAAAGGCGTTTac                          <  1:342411/228‑1 (MQ=255)
                                                                                                                                                                                                   tGGTTTACGACGCCACAAAAGTGTTTGAACACGGGTCCCGTAAGCTGGTGGAAGTTCCCCGTCGGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGTGATATCCGCGACCCGATGGAGCTGCTCGATGAAGTTGAGAACGAGCTGAAAATCGGCTGTGCGCCGATCACCTGGCCGATTGGCTGCGGCAAGCTGTTTAAAGGCGTTTACCACCTTtata                <  2:391763/228‑1 (MQ=255)
                                                                                                                                                                                                       taacgccgccgcAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGGCGATCCTCCCCTTTATGAACAAACTTGACCGTGATATCCGCGACCCGATGGAGCTGCTCGATGAAGTTGAGAACGAGCTGAAAATCGGCTGTGCGCCGATCACCTGGCCGATTGGCTGCGGCAAGCTGTTTAAAGGCGTTTACCACCTTTATAAAGa            <  1:203084/222‑1 (MQ=255)
                                                                                                                                                                                                       tATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGTGATATCCGCGACCCGATGGAGCTGCTCGATGAAGTTGAGAACGAGCTGAAAATCGGCTGTGCGCCGATCACCTGGCCGATTGGCTGCGGCAAGCTGTTTAAAGCCGTTTACCACCTTTATAAAGAc           >  1:455166/1‑229 (MQ=255)
                                                                                                                                                                                                            acgtcccAAAAGGTTTTGAAGCTCGTACCCGTAAGCGGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGTGATATCCGCGACCCGATGGAGCTGCTCGATGAAGTTGAGAACGAGCTGAAAATCGGCTGTGCGCCGATCACCTGGCCGATTGGCTGCGGCAAGCTGTTTAAAGGCGTTTACCACCTTTATAAAGACGAAAc      <  1:709122/223‑1 (MQ=255)
                                                                                                                                                                                                              gccgcAAAAGGGGTTGAAGTTCGTACCCGTAGGCTGATGGAAGTTCCCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGTGATACCCGCGACCCGATGGAGCTGCTCGATGAAGTTGAGAACGAGCTGAAAATCGGCTGTGCGCCGATCACCTGGCCGATTGGCTGCGGCAAGCTGTTTAAAGGCGTTTACCACCTTTATAAAGACGAAACCt    <  2:515499/229‑1 (MQ=255)
                                                                                                                                                                                                              gccgcAAAAGGGGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTCCCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGTGATATCCGCGACCCGATGGAGCTGCTCGATGAAGTTGAGAACGAGCTGAAAATCGGCTGTGCGCCGATCACCTGGCCGATTGGCTGCGGCAAGCTGTTTAAAGGCGTTTACCACCTTTATAAAGACGAAACCt    <  1:719953/229‑1 (MQ=255)
                                                                                                                                                                                                              cccgcAAATGGTGTGGAAGCTCGTCCCCGTAATCTGATGGAAGTTACCCGTCTGCGCGACAGGCCGATCCTCCCCTTTATGAACAAACTTGACCGTTATAGCCGCGACCCGAGGGAGCTGCTCGATGAAGTTGAGAACGAGCTGAAACTCGGCTGTGCGCCGATCACCTGGCCGATTGGCTGCGGCAAGCTGTTTAAAGGCGTTTACCACCTTTATAAAGACGAAACCt    <  2:72006/228‑1 (MQ=255)
                                                                                                                                                                                                                cgcAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCACCTTTATGACAAAACTTGACCGTGATATCCGCGa                                                                                                                              >  2:76761/1‑105 (MQ=255)
                                                                                                                                                                                                                cgcAAAAGGTGTGGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGTGATATCCGCGa                                                                                                                              <  1:76761/105‑1 (MQ=255)
                                                                                                                                                                                                                cccAAAAGGTGTTGAAGACCGTACCCGTAAGTTTAGGGAATTTCCCCGTCTGGGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGTGATATCCGCGACCCGATGGAGCTGCTCGATGAAGTTGAGAACGAGCTGAAAATCGGCTGTGCGCCGATCACCTGGCCGATTGGCTGCGGCAAGCTGTTTAAAGGCGTTTACCACCTTTATAAAGACGAAACCTAt  <  2:119019/227‑1 (MQ=255)
                                                                                                                                                                                                                          |                                                                                                                                                                                                                            
TAAAAGGCCGTGGTTCCAACCAGCACGCTAAGTCGGACTGGATGGAGATGGAAAAGCAGCGTGGGATCTCCATTACTACGTCTGTGATGCAGTTTCCGTATCACGATTGCCTGGTTAACCTGCTCGACACCCCGGGGCACGAAGACTTCTCGGAAGATACCTATCGTACCCTGACGGCGGTGGACTGCTGCCTGATGGTTATCGACGCCGCAAAAGGTGTTGAAGATCGTACCCGTAAGCTGATGGAAGTTACCCGTCTGCGCGACACGCCGATCCTCACCTTTATGAACAAACTTGACCGTGATATCCGCGACCCGATGGAGCTGCTCGATGAAGTTGAGAACGAGCTGAAAATCGGCTGTGCGCCGATCACCTGGCCGATTGGCTGCGGCAAGCTGTTTAAAGGCGTTTACCACCTTTATAAAGACGAAACCTAT  >  NC_000913/4609551‑4609987

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