Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 71,824 A→T P158P (CCA→CCT yabI → DedA family protein YabI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091371,8240AT96.0% 208.5 / 2.0 76P158P (CCA→CCTyabIDedA family protein YabI
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base T (30/43);  total (31/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.08e-01

GGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGT  >  NC_000913/71582‑72068
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
ggctGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAATCAATACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGTCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCt                                                                                                                                                                                                                                                <  1:398158/248‑1 (MQ=255)
ggctGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCt                                                                                                                                                                                                                                                <  1:971861/248‑1 (MQ=255)
 gTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTg                                                                                                                                                                                                                                               <  1:308489/251‑1 (MQ=255)
  gcaaccatttAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGGTCTGCCGGTGGCACCATTCATTACACCGAATATTATCGGCTGCCCGCTGTGGCCCCCGTTTTACTTCCTGCCTGGGATTCTgg                                                                                                                                                                                                                                              >  2:505415/8‑251 (MQ=255)
        ttttAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTggcgggcg                                                                                                                                                                                                                                       >  1:336791/1‑251 (MQ=255)
         tttAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTggcgggcg                                                                                                                                                                                                                                       <  1:691650/250‑1 (MQ=255)
         tttAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGgcgc                                                                                                                                                                                                                                      >  1:715484/1‑251 (MQ=255)
              aaaGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGGAGTCTGGATCTGCCGGTGGCAAAATTCATTACCCCGAATATTATCGGCCGCCTGCTGTGGCCCCCGTTTTACTTCCTGCCTGGGagtcgtgcgggggcgggcg                                                                                                                                                                                                                                  >  2:356411/1‑233 (MQ=255)
               aaGCCGTTGCATCGCTGGTCATTTCTGTAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCCCCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGcgat                                                                                                                                                                                                                                <  2:677797/251‑1 (MQ=255)
                 ggcGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCTGGCGCGGcgatcg                                                                                                                                                                                                                              <  1:585811/249‑1 (MQ=255)
                         atcccTGGTCATTTTTTAAGAAAAACAAAGCATTACCCGAGAAAACTGAACATGCGTTGCACTACCACAGCATGTTCACCATCCGGGTCGGTCGCTGTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTAGCGGCTGCCTGCTGTGGCCTCCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTcca                                                                                                                                                                                                                     <  2:592820/247‑2 (MQ=255)
                           cGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGGACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGACCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGCTGGCAAAATTCCTTACGCCGAATATTATCGGCTGCCTGCTGTGGCCCCCGTTTTACTTCCTGCCTGGGATTCTGGGGGGCGCGGGGATCGATAtttccgc                                                                                                                                                                                                                   >  2:132311/1‑246 (MQ=255)
                            acTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGACGCTTTGTTGGCCCGACGCGTCCGCTGGTGCAAATGGTGGCGGGAATGCTGGATATGACGGTAGCAAAATTCATTACGCCGCATATTATAGGCTGCCTTCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTggcggggctgggggtcgatcattcagcc                                                                                                                                                                                                                    >  2:508318/2‑229 (MQ=255)
                               gggcATTTCGGAAGAAAAAAAAAGCATTACTCGATAAAACGGAACAGGCGTTGCATCAACACAGCATGGTCACCATTCGGGGCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGGGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCg                                                                                                                                                                                                                 <  2:285485/246‑1 (MQ=255)
                               ggTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGGTATTCCCGCCggg                                                                                                                                                                                                               >  1:114978/1‑250 (MQ=255)
                                  cATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGGTCGATATTCCAGCCGGAATg                                                                                                                                                                                                            >  1:232588/1‑251 (MQ=255)
                                   aTTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGTATGTTGGGTCTGCCGGTGGCCACATTCGTTACGCCGGATATTATCGGCGTCCTGCTGTGGCCACCGTTTTAACTCCTTCCTGGGCATCTGGCGGGCGGGGcgagcgggttttcagccgggaagg                                                                                                                                                                                                            >  2:16574/1‑230 (MQ=255)
                                   aTTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGGGCGATATTCCAACCGGAATgg                                                                                                                                                                                                           >  1:648349/1‑250 (MQ=255)
                                     ttCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGGTCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTTCCTGGGATTCTGGCGGGCGCGGCGATCGGTATTCCAGCCGGGATGc                                                                                                                                                                                                           >  1:635148/1‑249 (MQ=255)
                                      cccgccgcccttccgatctttactCGATAAAACTGAACATGCGTTGCATCAACCCAGCATGTTCCCCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCaga                                                                                                                                                                                                        <  2:503087/231‑1 (MQ=255)
                                          aagaaAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCgg                                                                                                                                                                                                    <  1:572563/251‑1 (MQ=255)
                                            gaaAAACAAAGCATTACGCGATAAAACTTAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGGCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCggcg                                                                                                                                                                                                  <  1:132311/251‑3 (MQ=255)
                                             aaaaaCAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGGTTCTGGCGGGCGCGGCGCTCGCTATTTCAGCCGGAATGCAGAGCggcgg                                                                                                                                                                                                 >  1:674943/1‑248 (MQ=255)
                                             aaaaaCAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGGGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGGTATTCCAGCCGGAATGCAGAGCggggg                                                                                                                                                                                                 >  1:891673/1‑248 (MQ=255)
                                              aaaaCAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGCCGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGAGTCTGGCGGGCGCGGCGGTCGAATTTTCAGCCGGAAATCAGAGCggcggg                                                                                                                                                                                                >  1:97749/1‑247 (MQ=255)
                                                  cAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCggcgcgttt                                                                                                                                                                                             >  1:983021/1‑243 (MQ=255)
                                                       cattactCGATAAAACTGAACATGCGATTCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTACGCTGGTGCCAATGGTGGCGGGAATGCCGGAACTGCCGGTGGCAAAATTCATTACGCCGAATAATATCGGATGCCCTCTATGGCCACCGTTTTAATTCCTACCTGGGATTCTGGCGGGCGCGGAGATCGAATATCCAGCCGGAATGAAGAGCgggaggattaaatgg                                                                                                                                                                                       >  1:917214/4‑238 (MQ=255)
                                                        attactCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGCATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGGGGGGGGGTCGATATTCCCGCCCGAATgacgaggggcgaggtttaatggt                                                                                                                                                                                      >  2:307788/3‑229 (MQ=255)
                                                         ttactCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCCGCCGGAATGCaga                                                                                                                                                                                                        >  1:503087/2‑232 (MQ=255)
                                                          tacccgATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGtt                                                                                                                                                                                     <  1:920128/245‑1 (MQ=255)
                                                           actcgATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACACCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTtgc                                                                                                                                                                                   <  1:505415/247‑1 (MQ=255)
                                                             tcgATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCCGCCGGAATGCCGAGCGGCGGGTTTAAATgggttctg                                                                                                                                                                                 >  1:677428/3‑245 (MQ=255)
                                                                aTAAACCTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGctgctg                                                                                                                                                                              <  2:336791/251‑1 (MQ=255)
                                                                aTAAAACTGAACATGCGTTGCATCAACACAGCATTTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGctgctg                                                                                                                                                                              <  1:191263/251‑1 (MQ=255)
                                                                 aaaaactaaaccggcgctGCATCCACACTGCAGGTGCACCATTCGGGGCGGGCGCTTTGTTGGCCCGACGCTGCCGCGGGTGCCAATGTGGGCGGGAAGGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTgg                                                                                                                                                                             <  2:388327/234‑1 (MQ=255)
                                                                    aaCTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCCGCCGGAATGCAGAGCGGGGGGTTTAAATGGTTGCTGCTGGcgg                                                                                                                                                                          >  1:708062/1‑249 (MQ=255)
                                                                        gAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATg                                                                                                                                                                                                            <  1:630070/213‑1 (MQ=255)
                                                                         aaCATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGGGGGCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCCGCCGGAATGAAGAGCGGCGGGGTTAAATGGTTGGTGCTGGcggaggcg                                                                                                                                                                     >  1:142919/1‑244 (MQ=255)
                                                                            aTGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTAGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGATGGGTATGCCGGTGGCAAAATTAATTAAGCCGAATATTATCGGCTGACTGGTGTGGACACCGTTTTAATTCCCGCCTGGGGTTCTGGGGGGGCGCGGGGTCGATATTCCAGGCGGCATCCAGAGCGGGGGGTTTTAATGGTtgctggcggggacggaggc                                                                                                                                                                   >  1:415293/1‑236 (MQ=255)
                                                                              gctttGCATCAACCCAGCAGTTTCCCCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGgtat                                                                                                                                                                 <  2:114978/247‑3 (MQ=255)
                                                                              gCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGAGATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGgtatt                                                                                                                                                                <  1:637119/251‑4 (MQ=255)
                                                                                  tGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGCATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATTTTCCAGCCGGAATGCAGAGCGGCGAGTTTTAATGGTTGCTGCTggggacggcggtatttgtg                                                                                                                                                            >  1:81655/1‑234 (MQ=255)
                                                                                             aGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGGGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGGATGCAAAGCGGCGAGTTTAAATGGTTTCTTACGGCGACGGCGGTATTTGTGTGGGTttgtgg                                                                                                                                                 >  1:611649/1‑251 (MQ=255)
                                                                                                   ttCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGGTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTTGTggctggct                                                                                                                                           >  1:269377/1‑251 (MQ=255)
                                                                                                           tCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTccag                                                                                                                                                                                                                    <  1:651795/170‑3 (MQ=255)
                                                                                                               gacggtcgCTTTGTGGGCCCGCCGCGTCCGCGGGTGCCAATGGTGGCGGGAATGGTGGATCTGCCGGTGGCAAAATTCATTGCGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCGGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGtt                                                                                                                               <  2:203310/249‑1 (MQ=255)
                                                                                                                 cggtcgCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGGTCGATATTCCAGCCGGGATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGGCGGCGGGATTTGTGTGGGGTGGTGGCTGGCTGTGCTGGCgggttt                                                                                                                             >  1:603635/1‑247 (MQ=255)
                                                                                                                  ggtcgCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGGATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGGATTTGTGTGGGTTGGGGGGTGGCTGTGCTGGCgggtgtt                                                                                                                            >  1:230587/1‑246 (MQ=255)
                                                                                                                        tttGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCCGCCGGGATGCAGAGCGGCGGGTTTAAATGGTTGCTGCTGGCGACGGCGGGATTTTGGTGGGGTGGGGGCTGGCTGtggtgggggggtttgggcgg                                                                                                                      >  1:117252/1‑237 (MQ=255)
                                                                                                                        tttGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGCGTTTAAATGGTTGCTGCTGGCGACGGCGGGATTTGTGTGGGTTTGTGGCTGGCTGTGCTGGCgggtttggggggg                                                                                                                      >  1:92557/1‑240 (MQ=255)
                                                                                                                          tggtgtcccggcgggtcctcgggTGCAAATGGTTGCGGGATTGTTGGATCTGCCGGTGGCAAAATGCATTACGCCGAATATTATCGTCTGCCGGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCTGCGATTGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGTTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGc                                                                                                                     <  2:64203/229‑1 (MQ=255)
                                                                                                                          ggtttgcccaccgcgTCCGCTGGTGCCAATGGTGGCGGGAAGGCTGGCTCTGCCGGTGTCCAAATTCACTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGGTTGACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGGTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGc                                                                                                                     <  2:394733/239‑1 (MQ=255)
                                                                                                                               gCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTTCCGCCGGGATGCAGAGCGGCGAGGTTAAATGGGTGCTGGTGGCGGCGGCGGTTTTTTTGTGGGGTGGGGGCTGGCGGTGCTTGCgggttttgtggtccgggg                                                                                                                 >  2:307855/1‑233 (MQ=255)
                                                                                                                                      gcgTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGcaaccg                                                                                                        <  1:473692/251‑6 (MQ=255)
                                                                                                                                      gcgTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGcaaccg                                                                                                        <  1:149110/251‑6 (MQ=255)
                                                                                                                                         cccGCTGGTGCCACTGGTGGCGGGAATGCTGGAGCTGCCGGTGGCAAAAGTCATTTCGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGcaaccgatc                                                                                                     <  2:427747/250‑9 (MQ=255)
                                                                                                                                                        ggtggCGGGAATGCTGGATCTGCCGGTGGCAAAAGTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTATTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCattat                                                                                       <  1:993768/250‑1 (MQ=255)
                                                                                                                                                              gggAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCTGTAAGTTTAAATGGTTGCTGCTGGCAACAGCgg                                                                                                                                                                    <  2:118778/167‑1 (MQ=255)
                                                                                                                                                              gggAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCTGTAAGTTTAAATGGTTGCTGCTGGCAACAGCgg                                                                                                                                                                    >  1:118778/1‑167 (MQ=255)
                                                                                                                                                                 aaTGCTGGAGCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTcgcgcg                                                                             <  1:865995/251‑1 (MQ=255)
                                                                                                                                                                  aTGCTGGATCTGCAGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCgg                                                                            <  1:783426/251‑1 (MQ=255)
                                                                                                                                                                        gttctgtgtgtgtcaaaagtcattaCGGCGAATATTATCGGCTGCCGGCTGTGGTATTCGTTTTTCTTCATGACTGCGATTCTGGCGGGCGCGTCGATCGATGGTCCAGCCGGAATTCAGAGTGGCGAGTTGAAATGTTTGCTGATGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGAGGTAATGCAACCGATCGTTTAAGTCATGATTTGTCGCGCGGCCGttt                                                                      <  2:756836/230‑1 (MQ=255)
                                                                                                                                                                             gCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTTGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGTCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTttg                                                                     <  1:242311/247‑1 (MQ=255)
                                                                                                                                                                             gCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTttgt                                                                    <  1:137733/248‑1 (MQ=255)
                                                                                                                                                                              ccGGTGGCAAAATTCATTACGCCGAATATTTTCGGCTGCCTGCTGTGGCCACCGTTTTTCTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTtgtg                                                                <  2:648349/251‑1 (MQ=255)
                                                                                                                                                                               cGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTgg                                                               <  1:4234/251‑1 (MQ=255)
                                                                                                                                                                                   ggcaaaaTTCATTACGACGAAGAGTATCGGCTGGCTGCTGTGGCCATCGTTTTACTGCCGGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTg                                                            <  2:183427/246‑1 (MQ=255)
                                                                                                                                                                                                    cGAATATTATCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGcc                                          <  1:553288/251‑1 (MQ=255)
                                                                                                                                                                                                        gattagCGGCTGCCTGCTGGGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATAGTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACTGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTAGGGCGGAGTGGTAAAGCAACCGATCGTTTAAGTCATTATTTATCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCg                                      <  1:616216/250‑1 (MQ=255)
                                                                                                                                                                                                          ttaccGGCTGCCTTCTGTGGCCACCTTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATTGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCTATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCgg                                     <  2:830481/246‑1 (MQ=255)
                                                                                                                                                                                                           tatCGGCTGCCTGCTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGc                                    <  1:73009/250‑1 (MQ=255)
                                                                                                                                                                                                                   ggctgtggtgGCCTCGGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGgcg                           <  2:81655/244‑1 (MQ=255)
                                                                                                                                                                                                                   gcctgcTGTGGCCCCCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGgc                            <  2:344918/250‑1 (MQ=255)
                                                                                                                                                                                                                   gcctgcTGTGGCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCTGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGtggtgg                             <  1:924854/249‑1 (MQ=255)
                                                                                                                                                                                                                            ggCCACCGTTTTACTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGgtggtg                  <  1:841082/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aCTTCCTGCCTGGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCa       <  2:708062/250‑1 (MQ=255)
                                                                                                                                                                                                                                            cctgcctgGGATTCTGGCGGGCGCGGCGATCGATATTCCAGCCGGAATGCAGAGCGGCGAGTTTAAATGGTTGCTGCTGGCGACGGCGGTATTTGTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCAACCGATCGTTTAAGTCATTATTTGTCGCGCGGCCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGt  <  1:663762/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
GGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCCGTTTTACTTCCTGCCAGGGATTCTGGCGGGCGCGGCGATCGATATTCCTGCCGGAATGCAGAGCGGTGAGTTTAAATGGTTGCTGCTGGCAACAGCGGTGTTTTTGTGGGTTGGTGGCTGGCTGTGCTGGCGGTTATGGCGCAGCGGTAAAGCGACTGACCGTTTGAGTCATTATTTGTCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGT  >  NC_000913/71582‑72068

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