Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 249,538 T→G 26.4% pseudogene (533/1713 nt) lfhA ← pseudogene, flagellar system protein, promoterless fragment, flagellar biosynthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913249,5380TG26.4% 74.8 / 3.8 77pseudogene (533/1713 nt)lfhApseudogene, flagellar system protein, promoterless fragment, flagellar biosynthesis
Reads supporting (aligned to +/- strand):  ref base T (22/26);  new base G (19/7);  total (44/33)
Fisher's exact test for biased strand distribution p-value = 2.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGGTTACCCTGTGCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGGCGCTCAGGTT  >  NC_000913/249308‑249773
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                
gCGGGTTACCCTGTGCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGTTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCAt                                                                                                                                                                                                                                         >  1:435091/1‑238 (MQ=255)
 gctgtttcccttgtcTTTGCCCCCCAGCGCCACCTTTTTGTAACCGGGGCTTAAACTGTTGGGCGCGATCAGGGGAATGGTTTCCCAACTGCCCTGCTGCTCGCTGGTTTGGGTGGTGGTGCGGGTCAGCGTTTCGGGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCACCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCAtt                                                                                                                                                                                                                                        <  2:629212/224‑1 (MQ=255)
     ttACCCTGTGCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGGGGCATTCCCgg                                                                                                                                                                                                                                  >  1:474013/1‑239 (MQ=255)
             tgCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCCCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTGAAACCAAGCAGGGCGCTGAGCAGAAAAAACGGCCGGGGGGGGACTCCCgcggcaacc                                                                                                                                                                                                                           >  2:107618/1‑230 (MQ=255)
              gCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCCAAACAGGCAGGGGGGGCATTCCCGGTacccccc                                                                                                                                                                                                                         >  1:1190345/1‑237 (MQ=255)
              gCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGGaccaccc                                                                                                                                                                                                                         >  1:755586/1‑237 (MQ=255)
              gCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCCGTCCCGGTCACACCg                                                                                                                                                                                                                         >  1:152174/1‑238 (MQ=255)
              gCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCACTGGTTTCGGTGATGGTGCGGGTCAGGGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGGAAAAACAGCAGGCGCGGCATTCCCGGGACCCCCg                                                                                                                                                                                                                         >  1:1228219/1‑238 (MQ=255)
                  tGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTTCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAAAAGCAAAAACGGCAAGGGAGGCGTTTCCCGTACCACCGCCAg                                                                                                                                                                                                                     >  1:74636/1‑238 (MQ=255)
                  tGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGGCCTGACCAGCAAAAACGGCAGGGTGGGGCGTTCCCggaccccccccccc                                                                                                                                                                                                                      >  2:1150788/1‑227 (MQ=255)
                  tGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAGAAACGGCAGGTGCGGCATTCCCGGTACCCCCGcccg                                                                                                                                                                                                                     >  1:1011975/1‑236 (MQ=255)
                  tGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGCCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGc                                                                                                                                                                                                                    >  1:234142/1‑239 (MQ=255)
                  tGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGGTTGCTCATCCGCCAGCCGGTAAAACCCAGCAGGGCGCTGAACCGCAAAAAAGGACGGGGCGGCATTTCCGTGACCACCCCCac                                                                                                                                                                                                                     >  1:193031/1‑237 (MQ=255)
                    tCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCCGCCGGTAAAACCAAGCAGGGCGCTGAACAGCCAAAACGGCAGGGGGGGCCTTTCCGGGACCCCCGCCaccc                                                                                                                                                                                                                   >  1:614381/1‑235 (MQ=255)
                       accaGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGTACCACCCCCCGCACGa                                                                                                                                                                                                                >  1:1373424/1‑238 (MQ=255)
                        ccccccccccctgTTTGTAACCGAGGGTTTAAACGATTGGCTTCGTCGGGGGGATTGTTTCCCCGATGGCCCGGTGCTCGGTTGTTTCGGGGGATGGGGCGGTTCGGGGTTCGGGGGTTTTTCTCTTCCCCCCCGGCGCCCAGCGGCTGTTTGCCCACCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCGGGTGCGGCATTCCCGGTCCCACCGCCAGCACGaa                                                                                                                                                                                                               <  2:74874/227‑1 (MQ=255)
                            cgcCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCCGGCGCGGCATTCCCGGGACCACCGCCAGCACCAACATa                                                                                                                                                                                                           >  1:394942/1‑238 (MQ=255)
                              ccaccaGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCTCCGGACACCCGCCAGACCAAACACAAt                                                                                                                                                                                                         >  1:1314589/1‑238 (MQ=255)
                              ccaccaGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGGACCACACCCAGCACGAACATaaa                                                                                                                                                                                                         >  1:853514/1‑237 (MQ=255)
                              ccaccaGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCCGTAAAACCAAGCCGGGCCCTGAACAGCAAAAACGGCCGGTGCGGCAATCCCGGTACCACCACCACCACCAACCTAAt                                                                                                                                                                                                         >  1:227050/1‑238 (MQ=255)
                                     tttGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGTGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCCTCCGCCAGCCGGTACAACCCAGCACGGCGCTGAACAGCAAAAACGGCAGGGCGCGGAATCCCCGTACCACcccccgacccaaccaaaaccccgagg                                                                                                                                                                                                 >  2:215854/1‑217 (MQ=255)
                                     tttGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTgg                                                                                                                                                                                                 >  1:281267/1‑239 (MQ=255)
                                               aGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACTGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGGCCGGACATAATCCCGGGCGCGgtgtgagg                                                                                                                                                                                       >  1:598766/1‑235 (MQ=255)
                                               aGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGAGGCCTTTCCCGTACCAACGCCAGACAGAACCTAATACCGGTGgggggggtag                                                                                                                                                                                        >  1:91351/1‑232 (MQ=255)
                                               aGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGAAGGGCGCTGAAAAGCAAAAACGGCCGGGGCGGGATTTCCCGGACCCCACCCCGCACCAACACAAGACCgcgggcgggggaaag                                                                                                                                                                                       >  2:824480/1‑232 (MQ=255)
                                                     aaatggaggggctgaatcgggggAATGGTTTCCCAACTGACCTGCGTCTCGTGGTTTCCGGTGGTGGGGCGGGTCAGCGTTCCGAGGCTTTTCTCTTCCGCCGCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGa                                                                                                                                                                                 <  1:788703/218‑1 (MQ=255)
                                                           aTAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCCGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACAACCGCCAGGACGCACATAATCCCGGTGGCGGGTGAGAGCACCCACGGGCt                                                                                                                                                                           >  1:1060534/1‑239 (MQ=255)
                                                           aTAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTTAACAGCAAAAACGGCAAGGGCGGGCATTCCCGGACCACCCACCACACGAACATAAACCCGGGTGCgggggtgaagacccacggggg                                                                                                                                                                           >  2:1345501/1‑220 (MQ=255)
                                                           aTAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGAAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGGATTTCCGGTACGCACGCCAGCACGGAATAAAACACCGTGGCGGGTGAGAGCAccacaggg                                                                                                                                                                             >  1:200541/1‑231 (MQ=255)
                                                           aTAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGCTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCCGTACCACCGCCACCCCGAACATAATCCCCGTGGCGGTGTAGAGCACCCACGGGCt                                                                                                                                                                           >  1:950329/1‑239 (MQ=255)
                                                                cTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCCGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCGTTCCCGGTCCCACCGCCAGCAACAAAAAAATCCCgggggggggggacggccaccacgggccggcc                                                                                                                                                                       >  2:730804/1‑215 (MQ=255)
                                                                   gATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCACGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTGCCGGTACCCCCGCCACCACGAACACAATCCACGTGGCGGTGTAGAGCACCGACGGGCTTGCcaaca                                                                                                                                                                   >  1:39885/1‑239 (MQ=255)
                                                                   gATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTAACACCGCCAGCACGAACATAATCCCGGGGGCGGGGTAGAGCACCGACGGGCTTGCcagc                                                                                                                                                                    >  1:1131021/1‑238 (MQ=255)
                                                                    aTCAGGGGAATGGTTGCCCACCTGCCCGGCTGCTCGCTGGTTCGGGTGAGGGGGCGGTTCAGCTTTTCGGGGCTTTTCTCTCCCGCCCCCGCCGCAGGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCcagca                                                                                                                                                                   <  1:911307/238‑1 (MQ=255)
                                                                     tCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGGGGCGGTGTaga                                                                                                                                                                                        >  1:714209/1‑216 (MQ=255)
                                                                             agggtcaccaatcTGCCCGGCTGCTCGTTGGTGTCGGTGTTGGTGGGGGTCAGCTTTCCGAGTCTTTTCTCTCCCCCCTCCCCCCCCAGCGGCTTTTTGCTCACCCGCCAGCCGGTAAACCCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACg                                                                                                                                                         <  1:730804/227‑1 (MQ=255)
                                                                                   tCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTTCCGGTACCCCCGCCCGCACGAACATAATCCCGGTGGCGGGGTAAGGCACCGACGGGCTTGCCAGCAGCTGGGTACGCACGt                                                                                                                                                    >  1:1032700/1‑238 (MQ=255)
                                                                                   tCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGGGGCGGGGTAGAGCACCGACCGGCTTGCCAGCAGCCGGTTACCCCCGt                                                                                                                                                    >  1:1021769/1‑238 (MQ=255)
                                                                                   ccccACCTGCCCTGCTGCTCGCTGTTTGCGGTGAGGGTCGGGGTCAGCGTTTCGAGGTTTTTCTCTTCCGCCTCCGCCGCCGGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTc                                                                                                                                                   <  1:579205/238‑1 (MQ=255)
                                                                                         cTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGGACCACCGCCAGCCCGAAAATAAAACCCGGGGCGGGGGAGAGGACCGCCGGGGTTGCCCACAGGCGGGTAAGCACGGCaaaggg                                                                                                                                             >  2:1363074/1‑234 (MQ=255)
                                                                                          tGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCCCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCCTTCCCGCTCACACCGCCCGCACGGACACAAACCCCGGGGCGGGGGAAAGACCCCACCGGCCTGCCCGCAACCGGTGTCGCAcccccaggggg                                                                                                                                            >  2:1358756/1‑229 (MQ=255)
                                                                                                       gCTGGTATCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGTCGGCCTTCCCGGTACCCCCGCCAACCCGAACACAAACCCGGTGGGGGTGTGGAGCACCCACCGGCTTTCCAGCAGCCGGGTACGCACGCCCTGGGCGGTCACGCCCTTGt                                                                                                                              >  2:1318950/1‑239 (MQ=255)
                                                                                                                       gggggggggTCCGGGTTTCGAGGTTTTTCTCTCCCCCCTCCGCCGCCAGGGGCTTTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGt                                                                                                                 <  1:1075235/232‑1 (MQ=255)
                                                                                                                            cGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCCTTCCCGGTACCCCCGCCACCCCGAACATAATCCCGGGGGCGGGGTAGGGGCCCCACCGGCTTGCCCCCACCTGGGTACGGCCGCCAGTGGGGGATTTCCCCGTGTCGCCGAccgggggggggga                                                                                                           >  2:1124541/1‑227 (MQ=255)
                                                                                                                              ggTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGGGTAGGGGACCGACGGGCTTGCCAGCAGCTGGTTACCCACGTCATGGGCGGAATCCCCGTTGTGGCTGACGCGGGGGACGGTaa                                                                                                         >  1:1266286/1‑239 (MQ=255)
                                                                                                                                     gTTGCGGGGCTTTTCTCCGCCGCCTCCGCCGCCAGCGGCTTTTGCCCCATCCGCCACGCGGTAAACCCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTCCCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATcgc                                                                                                     <  1:107618/236‑1 (MQ=255)
                                                                                                                                              gttttCTCCTCCGCCCTCTCCGGCACGGGCTGTTTTTGTCTTCCCCCCGGCTGGAAAATCACGCCAGGGCGTGGAACAGCAAACAGGCCGGGTGGGCCTTCCCGGGTACCACCGCCAGCACGAACATAAGCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTg                                                                                           <  2:323308/238‑1 (MQ=255)
                                                                                                                                               ttttCTCTTCCCTCTCCGCCGCCAGCGGATGTGTGCGCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCGGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGga                                                                                         <  1:1044417/238‑1 (MQ=255)
                                                                                                                                               ttttCCCTCCCGCCCCCGGCGCCCGGGGTTGTTGACATACCCCCCAGCCGGAAAAACCAAGCAGGGCGCTGAACAGCAAAACCGGCAGGTGCGGCATTCCCGGTCCCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGgag                                                                                        <  1:1174470/239‑1 (MQ=255)
                                                                                                                                                 tcctctCCCCCCTCCGCCCCAGGCGGTTGTTTCCACCACCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGagag                                                                                      <  1:824480/237‑1 (MQ=255)
                                                                                                                                                    tccccccccccccccgccTGCGGTTTTTTCCCCCCCCCCCCGCCGGTAAAACCAGGCGGGGCGCTAAACGGAAAAACGGGGAGGTGGGGCATCCCCGGTCCCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGagcag                                                                                   <  2:474013/226‑1 (MQ=255)
                                                                                                                                                        ccgcctccgccGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTa                                                                                                                                                                                          <  1:1299122/133‑1 (MQ=255)
                                                                                                                                                        ccgcctccgccGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTa                                                                                                                                                                                          >  2:1299122/1‑133 (MQ=255)
                                                                                                                                                                            ttgttgcatcacccccaaccgggTAAACCAAAGCGGGGCGCTAAACAGAAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGccacc                                                             <  1:1193238/217‑1 (MQ=255)
                                                                                                                                                                              tttcgcccATCCGCCCGGCGGTAAAACCAGGCAGGGCGTTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTCCCCCCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAgg                                                          <  1:791305/232‑1 (MQ=255)
                                                                                                                                                                                tggCTCACGCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAAAGGAAAAAACGGCAGGTGGCGCATTCCCGGTGCCCGCGCCAGCACGAACATAATCCCGGGGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGcc                                                        <  1:15525/237‑1 (MQ=255)
                                                                                                                                                                                  tccccccgcccgggccgTAAAACCAGGGGGGCGCGACACAGCAAAAAACGGCGGGGGCGGGTTTCCCGTACCCACCGCCAGCACGAAAATAAGCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGt                                                      <  2:950329/223‑1 (MQ=255)
                                                                                                                                                                                  gCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCCGCGCCACCAGGccc                                                       >  1:330389/1‑237 (MQ=255)
                                                                                                                                                                                       cccGCCAGCCGGTAAAACCAAGCAGGGCGCAGAAAAGAAAAAACGGAGGGTGCGGCTTTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCg                                                 <  1:1150788/238‑1 (MQ=255)
                                                                                                                                                                                         cgcgcaccggggaaaCCCAGGCGGGCGCTTGAACGACAAAAGCCGCGGGTGCGGCTTCCCCGGCCCCCCCGCCGGCACGAACATAATCCCGGTGGGGGTGTAGGGCACGGACGGGCTTGCCAGCAGCGGGTTTCGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGa                                                <  2:545634/227‑1 (MQ=255)
                                                                                                                                                                                          gaccgccgtaaaaaCCAATGGGGTCGGTGAACAGAAACAGCGGCAGGTGCGGCCTTCCCGGCACCACCGCCAGCACCAACATAATCCCGGTGGCGGGGTAGAGCACCGACGGGGTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATg                                              <  2:91351/229‑1 (MQ=255)
                                                                                                                                                                                            ccgCCGGTAACCCCCAGCGGGGCCATGAAACACAAAACCGGGGGGGGGGGCTTTCCGGGCCCCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGt                                            <  2:1266286/237‑1 (MQ=255)
                                                                                                                                                                                              gccgtaaaaaCCAGCCGGGGCCCTGACAGAAAAAACCGGCAGGGGCGGCACTCCGGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGGCGGGCTTGCCAGCAGGTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGtca                                          <  2:1032700/233‑1 (MQ=255)
                                                                                                                                                                                               cccggtctaccctggcggggCCCTGACAGCAAAAATGGGCCGGTGCGGATGACCGGTGCCCCCCGCCAGCGCCAACATAACTCCGGTGGCGGTGTAGAACTCCGGCGGGCTTGCCAGCAGCTTGTTACGCACGTTGTGGGCGATTTCGCGGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGtca                                          <  2:634469/222‑1 (MQ=255)
                                                                                                                                                                                               ccGGTAAATCCGAGCGGGGCGCTGAAACGAAAACGCGGCGGGTGCGGCTGTCCCGGTACCACCGCCCGCCCGAACATAAGCCCGGGGGGCGTGGAGAGCACCGGCGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGtcat                                         <  2:302256/239‑1 (MQ=255)
                                                                                                                                                                                                cGGTAAACCCAAGCAGGGCGCTAACAAGAAAAACCGGGAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCCCCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGtcat                                         <  1:1085502/238‑1 (MQ=255)
                                                                                                                                                                                                cGGTAAAACCAGAGAGGCCGCGAAACAGAAAAAACGGCAGGGGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGtcat                                         <  1:1318950/238‑1 (MQ=255)
                                                                                                                                                                                                cGGAAAACCCAAGCGGGGGGAGGAACAAAAAAACGGGCGGGGGCGGCTTTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGtcat                                         <  2:10916/238‑1 (MQ=255)
                                                                                                                                                                                                  gTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCCCGGTGGAGAGCAGCAGGGAAGGGGTCTGCGCCCCCAGGCCGTCGCCGCTGGgcatca                                       >  1:10916/1‑238 (MQ=255)
                                                                                                                                                                                                      acccAAGCGGGGCGCTAAACAGCAAAAGCGGCAGGTGCGGCATTCCCGGTACCTCCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGcaca                                  <  1:182047/237‑1 (MQ=255)
                                                                                                                                                                                                         caggCAGGGCGCGGAACAGCAAAACGGGCGGGGGGGGCTTCCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATAc                               <  1:665595/236‑1 (MQ=255)
                                                                                                                                                                                                                   cgcTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATcgc                                                                                                     <  1:233452/160‑1 (MQ=255)
                                                                                                                                                                                                                                aaaaaCGGAGGGTGGGCGTTTCCGGTCACCCCCGCCAGCACGAACAGAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCgg            <  2:1079991/238‑1 (MQ=255)
                                                                                                                                                                                                                                    aCGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGa                                                                                                            >  1:1103258/1‑138 (MQ=255)
                                                                                                                                                                                                                                    aCGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGGCGCTc       >  1:260491/1‑239 (MQ=255)
                                                                                                                                                                                                                                          aGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCCCCTACTGCTGGGAGGCAGCAACGGCGCTCAgggt  >  1:1079991/1‑237 (MQ=255)
                                                                                                                                                                                                                                            ggTGCCGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGGTGCGGGGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGGCGCTCAGGtt  <  1:413812/238‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                
GCGGGTTACCCTGTGCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGGCGCTCAGGTT  >  NC_000913/249308‑249773

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: