Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,371,401 C→T A126V (GCT→GTT)  groL → chaperonin GroEL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,371,4010CT98.7% 253.6 / NA 76A126V (GCT→GTT) groLchaperonin GroEL
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (36/38);  minor base G (0/1);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GGTGCACCGACCATCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGC  >  NC_000913/4371157‑4371646
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
ggTGCACCGACCATCACCAAAGATGGTGTTTCCGTTGCTCGTGAACTCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGccgt                                                                                                                                                                                                                                                 >  2:30882/1‑248 (MQ=255)
 gTGCATCGACAATCCCCAAAGATGGTGTTTCCTTTGCTCGTGAAAACGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCAACACTGAAGGTCTTAAAGCTGTTGCTGCGGGCATGAACCCGATGGATCTGAAAAGTGGTATCTACAAAGCGGTTACCGTTGCAGtt                                                                                                                                                                                                                                                <  2:343462/251‑1 (MQ=255)
   gCACCGACCATCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTg                                                                                                                                                                                                                                               <  1:357399/250‑1 (MQ=255)
       cGACCATCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTgaagaa                                                                                                                                                                                                                                          <  1:390682/251‑1 (MQ=255)
        gACCATCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACCAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGTGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAAc                                                                                                                                                                                                                                         >  1:288845/1‑251 (MQ=255)
             tCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGCCCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGaa                                                                                                                                                                                                                                     >  2:414533/1‑250 (MQ=255)
               accaacgctggagtctCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGCGCAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAgc                                                                                                                                                                                                                                  >  2:413372/16‑251 (MQ=255)
                      aTGGTGTTTCCGTTGCTCGTGAAAGCGAACTGGAAGACACGTTCGAAAAGATGGGTGCGCTGATGGGGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGCCGGTACCACCACTGCAACTGTACGGGCTCAGGCTATCGTCACTGAAGGTCTGAAAGCTCGTGCCGCGTGCTTTAATCCGTGGTACCTGAAACGTCGTAGCGACAAAGCGGTTGCCTTTGCAGTTGAAGAACTGAAATCTCTGTcc                                                                                                                                                                                                                           <  2:430080/251‑1 (MQ=255)
                      aTGGTGATTCCGTTGCACGTGACATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCACAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGATGTcc                                                                                                                                                                                                                           >  2:88520/1‑251 (MQ=255)
                       tGGTGTTTCCGTTGCTCGAGAAATCGAACTGGAAGACAAGTTCGAAAATTTGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGCCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCg                                                                                                                                                                                                                          >  2:80298/1‑251 (MQ=255)
                        ggTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGATGCGCAGATGGTGCAAGCCGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGCAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGt                                                                                                                                                                                                                         >  1:292462/1‑251 (MQ=255)
                                cGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGCGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACTCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGct                                                                                                                                                                                                                 <  1:74624/251‑1 (MQ=255)
                                    gCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCTTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTg                                                                                                                                                                                                              <  1:393064/250‑1 (MQ=255)
                                                 aaCTGGAAGATACGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTGTCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGTACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTTCCATGCTCTGACTCTAAAGCGAtt                                                                                                                                                                                                <  2:254368/251‑1 (MQ=255)
                                                       aaGACAAGTTCGAAAATATGTGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCAGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCGGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTTCCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTTACTCTAAAGCGATTGCTCAg                                                                                                                                                                                          <  2:2291/251‑1 (MQ=255)
                                                         gACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACt                                                                                                                                                                                                                                        >  1:79633/1‑203 (MQ=255)
                                                         gACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGGTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACt                                                                                                                                                                                                                                        <  2:79633/203‑1 (MQ=255)
                                                          aCAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGTGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGtt                                                                                                                                                                                       <  2:288845/251‑1 (MQ=255)
                                                           cAAGTGCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTg                                                                                                                                                                                      <  1:114517/251‑1 (MQ=255)
                                                             aGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGt                                                                                                                                                                                    >  1:48385/1‑251 (MQ=255)
                                                                    aaaTATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCAtct                                                                                                                                                                             <  1:428537/251‑1 (MQ=255)
                                                                        atGGGAGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGACGAACTGAAAGCGCTTTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGc                                                                                                                                                                         >  2:376611/1‑251 (MQ=255)
                                                                              gcgcAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACt                                                                                                                                                                    <  2:422282/250‑1 (MQ=255)
                                                                                 aaGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCcga                                                                                                                                                                <  2:377375/250‑1 (MQ=255)
                                                                                  aGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCcgac                                                                                                                                                               <  1:293722/251‑1 (MQ=255)
                                                                                  aGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCcgac                                                                                                                                                               >  1:179688/1‑251 (MQ=255)
                                                                                        tGAAAGAAGTTGCCTCCAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGAGGTATCGACAAAGCGGTTACCGTTGCAGTTGACGAACTGAAAGCGCGGTCCGTACCATGCTCTGACTCTAAAGCGTTTGCTCAGGTTGGGACCATCTCCGCTACCTCCGCCGAAAcc                                                                                                                                                         >  2:18232/1‑251 (MQ=255)
                                                                                          aaaGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCg                                                                                                                                                        >  2:114426/1‑250 (MQ=255)
                                                                                                    ccTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTg                                                                                                                                             <  1:57450/251‑1 (MQ=255)
                                                                                                     ctctAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTg                                                                                                                                             >  2:139600/1‑250 (MQ=255)
                                                                                                     ctctAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTg                                                                                                                                             >  1:465396/1‑250 (MQ=255)
                                                                                                         aaaGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGc                                                                                                                                        >  2:88903/1‑251 (MQ=255)
                                                                                                            gCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTAATGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGGTGCAGTTGAAGAACTGAAAGCGCTGGCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGa                                                                                                                                     <  1:60749/251‑1 (MQ=255)
                                                                                                              aaccgacgCTGAAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGcaga                                                                                                                                     >  2:286608/4‑246 (MQ=255)
                                                                                                              aaaCGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGa                                                                                                                                     >  1:20986/1‑249 (MQ=255)
                                                                                                               aaCGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGc                                                                                                                                  >  1:343990/1‑251 (MQ=255)
                                                                                                                     gCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATgg                                                                                                                             >  2:110209/1‑250 (MQ=255)
                                                                                                                         cAGGCGACGGAACCACCACTGCAACGGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGAACTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGCTTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGACACCGTAGGTAAACTGATCGCTGAAGCGACGGACaa                                                                                                                         >  2:71365/1‑250 (MQ=255)
                                                                                                                                cGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTa                                                                                                                 <  2:465396/251‑1 (MQ=255)
                                                                                                                                   tACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTaacg                                                                                                              >  2:372496/1‑249 (MQ=255)
                                                                                                                                    accaccACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAaag                                                                                                              <  1:248190/250‑1 (MQ=255)
                                                                                                                                        ccacTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGc                                                                                                         <  1:153347/251‑1 (MQ=255)
                                                                                                                                             gCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTAt                                                                                                    >  1:266520/1‑251 (MQ=255)
                                                                                                                                                 ccGTCCTGGCTCAGGCTATCATCACTGAAAGTCGGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGTTTGGTACCATCTACGCTAACTCCTACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAATAAggttttttcacc                                                                                                <  2:308200/251‑11 (MQ=255)
                                                                                                                                                     aCTGGCTCAGGCTAGCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGGCGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGtt                                                                                             <  2:29708/250‑1 (MQ=255)
                                                                                                                                                       ttgCTCACGCTATCATCACTGAAAGTCTGCAAGCTGTAGCTGCAGGCCTGAACCCGATGGACCTGCAACGTAGTATCGACAAAGCGGTTACCGTTGCCGTTGAAGAACTGAAAGAGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTACACAGATCGCTGAAGCAATGGCCAAAGACGGTACCGAAGGCGTTATCACCGTTGa                                                                                           >  2:304152/3‑250 (MQ=255)
                                                                                                                                                        ggCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACAGTTGa                                                                                           >  1:222200/1‑249 (MQ=255)
                                                                                                                                                            cAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACgg                                                                                     <  1:110209/251‑1 (MQ=255)
                                                                                                                                                                 tATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCg                                                                                >  1:428016/1‑251 (MQ=255)
                                                                                                                                                                     atcaCTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCt                                                                            <  1:191623/251‑1 (MQ=255)
                                                                                                                                                                           gAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACTATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCTTTATCACCGTTGAAGACGGTACCGGTCTGCAg                                                                        <  1:369351/249‑1 (MQ=255)
                                                                                                                                                                             aGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACg                                                                    <  1:215216/251‑1 (MQ=255)
                                                                                                                                                                                   gAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAAc                                                                 <  1:396292/248‑1 (MQ=255)
                                                                                                                                                                                   gAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACt                                                                <  1:350518/249‑1 (MQ=255)
                                                                                                                                                                                     aaGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGAc                                                            <  1:29001/251‑1 (MQ=255)
                                                                                                                                                                                       gctgaggcggcGGGCAAGCACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTACAGAGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACGGTAGGTAAACTGATCGCTGAAGCGATGGACAACGTCGGAAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACAGGACGt                                                          >  2:278313/10‑251 (MQ=255)
                                                                                                                                                                                         tGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTcc                                                                                                                                                                  >  2:389370/1‑145 (MQ=255)
                                                                                                                                                                                         tGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTcc                                                                                                                                                                  <  1:389370/145‑1 (MQ=255)
                                                                                                                                                                                               ggcGGGCATGAACCCGATGGACCTGAAACGTGGTATCGCCAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCGCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTGGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAg                                                  <  2:188443/250‑1 (MQ=255)
                                                                                                                                                                                                 ccggCATGAACCCGATGGACCTGAAACGTGGTATCGCCAAAGCGGTTACCGTTGCAGTTGAAGCACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGCTCGCTGAAGGGATGGACACAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACAGGTCTGCAGGACGAACTGGACGGGGTTGAAGGt                                                >  2:206124/3‑251 (MQ=255)
                                                                                                                                                                                                 cGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGt                                                <  2:48385/251‑1 (MQ=255)
                                                                                                                                                                                                  gggCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGt                                                >  1:406673/1‑250 (MQ=255)
                                                                                                                                                                                                      aTGAACCCGATTGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACGCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAGCCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGGTATCACCGTTGAAGTCGGTACCGGTCTGCAGGACTAACTGGACGTGGTTGAAGGTATGCa                                           <  2:278633/251‑1 (MQ=255)
                                                                                                                                                                                                           cccGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGGGGTTGAAGGTATGCAGTTCg                                      <  2:179688/251‑1 (MQ=255)
                                                                                                                                                                                                            ccGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTgg                                                                                                                                                                                     >  1:37477/1‑107 (MQ=255)
                                                                                                                                                                                                            ccGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTgg                                                                                                                                                                                     <  2:37477/107‑1 (MQ=255)
                                                                                                                                                                                                               aTGGACCTGAAACGTGGTATCGACAAAGTGGTTACCGTTGCAGTTGGAGAACTGAAAGCGCTGTCCGTACCCTGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCg                                  >  1:185169/1‑251 (MQ=255)
                                                                                                                                                                                                                tGGACCTGCACCGTGGTATCGACAAAGCGGTTCCCGATGCAGTTGAAGACATGAAAGAGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGAAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGt                                 >  2:442344/1‑251 (MQ=255)
                                                                                                                                                                                                                 ggACCTGAAACGTGGTTTCGACAAAGCGGTTACCGTTGCAGTTGAAGAACGGAAAGCGCTGTCCGTACCAGGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTGCGACCGTg                                <  2:315319/251‑1 (MQ=255)
                                                                                                                                                                                                                 ggACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGt                                 <  1:114426/250‑1 (MQ=255)
                                                                                                                                                                                                                           aCGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGt                       >  1:35978/1‑251 (MQ=255)
                                                                                                                                                                                                                                  aTCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGGTATCACCGTTGAATACGGTACCGGTCTGCAGGACGAACTGGACGTGGGTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCtt                 <  2:80363/250‑1 (MQ=255)
                                                                                                                                                                                                                                  aTCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCGAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGTTAAACTGATCGCTGAAGCGAGGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGGATACGTTACCGGTCTGCATGACGAACTTGACGTGGGTGAAGGTATTCAGTGCGACCGGGGCTTCCTGTCTCCTTa                <  2:240509/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACtt             >  1:213082/1‑251 (MQ=255)
                                                                                                                                                                                                                                       cAAAGAGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGCTCGCTGACGCGCTGGACCAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTtc            >  2:235036/1‑250 (MQ=255)
                                                                                                                                                                                                                                                tACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGc  >  2:430431/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
GGTGCACCGACCATCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAACTGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGC  >  NC_000913/4371157‑4371646

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