Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,403,010 C→T R325C (CGT→TGT)  hflK → regulator of FtsH protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,403,0100CT97.2% 222.0 / NA 73R325C (CGT→TGT) hflKregulator of FtsH protease
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (30/39);  minor base G (1/1);  total (32/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGC  >  NC_000913/4402764‑4403255
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
ccgccgGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCt                                                                                                                                                                                                                                                   <  2:279292/251‑1 (MQ=255)
  gccgGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGt                                                                                                                                                                                                                                                 <  1:428498/251‑1 (MQ=255)
   ccgGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATGCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGta                                                                                                                                                                                                                                                <  2:124562/251‑1 (MQ=255)
       aagaagTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGCGTGTCTGTATATc                                                                                                                                                                                                                                            >  1:26650/1‑251 (MQ=255)
       aagaagTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAAGATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATc                                                                                                                                                                                                                                            <  1:170319/251‑1 (MQ=255)
        agaagaAAAAGCGGCGTTTGACGCTGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGCCGCGTATAACCACGAAGTTCAGCCGCGTGCGAACGGCCAGGCGAAACGTATCCTCGAAGAGGCGCGTGCGTACAACGCCCAGACAATCCTGGACGCTCAGGGTGATGGGGCGCGCTTTGCTAAACTTCTGCCCGAATATAAAGACGCGCAGGAAAATACTAGCGAGTGTCTCTATTTCg                                                                                                                                                                                                                                           >  2:49907/1‑251 (MQ=255)
          aagTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCgag                                                                                                                                                                                                                                         <  1:312286/251‑1 (MQ=255)
               aaaGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGACGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGAt                                                                                                                                                                                                                                    >  1:344765/1‑251 (MQ=255)
               aaaGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACgct                                                                                                                                                                                                                                    >  1:57923/1‑249 (MQ=255)
               aaaGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTGCTGCCTGAATATAAAGCCGCGCCGGAAATTTCTCGCGAGTGTCTGTATATCGAGACGAt                                                                                                                                                                                                                                    <  1:93049/251‑1 (MQ=255)
                aaGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATg                                                                                                                                                                                                                                   <  1:169363/251‑1 (MQ=255)
                   cggcgTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAAGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGaa                                                                                                                                                                                                                                <  2:319384/251‑1 (MQ=255)
                     gcgTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGaaaa                                                                                                                                                                                                                              <  1:51198/251‑1 (MQ=255)
                             cGATGCGATTGCCGCGCGAGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTgg                                                                                                                                                                                                                      >  2:502192/1‑251 (MQ=255)
                               aTGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTGTCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCGTTGCTAAACTTCTGCCGGAATATAATGCCGCGTCGGAAATTGCTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGt                                                                                                                                                                                                                    <  2:273654/251‑1 (MQ=255)
                                  cGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAac                                                                                                                                                                                                                 >  2:543535/1‑251 (MQ=255)
                                     aggccgctcgggaaagcgatctgCAATACATTCGTGAAGCAGAAGCGTATACCAACGCAGTGCAGCCGCGTGCGAACGGTCAGGCGCCACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGGGAAGTGGCGCGCTTTGCTAAACTGCGGCCGGAATAGAAAGCCGCGCCGGAAATTACTCGCGATGGTCTGTAGATCGAGACGATTGAAAAAGTGTTGGGTAACAcc                                                                                                                                                                                                              <  2:75798/229‑1 (MQ=255)
                                      tGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACAccc                                                                                                                                                                                                             >  1:371329/1‑251 (MQ=255)
                                               tGAAAACGAACAGCAATACATTAGTGACGCAGAAGCGTGTACAAACGAAGTTCAGCCGCGTGCGAACGGACAGGTGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTAAGGGTGAAGTGGAGCGATTTGCTAAAATACTGACGGAATATAAAGCCGCGGCGGAAAATCTTCGCGCGTGTCTGTATCTCGAGAAGATGGAAAACATGTTGGGTAACAACCGcctagggc                                                                                                                                                                                                    >  2:244127/1‑244 (MQ=255)
                                                gAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGGTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCt                                                                                                                                                                                                   <  1:249647/251‑1 (MQ=255)
                                                 aaaaCGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTCTATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTg                                                                                                                                                                                                  >  1:114341/1‑251 (MQ=255)
                                                      gAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTaa                                                                                                                                                                                             <  1:543535/251‑1 (MQ=255)
                                                       aaCAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAAc                                                                                                                                                                                            <  2:213797/251‑1 (MQ=255)
                                                          aGCAATACATTCGTGACGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGAt                                                                                                                                                                                         >  1:417081/1‑251 (MQ=255)
                                                             aaTACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCACGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATaa                                                                                                                                                                                       >  2:459373/1‑250 (MQ=255)
                                                                aCATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAggg                                                                                                                                                                                   >  1:265445/1‑250 (MQ=255)
                                                                 cATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTg                                                                                                                                                                                  <  1:426225/251‑1 (MQ=255)
                                                                     cGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCa                                                                                                                                                                               >  2:105228/1‑250 (MQ=255)
                                                                        gAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTGAACGATAAAGGTGGCAACCt                                                                                                                                                                           <  1:217794/251‑1 (MQ=255)
                                                                            cAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATg                                                                                                                                                                       <  1:459373/251‑1 (MQ=255)
                                                                            cAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATg                                                                                                                                                                       <  1:71536/251‑1 (MQ=255)
                                                                                   gTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGAAAAAAGTGTTGGGTAACACCCGCAAAGTACTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGc                                                                                                                                                                >  1:156513/1‑251 (MQ=255)
                                                                                     ataCCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGGATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCg                                                                                                                                                              >  2:316423/1‑251 (MQ=255)
                                                                                     ataCCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCTCTCCGGAAATTCCTCGCGCGTGTCTGTCTATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGGGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCg                                                                                                                                                              >  2:258539/1‑251 (MQ=255)
                                                                                     ataCCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGCATCCTCGAAGAGGCGCGTGCGTCCAAAGGCCAGACAATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGGGTCTGTAGATCGAGACGCTGGAAAAAGAGTTGGGGAACAACCGCCCAGTGATGGTTAACGAAAAAGGAGGCAACCTGCTGTGTCTGCCg                                                                                                                                                              >  2:249967/1‑251 (MQ=255)
                                                                                       aCCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGtt                                                                                                                                                            <  1:293437/251‑1 (MQ=255)
                                                                                        ccAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGCGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGtt                                                                                                                                                            >  1:283400/1‑250 (MQ=255)
                                                                                             gAAGGTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCa                                                                                                                                                      <  2:371329/251‑1 (MQ=255)
                                                                                                   cAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGc                                                                                                                                                 <  1:66678/250‑1 (MQ=255)
                                                                                                   cAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCt                                                                                                                                                <  2:10109/251‑1 (MQ=255)
                                                                                                              gAACGGTCGGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTGTGCTAAACTGCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATAGCGAGACGATGGAAAAAGTGTTGGGTAACATCCGCATAGTTCTGGTTAACGATAAAGGGGGCAACCTGATGGTTCTGCCGTTAGACCAGATTCTGAAAGGTGGTa                                                                                                                                     <  2:546650/251‑1 (MQ=255)
                                                                                                                   gTCAGGCGCAACGGATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGcc                                                                                                                                <  2:126868/251‑1 (MQ=255)
                                                                                                                     cAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGcccc                                                                                                                              <  1:265695/251‑1 (MQ=255)
                                                                                                                      aGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGCAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCt                                                                                                                             >  1:541327/1‑251 (MQ=255)
                                                                                                                        gcgcAACGTATCCGCGAAGAGGCGCTGGCGTACAAGTCCTATACCATCCTGGAAGCTCAGGGTGAAGTGGCGCTCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCTCTGGAAATTTCTCGCGACTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACTCCCTCAAAGTGCTTGTTAACGATAAAGGTGGCATCCTGATGGTTTTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTgc                                                                                                                           <  1:329443/251‑1 (MQ=255)
                                                                                                                         cgcAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTgcg                                                                                                                          >  2:446287/1‑251 (MQ=255)
                                                                                                                            aaCGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTGCTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTgcggcg                                                                                                                       <  1:178897/251‑1 (MQ=255)
                                                                                                                            aaCGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTgcggcg                                                                                                                       <  1:119342/251‑1 (MQ=255)
                                                                                                                             aCGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGCAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTgcggcg                                                                                                                       >  1:331763/1‑250 (MQ=255)
                                                                                                                             aCGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGa                                                                                                                      <  1:105228/251‑1 (MQ=255)
                                                                                                                                 aTCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAagag                                                                                                                  <  1:7295/251‑1 (MQ=255)
                                                                                                                                      cGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATa                                                                                                             <  2:313775/251‑1 (MQ=255)
                                                                                                                                              cgcgTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGcc                                                                                                     <  1:258476/251‑1 (MQ=255)
                                                                                                                                                     gtACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGGGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATc                                                                                              <  2:26650/251‑1 (MQ=255)
                                                                                                                                                       aCAAGGCCGAGACCATCGTGTAAGCTCAGGTTTAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGTAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGTTGGTTAAGGATAAAGGTGGCAGCCTGATGGTGCTGGCGTTAGACCAGATGCTGTATGATGGTATCGCTCCTGCGTTGAAGAGCGATAACGGTGCCAGCAATct                                                                                             <  2:31691/250‑1 (MQ=255)
                                                                                                                                                          aGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATctgctg                                                                                         <  1:258539/251‑1 (MQ=255)
                                                                                                                                                                    cATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGCAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGCGCGATAACGGTGCCAGCAATCTGCTGCGTCTgccgc                                                                               >  1:25827/1‑251 (MQ=255)
                                                                                                                                                                     ctCCTGGACGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAAAGGTGCCAGCAATCTGCTGCGTCTgccgcc                                                                              >  2:316928/2‑251 (MQ=255)
                                                                                                                                                                           gAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCctc                                                                        >  2:21620/1‑251 (MQ=255)
                                                                                                                                                                             agcgcAGGGTGAAGGGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTAGATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCtt                                                                      <  1:433071/247‑1 (MQ=255)
                                                                                                                                                                                    ggTGCAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAAATACTCGCGAGCGTCTGTATATCGAGACGATGGAAACAGTGACGGGAAACACCCGCCAAGCGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGCAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATCACGGTGCCAGCAATCTGCTGCGTCTGCCGCAAGCCTCTCCCTCcac                                                               >  2:299839/1‑251 (MQ=255)
                                                                                                                                                                                                      gCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGagcaagca                                              <  1:490403/250‑1 (MQ=255)
                                                                                                                                                                                                        tAAACTTCTGCCGGGATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGGGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGAGTCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTg                                                                                                       <  1:353841/191‑1 (MQ=255)
                                                                                                                                                                                                           aCTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACAcgt                                        <  1:502192/251‑1 (MQ=255)
                                                                                                                                                                                                                tGCCGGCATATAAAGCCGCGCCGGCAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCACAGTGCTGGTTAACGATACCGGTGGCAACCTGATGGTTCTGCCGTTAGCCCAGATGCTGAAAGGAGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCACGCAACACGTCGTcc                                   >  2:221377/1‑251 (MQ=255)
                                                                                                                                                                                                                  ccGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACCCCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCATCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTcca                                  <  1:333342/250‑1 (MQ=255)
                                                                                                                                                                                                                         ataAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTc                            >  2:153727/1‑249 (MQ=255)
                                                                                                                                                                                                                                    ccGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGAGTTGGGTAACACCCGCAAAGTGCTGGATAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCACCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGAtcttat               >  1:523997/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cGGAAATTACTCGCGGGAGTCTGTATAGCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCGCTTCCGCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATg              <  2:306458/251‑1 (MQ=255)
                                                                                                                                                                                                                                       gAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACAACCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATTGTTCTGCCGTTAGACCAGATGCTGACAGGTGGTAACGCCCCTGCGGCGCAGAGCGCTAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGAGATATTATgg             >  2:286682/1‑250 (MQ=255)
                                                                                                                                                                                                                                                cgcgAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGGTAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGCGCACGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAAcgccg   >  2:396338/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 gcgAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACCCCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCCGTCAGGGCGATATTATGGACCAAcgccgc  >  1:548183/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 gcgAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAAcgccgc  >  1:516182/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
CCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGC  >  NC_000913/4402764‑4403255

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