Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 890,913 Δ18,382 bp [ybjC][ltaE] 21 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 890913 909294 18382 83 [0] [1] 84 [ybjC]–[ltaE] [ybjC],nfsA,rimK,ybjN,potF,potG,potH,potI,ybjO,rlmC,artJ,artM,artQ,artI,artP,ybjP,ybjQ,amiD,ybjS,ybjT,[ltaE]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 8909120 (0.000)83 (1.060) 73/492 0.0 100% intergenic (‑159/‑1) grxA/ybjC reduced glutaredoxin 1/DUF1418 domain‑containing protein YbjC
?NC_000913 909295 = 0 (0.000)intergenic (‑1/+36) ltaE/poxB low‑specificity L‑threonine aldolase/pyruvate oxidase

CATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTATGCGCGCGATCGGTAAATTGCCTAAAGGCGTGTTGATACTGGAATTTATCGGAATGATGCTACTGGCGGTGGCGCTGCTGTCGGTAAGCGACTCCCTGTCGCTGCCTGAGCCATTTTCTCGGCCAGAAGTGCAGATTCTGATGATTTTTCTCGGTGTTTTGCTCATGCTTCCCGCTGCGGTGGTGGTTATTCTTCAGGTGGCAAAACGTCTTGCCCCACAGCTGATGAACCGTCCACCGCAATATTCA  >  NC_000913/890666‑891161
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                         
catcGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggc                                                                                                                                                                                                                                                        >  1:204720‑M1/1‑247 (MQ=255)
 ctcgCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCTAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTCACCTGTTGCATTACTTGCTACAAGCTATAACTGTTAAACACAATACCGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcat                                                                                                                                                                                                                                                      >  2:634091‑M1/2‑246 (MQ=255)
 atcGCGTGCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCTCTTTGCAGAGGGAAGCGTggcat                                                                                                                                                                                                                                                      <  2:590701‑M1/251‑6 (MQ=255)
   cgcgTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgt                                                                                                                                                                                                                                                    <  2:135017‑M1/251‑8 (MQ=255)
      gTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtcct                                                                                                                                                                                                                                                 <  1:637588‑M1/251‑11 (MQ=255)
         cATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttat                                                                                                                                                                                                                                              >  2:559204‑M1/1‑238 (MQ=255)
                 aaTTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgccgg                                                                                                                                                                                                                                      >  2:287393‑M1/1‑230 (MQ=255)
                  aTTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacggg                                                                                                                                                                                                                                     >  1:567678‑M1/1‑229 (MQ=255)
                  aTTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacggg                                                                                                                                                                                                                                     <  1:169154‑M1/251‑23 (MQ=255)
                         aGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgcc                                                                                                                                                                                                                              <  1:583862‑M1/251‑30 (MQ=255)
                          gCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgcc                                                                                                                                                                                                                              >  1:161334‑M1/1‑221 (MQ=255)
                            cAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccacc                                                                                                                                                                                                                           >  2:172201‑M1/1‑219 (MQ=255)
                                   tttGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccct                                                                                                                                                                                                                         <  2:377283‑M1/246‑35 (MQ=255)
                                    ttGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccacccttttta                                                                                                                                                                                                                    <  1:520865‑M1/250‑40 (MQ=255)
                                     tGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttac                                                                                                                                                                                                                   <  1:604455‑M1/250‑41 (MQ=255)
                                        acgcacACAGTAAGGGCAACCCGAACGACCAAAAATAACGGATTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttacctt                                                                                                                                                                                                                >  1:145760‑M1/1‑207 (MQ=255)
                                                       gCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcg                                                                                                                                                                                                <  2:204720‑M1/251‑60 (MQ=255)
                                                             cGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagtt                                                                                                                                                                                          <  1:104881‑M1/251‑66 (MQ=255)
                                                             cGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagtt                                                                                                                                                                                          <  1:248194‑M1/251‑66 (MQ=255)
                                                              gAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttc                                                                                                                                                                                         >  2:53307‑M1/1‑185 (MQ=255)
                                                                      aaaaaTAACGGGTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcactt                                                                                                                                                                                 <  2:570105‑M1/251‑75 (MQ=255)
                                                                       aaaaTAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttcttatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttc                                                                                                                                                                                >  1:215785‑M1/1‑176 (MQ=255)
                                                                            aaCGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaactgccaccctttttaccttagccagtttgttttcgccaggtcgatcacttcatcac                                                                                                                                                                           >  2:217052‑M1/1‑171 (MQ=255)
                                                                                 tttGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgt                                                                                                                                                                      >  1:570798‑M1/1‑166 (MQ=255)
                                                                                   tGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtcc                                                                                                                                                                    >  1:597182‑M1/1‑164 (MQ=255)
                                                                                   tGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtcc                                                                                                                                                                    >  1:144484‑M1/1‑164 (MQ=255)
                                                                                      attattTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgct                                                                                                                                                                 >  2:141699‑M1/1‑161 (MQ=255)
                                                                                         attTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgat                                                                                                                                                              <  2:605687‑M1/251‑94 (MQ=255)
                                                                                            tctctcCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgat                                                                                                                                                           >  1:146388‑M1/1‑155 (MQ=255)
                                                                                                   tcATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgc                                                                                                                                                    >  1:497040‑M1/1‑148 (MQ=255)
                                                                                                    cATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgca                                                                                                                                                   <  2:271331‑M1/251‑105 (MQ=255)
                                                                                                    cATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcg                                                                                                                                                     <  2:146388‑M1/249‑103 (MQ=255)
                                                                                                     ttaGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcggttgg                                                                                                                                                  >  2:368173‑M1/2‑146 (MQ=255)
                                                                                                         aTTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcata                                                                                                                                              <  2:567678‑M1/251‑110 (MQ=255)
                                                                                                         aTTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcat                                                                                                                                               >  1:439799‑M1/1‑142 (MQ=255)
                                                                                                                ccTGTAATGATCCCGCTAAAATGCATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAACAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACATTGCAGAGGGAAGCGTggcatgtacttattcggacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcctaaccgcgtccgctgatgcttgcgcgatgcatatacgggc                                                                                                                                       >  2:549814‑M1/1‑135 (MQ=255)
                                                                                                                ccTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatataccggc                                                                                                                                       >  2:44869‑M1/1‑135 (MQ=255)
                                                                                                                         aTCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctt                                                                                                                              <  1:287393‑M1/251‑126 (MQ=255)
                                                                                                                                       tATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgcaacctttggcctgttcgagt                                                                                                                >  2:475779‑M1/1‑112 (MQ=255)
                                                                                                                                          tCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatctccgcgtccgctgatgattgcgcgcagcatatacaggctgaaaccttgggcatgttcgagttt                                                                                                              <  2:399291‑M1/250‑142 (MQ=255)
                                                                                                                                           cGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttg                                                                                                             >  1:204328‑M1/1‑108 (MQ=255)
                                                                                                                                            gCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTACAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgat                                                                                                           <  1:313940‑M1/251‑145 (MQ=255)
                                                                                                                                             cTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGAATGAAGGCTCAATTTGCAGAGGGAAGCGTggcatgcctttattattacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgctgattccgctcagcatatccaggctgaaatccttggccggttcgagtttgttc                                                                                                          >  2:317643‑M1/1‑106 (MQ=255)
                                                                                                                                               gAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctg                                                                                                        <  1:533667‑M1/251‑148 (MQ=255)
                                                                                                                                                aaaGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgc                                                                                                       >  1:98178‑M1/1‑103 (MQ=255)
                                                                                                                                                aaaGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgc                                                                                                       >  1:534327‑M1/1‑103 (MQ=255)
                                                                                                                                                  aGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgctgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcgg                                                                                                     >  2:288446‑M1/1‑101 (MQ=255)
                                                                                                                                                  aGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcgg                                                                                                     >  2:423124‑M1/1‑101 (MQ=255)
                                                                                                                                                     aGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcaattcatcaacgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtgg                                                                                                  >  1:400514‑M1/1‑98 (MQ=255)
                                                                                                                                                      ggTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtgga                                                                                                 <  1:308062‑M1/251‑155 (MQ=255)
                                                                                                                                                       gTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtgga                                                                                                 >  1:161673‑M1/1‑96 (MQ=255)
                                                                                                                                                               tGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGCAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtggaatggctaac                                                                                        >  2:458488‑M1/1‑88 (MQ=255)
                                                                                                                                                                    cATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtggcatggctaactcttc                                                                                   >  2:51633‑M1/1‑83 (MQ=255)
                                                                                                                                                                            gCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtggaatggctaactcttctttggcg                                                                            >  2:118906‑M1/1‑75 (MQ=255)
                                                                                                                                                                              tAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagcaagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggatgaaacctttggcctgttcgagtttgatctgcggtggaatggctaactcttctttgtcgacc                                                                         >  1:518381‑M1/1‑73 (MQ=255)
                                                                                                                                                                               aaaaGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtggaatggctaactcttctttggcgaccc                                                                        >  1:406346‑M1/1‑72 (MQ=255)
                                                                                                                                                                                aaaGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcctcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtggaatggctaactcttctttggcgaccac                                                                       >  2:495827‑M1/1‑71 (MQ=255)
                                                                                                                                                                                          cTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatataccggctgaaacctttggcctgttcgagtttgatctgcggtggaatggctaactcttctttggcgaccaccacatccacc                                                             >  2:407209‑M1/1‑61 (MQ=255)
                                                                                                                                                                                          cTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagaatatacaggctgaaacctttggcctgttcgagtttgatctgcggtggaatggctaactcttctttggcgaccaccacatccacc                                                             >  1:596066‑M1/1‑61 (MQ=255)
                                                                                                                                                                                            gTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtggaatggctaactcttctttggcgaccaccacatccaccaa                                                           <  1:532316‑M1/251‑193 (MQ=255)
                                                                                                                                                                                               aaaCACAATACAGTGAAAAGTTTTAGACTGAAGGCTCCCTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtggaatggctaactcttctttggcgaccgccacatccaccaacac                                                        <  1:405178‑M1/251‑196 (MQ=255)
                                                                                                                                                                                                aaCACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacaggaaatgccaccctttttaccttagccagtttgttttctccagttcgatcacttcatcaccgcgtccgctgatgattcagagcagcatatatagcctgaaacctttggacttttcgagtttgatctgcggtggaatggctaactcttctttggcgaccaccacctccaccaacacc                                                       >  1:446808‑M1/1‑55 (MQ=255)
                                                                                                                                                                                                   acaATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtggaatggctaactcttctttggcgaccaccacatccaccaacaccgga                                                    <  1:340115‑M1/251‑200 (MQ=255)
                                                                                                                                                                                                      aTACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagattgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtggaatggctaactcttctttggcgaccaccacatccaccaacaccggaccg                                                 >  1:336488‑M1/1‑49 (MQ=255)
                                                                                                                                                                                                       tACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtggaatggctaactcttctttggcgaccaccacatccaccaacaccggaccgt                                                >  1:196744‑M1/1‑48 (MQ=255)
                                                                                                                                                                                                        aCAGGGAAAAGTTTTAGACTGCAGGCTCACTTTGCAGAGGGAAGAGTggcatggccttattatgacgggaaaggcaaccctttttaccttaaccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgatcgagtttgttctgcggtggaagggctcactcttctttcgcgtccaccacatccaccaacaccggaccgac                                               >  2:148621‑M1/1‑47 (MQ=255)
                                                                                                                                                                                                          aGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtggaatggctaactcttctttggcgaccaccacatccaccaacaccggaacgtcg                                              >  1:535418‑M1/1‑45 (MQ=255)
                                                                                                                                                                                                             gAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtggaatggctaactcttctttggcgaccaccacatccaccaacaccggaccgtcgatgg                                          >  1:16845‑M1/1‑42 (MQ=255)
                                                                                                                                                                                                              aaaaGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtggaatggctaactcttctttggcgaccaccaaatccaccaaccccggaccgtcgacgg                                          >  2:141029‑M1/1‑41 (MQ=255)
                                                                                                                                                                                                                 aGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtggaatggctaactcttctttggcgaccaccccatccaccaacaccggaccgtcgctggaga                                       >  2:156916‑M1/1‑38 (MQ=255)
                                                                                                                                                                                                                 aGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtggaatggctaactcttctttggcgaccaccacatccaccaacaccggaccgtcgatggagaa                                      <  1:559204‑M1/251‑214 (MQ=255)
                                                                                                                                                                                                                 aGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtggaatggctaactcttctttggcgaccaccacatccaccaacaccggaccgtcgatggaga                                       >  2:168016‑M1/1‑38 (MQ=255)
                                                                                                                                                                                                                      tAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtggaatggctaactcttctttggcgaccaccacatccaccaacaccggaccgtcgatggagaaggcg                                  <  1:624798‑M1/250‑218 (MQ=255)
                                                                                                                                                                                                                       aGACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtggaatggctaactcttctttggcgaccaccacatccaccaacaccggaccgtcgatggagaaggcgc                                 <  1:156916‑M1/250‑219 (MQ=255)
                                                                                                                                                                                                                        gACTGAAGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtggaatggctaactcttctttggcgaccaccacatccaccaacaccggaccgtcgatggagaaggcgcgt                               >  1:155740‑M1/1‑31 (MQ=255)
                                                                                                                                                                                                                         aCTGAAGGCTCACTATTCAGATGGAAGCGTggcagttccgtattatgacggtaaatgccaccctttttaccttagccagtttgttgtcgccaggtcgatcacgtcatctccgcgtccgctgatgattgcgcgcagcatatacaggctgaatcttttggcctgttcgagtttgatctgcgctgtaatggctatctcttctttggcgaccaccacatccacctacaccggaccttcgatggagaaggcgcgtt                              <  2:564537‑M1/251‑222 (MQ=255)
                                                                                                                                                                                                                              aGGCTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtggaatggctaactcttctttggcgaccaccacatccaccaacaccggaccgtcgatggagaaggcgcgttgcagg                         >  2:500639‑M1/1‑25 (MQ=255)
                                                                                                                                                                                                                                 cTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtggaatggctaactcttctttggcgaccaccacatccaccaacaccggaccgtcgatggcgcaggcgcgttgcagggct                      >  2:531705‑M1/1‑22 (MQ=255)
                                                                                                                                                                                                                                 cTCACTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtggaatggctaactcttctttggcgaccaccacatccaccaacaccggaccgtcgatggagaaggcgcgttgcagggct                      >  1:331226‑M1/1‑22 (MQ=255)
                                                                                                                                                                                                                                    aCTTTGCAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtggaatggctaactcttctttggcgaccaccacatccaccaacaccggaccgtcgatggagaaggcgcgttgcagggcttc                    >  2:564734‑M1/1‑19 (MQ=255)
                                                                                                                                                                                                                                          cAGAGGGAAGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtggaatggctaactcttctttggcgaccaacacatccaccaacaccggaccgtcgatggggaaggcgcgttgcagggcttcatcaact             >  2:567856‑M1/1‑13 (MQ=255)
                                                                                                                                                                                                                                                 aaGCGTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaacctttggcctgttcgagtttgatctgcggtggaatggctaactcttctttggcgaccaccacatccaccaacaccggaccgtcgatggagaaggcgcgttgcagggcttcatcaacttcagacg      <  1:53307‑M1/251‑246 (MQ=255)
                                                                                                                                                                                                                                                     gTggcatgtccttattatgacgggaaatgccaccctttttaccttagccagtttgttttcgccagttcgatcacttcatcaccgcgtccgctgatgattgcgcgcagcatatacaggctgaaaccttttgcctgttcgcgtttgatctgcggtggaatggctaactcttcttttgcgaccaccacatccaccaacaccggaccgtcgagggcgaaggcgcgttgcagggcttcatcaacttcaaacgcttt  >  2:385114‑M1/1‑2 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                         
CATCGCGTTCATTGCTCAATTTCTCAGCCAGATCTTTTGCACGCACACAGTAAGGGCAACCCGAACGACCAAAAATAACGGTTTGCATTATTTCTCTCCTCATAGATTTATGCCTGTAATGATCACGCTAAAATGTATTCGCTGAAAGTAGGTTTAACCTGTTGCATTAATTGCTAAAAGCTATAACTGTTAAACACAATACAGTGAAAAGTTTTAGACTGAAGGCTCACTTTGCAGAGGGAAGCGTATGCGCGCGATCGGTAAATTGCCTAAAGGCGTGTTGATACTGGAATTTATCGGAATGATGCTACTGGCGGTGGCGCTGCTGTCGGTAAGCGACTCCCTGTCGCTGCCTGAGCCATTTTCTCGGCCAGAAGTGCAGATTCTGATGATTTTTCTCGGTGTTTTGCTCATGCTTCCCGCTGCGGTGGTGGTTATTCTTCAGGTGGCAAAACGTCTTGCCCCACAGCTGATGAACCGTCCACCGCAATATTCA  >  NC_000913/890666‑891161

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