Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 69,672 +T coding (377/1701 nt) araB ← ribulokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091369,6721.T98.9% 296.1 / NA 92E126A (GAA→GCA) araBribulokinase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  major base T (58/32);  minor base G (0/1);  total (58/34)
Fisher's exact test for biased strand distribution p-value = 3.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGC  >  NC_000913/69427‑69916
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        
aTTTGTGCCCGGCGCTGCAACGTCCGGTACAAATACCCTGGGGGCGGGGGGTCCCGGAAAGCAGGGCGGGCCCCCAGTCGCACAGCTCAATCCACGATGCGGCAGTTTGCGCCCGGGCGCTGTCCTGGCGAGTCACATGCAGGTTTTTTCCCCAGAACAATTCGGTGGAATAAATCCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAAcc                                                                                                                                                                                                                                                     <  1:385782/251‑1 (MQ=255)
  ttaggCCCGGCGCGCCAACGCCCGCGACGAATACCCGGCGGGCGGGTGGTACCGAAAAGCAGAGCGGGCCCCCGGTCGCACAGCCCAACCCAGGAGGCGGAAGATGGCGCCACGGCGCTGTCCTGGCGAGTCACAGGCAGGATTTTTGCCCAGAACCATCCGCGGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCg                                                                                                                                                                                                                                                    <  1:439804/246‑1 (MQ=255)
  ttATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATATCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGc                                                                                                                                                                                                                                                   >  1:213160/1‑251 (MQ=255)
       cccGGCGCTGCACCGTCCGCGACGAATTTCCGGCGGGCGGGTGGTACGGAAAAGCAGAGGTGCCCCCCAGCCCCCCAGCTCAATCCCCGAGGCGCCAGATGGCGCCACGGCGCGGTCCTGGCGAGTCACATTCAGGATTTTTCCCCAAACCCATCCGCTGAAAAAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAgtgt                                                                                                                                                                                                                                              <  1:812523/251‑1 (MQ=255)
          ggCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGt                                                                                                                                                                                                                                          >  1:411179/1‑251 (MQ=255)
          ggCGCTGCAACGCCCGCGACAAATTCCCGGGGGGCGGGTGGTACGGAAAAGAAGAGCTGGCACCATGTCGCACAGCACAACCCACGTTGCGGCAGATTGCCCCCGGGCGCTGTCCGGGCGAGCCACAGGCAGGTTTTTGGCCCAGAACCATTCGCGGAAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGt                                                                                                                                                                                                                                          <  1:746724/251‑1 (MQ=255)
           gcgcTGCAAGGTCCGCGACGAATACCCGGGGGGGGGGTGGTCCCGGAAAGCAGAGCGGGCCCCCAGTCGAACAGCTAAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCAAGAACCATTCGCTGAAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTc                                                                                                                                                                                                                                         <  1:963112/251‑1 (MQ=255)
            cgcTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCt                                                                                                                                                                                                                                        >  2:422019/1‑251 (MQ=255)
              cTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCACCCGCAGTGTGGTCttt                                                                                                                                                                                                                                      >  2:399096/1‑251 (MQ=255)
                   cggccggcgacaattccccggggggcggtTGGTGCCGGAAGGCGGAGCTGGACCCCGTGCGCACATCTCAATCCCCGATGCGGAAGATGGCGCCACGGCGCTGTCGGGGCTAATCACTTGCAAGATTTGTGCCCAAACCCTTTCGCGGAAAAAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTGCGCTTCTTTCAACCGCAGTGTGGTCTTTCcaca                                                                                                                                                                                                                                 <  2:714051/229‑1 (MQ=255)
                   aCGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCcac                                                                                                                                                                                                                                  >  1:424484/1‑250 (MQ=255)
                   aCGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGAAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCcac                                                                                                                                                                                                                                  >  1:296258/1‑250 (MQ=255)
                     gTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGCTTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCACCCGCAGTGTGGTCTTTCCACaa                                                                                                                                                                                                                                >  2:398015/1‑250 (MQ=255)
                     gTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCAGGCGCGTGGCACAACCGGGAAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCcaaaa                                                                                                                                                                                                                                >  1:784929/1‑247 (MQ=255)
                          ggacgaATAGCCTGCGGGCGGGTTGTCCCGAAAAGCAGTGCTGCCACCCCGCCGCCCTTCTAAATCCATGATCCGGCAGATTGCGCCACGGCGTTGTCCGGGCAAGTCACATGCAGTATTTTTCCCCATAACCATTCGCTTGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGaa                                                                                                                                                                                                                         <  1:679925/250‑1 (MQ=255)
                            acgaATACCCGGCGGGGGGGTGATACCGGACAGCAGAGCGGGCCCCCAGGCGCACAGCCAAACCCACGATGCGGCAGATTGCGCCCCGGCGCGGTCCTGGCGAGTCACATTGAGGATTTTTGCCGAGACCCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACa                                                                                                                                                                                                                       <  1:377864/251‑1 (MQ=255)
                              gaATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATc                                                                                                                                                                                                                     >  1:975851/1‑251 (MQ=255)
                                   accTGCGGGGGGGTGGTACCGGAAAACAGGGCGGGCCCCCCGGCGCAAACCACTACCCCCGATGGGGCAGATTGCCCCCCGGCTCTGCCCGGGCGAGTCACATTCAGGTTTTTTGCCCAAAACCATTCGCTGGAAAAAATCCCGCCAATTTAGCGGGAGTAGTCAACATTGCCCGGCGCCTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGt                                                                                                                                                                                                                 <  1:706884/249‑1 (MQ=255)
                                     cTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGCATAAATACCGCCAATATAGCGGGCGTAGGCAACATTGCCCGGCGCATGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGGGTTGTCTTTCCACAATACGATCACCGCGTTcc                                                                                                                                                                                                              >  2:235364/1‑250 (MQ=255)
                                        cgggcgggTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCATGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCCCTTCGCTTCTTTCAACCGCAGAGTGGGCTTTCCACAATACGAACGTCGCGTTCGGGt                                                                                                                                                                                                           >  1:565438/1‑251 (MQ=255)
                                          ggggggTGGAACCGGAAAGCAGGGGTGGCACCCAGGCGCACAGCTCAAGCCACGATGCGGCAGATTGCGCCCCGGCGCTGTCCTGGCGAGTCACAGGCAGGATTTTGGCCCAGAACCATTCGCTGAAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGttt                                                                                                                                                                                                         <  1:412543/248‑1 (MQ=255)
                                           gcgggTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATCCGAACATCGCGTTCGGGtttt                                                                                                                                                                                                        >  2:953606/1‑251 (MQ=255)
                                               gtggtACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGc                                                                                                                                                                                                    >  1:629577/1‑251 (MQ=255)
                                                 ggtACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCa                                                                                                                                                                                                   >  1:667874/1‑250 (MQ=255)
                                                        gAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCgg                                                                                                                                                                                           >  2:293988/1‑251 (MQ=255)
                                                         aaaGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCggg                                                                                                                                                                                          >  1:295173/1‑251 (MQ=255)
                                                            gCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGCACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATCCGAACATCGCGTTCGGGTTTTCGGCAAAATCCGTgcgc                                                                                                                                                                                       >  2:763796/1‑251 (MQ=255)
                                                              agagCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGgcgc                                                                                                                                                                                       >  1:993362/1‑249 (MQ=255)
                                                                   tGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGcc                                                                                                                                                                                 >  1:548974/1‑250 (MQ=255)
                                                                        cccAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACg                                                                                                                                                                           >  1:690433/1‑251 (MQ=255)
                                                                                 cacaGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTAGCTTCTTTCAACCGCAGTGTGGTCTTTCCACACTACGAACATCGCGTTCGGGTTTTCGGCACACTCCGGGCGCAGCGCCAGCACGTTACCGTc                                                                                                                                                                   >  1:147477/1‑251 (MQ=255)
                                                                                  acaGCCCAATCCACGAGGCGGCAGATGGCCCCCCGGCGCTGTCCGGGCGAGTCACATGCAGGATTTTTGCCAAAACCCATCCGCGGAAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCg                                                                                                                                                                  <  1:293988/251‑1 (MQ=255)
                                                                                           tCCACGATGCGGCAGATTGAGACACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATACATACCGCCAATATAGCGGGCGTAGTCAACATTGCCCGGCGCGTGGCACAAAAGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGGTTTCGGCAAACTACGGGCGCAGCGACAGCACggtaccggaggcctccatc                                                                                                                                                         >  2:897329/1‑233 (MQ=255)
                                                                                                gATGCGGCAGAATGAGAAACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATGCGCTGGAATAAATACAGCCGATATAGCGGGAGTAGTCAACATTGCCCGGGGCGAGACACAAACGGTTCATCTATTCGCATCTTTCAACAGAAGTGTGCACTTTTCCCAAAACGAAAATCGTGTTAGGTGTTTCGTGAAACTATGGACGCAGCGCCAGAACGTTTACGTCGGGCTCACTCGGTGc                                                                                                                                                    >  2:745035/1‑251 (MQ=255)
                                                                                                      gCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACAGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCATCACGTTACCGTCGGCATCAATCGGTGCGGGCgg                                                                                                                                              >  2:530673/1‑250 (MQ=255)
                                                                                                       cAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCTTCAATCGGTGCGGGCGTc                                                                                                                                            >  2:93996/1‑251 (MQ=255)
                                                                                                              ggccgcGGCGTTGTCGTGGCTGGTCACAGGCAGGTTTTTTGCCCAGACCCATCCGCGGAAATAAATCCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCGCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCgg                                                                                                                                     <  1:789772/246‑1 (MQ=255)
                                                                                                                      cgcTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGAGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATAAATCGGTGCGGGCGTCGAGCCGGTACTGTc                                                                                                                              >  1:503330/1‑250 (MQ=255)
                                                                                                                       gcTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGAGTGGCACAAACGGGTAATCTCTTCGCGTCTTTCAACAGCAGGGTGGGATTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCAAGTTACCGTGGGAATCAATCGGGCCGGGCGTCTAGCCGGTCCTggca                                                                                                                             >  2:209287/1‑247 (MQ=255)
                                                                                                                          gTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACAGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACg                                                                                                                         <  1:530673/250‑1 (MQ=255)
                                                                                                                              tggcggggcacaTGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCGATATAGCGGGAGTAGACAACCTTGCCCGGCGCGTGGAACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGGGTTCGGGTTTTCGGCAAACTCAGGGCGCAGCGCAAGCAGGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGAACTGTCAACGcgca                                                                                                                     >  2:41344/9‑247 (MQ=255)
                                                                                                                                 cGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATccc                                                                                                                 >  1:887059/1‑251 (MQ=255)
                                                                                                                                        aaTGCAGGATTTTTGCCCAGAACCATTCGCTGGAATCAATACCGCCAATCTAGCGGGAATAGTCAACATTGCCCGGCGCGAGGCACAAACGGGTCATCTCTTCGCATCTTTAAACCGCAGTGTGGTCTTTCCCCAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCCTCAATCGGTGCGTGCGTCGAGCCGGAACAGTCAACTCCAATCCCGCCcac                                                                                                           >  2:849121/2‑250 (MQ=255)
                                                                                                                                          tGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCAGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGc                                                                                                        >  1:957688/1‑251 (MQ=255)
                                                                                                                                             aGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTgc                                                                                                     >  1:787385/1‑251 (MQ=255)
                                                                                                                                                  ttttGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCt                                                                                                >  1:647870/1‑251 (MQ=255)
                                                                                                                                                      tcccaaaacctgtcgcggaaaaaaataCCGCCAATATAGCGGGAGTAGGCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGTCGTCGAGCCGGTACTGGCAAAGCGAATCCCGACCACAGCTGCGCGCTGtt                                                                                             <  1:803287/228‑1 (MQ=255)
                                                                                                                                                            aaCCATTCCCTGGACTAAAAACCGCCAATATAGCGGGAGTAGTCAACATTCACCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTACACAATCCGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGCCGCAGCGCCAGCCCGGTACCGTCGGCATCCACAGGCGAGGGCGTAGCGCCGGTACTGTACACACCAATCCCGCCcaacgatgcgcgcagtacgaagcc                                                                                      >  2:458370/1‑229 (MQ=255)
                                                                                                                                                              ccttctgttgaaaaaaataCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTaa                                                                                    <  1:978456/237‑1 (MQ=255)
                                                                                                                                                              ccATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTa                                                                                     >  2:756455/1‑250 (MQ=255)
                                                                                                                                                                aTTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCAAGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGCTCGACGCTAAGc                                                                                  >  1:745034/1‑251 (MQ=255)
                                                                                                                                                                             aatcccgccaaaatatcGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGTTAATTCCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCTCAGCGCCAGCACGTTACCGTCGGCATCAATCGTTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGTCGCTAAGCTCTGCAAGCACgg                                                                     <  2:754382/235‑1 (MQ=255)
                                                                                                                                                                             aatcccgccaATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACgg                                                                     <  1:347679/244‑1 (MQ=255)
                                                                                                                                                                             aataCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACgg                                                                     <  1:668023/251‑1 (MQ=255)
                                                                                                                                                                                accgccaATATAGCGGGCGGAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATAGCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCCCAAGACGAACATCGGTTTCGGATTTTCGGAAAACTCGGTGCGCAGCCCAAGCACGGTACAGTCGCCTTCAATCGACGCGGCGATCGAGCCGGAACTGCACACGCCAGTCCCTACCACATCGGCGCCCGGTTCGCGGCTAAGCTCAGCAAACACGGttt                                                                  >  2:587592/5‑251 (MQ=255)
                                                                                                                                                                                 ccgccaATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGtttt                                                                 >  1:691152/4‑251 (MQ=255)
                                                                                                                                                                                  cgccaATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTc                                                                >  1:270450/3‑251 (MQ=255)
                                                                                                                                                                                      aagagagcgggggtatccAACATTGCCCGGCGCGGTGCCCAAACGGGTAATCTCTTCGCGTCTTTCAACCGCAGGGTGGTCTTTCCACAATCCGAACACCGCGTGCTGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCCCGTTACCGTCGGCATAAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGTCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCGCTGCAAGCACGGTTTTCAGTg                                                            <  2:975851/234‑1 (MQ=255)
                                                                                                                                                                                       atatAGCGGGAGTAGGCAACATTGCCCGGCGCGTGGCAAAAACGGGTAATCTCTTCGTTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGc                                                           <  1:552/248‑1 (MQ=255)
                                                                                                                                                                                        tatAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGcc                                                          >  1:30865/3‑251 (MQ=255)
                                                                                                                                                                                         atAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCg                                                         <  1:728413/250‑1 (MQ=255)
                                                                                                                                                                                         atAGCGGGAGTAGACAACATTGCCCGTCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACGCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGcc                                                          <  2:667874/249‑1 (MQ=255)
                                                                                                                                                                                         atAGCGGGAGTAGACAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCACCACCAATCCCGACCACCGCTGCGCGCTGTTCGACGCTAAGCTCTGCTAACACGGTTTTCAGTGccc                                                         >  2:931220/2‑250 (MQ=255)
                                                                                                                                                                                         agagggggAGTAGTCAAAATTGGCCGGGGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCGTTCCACAATACGAACCTCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCGGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCg                                                         <  2:760786/246‑1 (MQ=255)
                                                                                                                                                                                            gCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCtt                                                      >  1:845617/1‑251 (MQ=255)
                                                                                                                                                                                              gggAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGAAATCTCTTCGCTTCTTTGAACCGCAGTGTGGTCTATCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGATACGGGCGTCGAGCCGGAACTGTCAACGACAATCACGACCACAGCTGCGCGCTATTCGACGCTAAGCTCTGCCAGCACGGTTTTCGGTGCCGCttac                                                    >  1:919566/1‑249 (MQ=255)
                                                                                                                                                                                              gggAGTAGTCAACATTGCCAGGAGCGGGGACCAAACGGGTGATCTCTTAGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGGTACCGTCGGCATCAATCGGTACGGGCGTCGAGCCGGTACTGTCAACGCCAGTCCCGCCCACAGCAGCGAGCTGTTCGACGCTAAGCTCTGCCAACACGGTTTTCAGAGCCGCTTcc                                                    >  2:181629/1‑251 (MQ=255)
                                                                                                                                                                                               ggAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGCACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGAATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCCAGGTGTTCAGTGCCGCTTccc                                                   >  2:814803/1‑250 (MQ=255)
                                                                                                                                                                                                gAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCGGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCAt                                                  >  2:265623/1‑251 (MQ=255)
                                                                                                                                                                                                  gtagtCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTg                                                <  1:771548/251‑1 (MQ=255)
                                                                                                                                                                                                     gtCAACATTGACCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGATTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGGCGAGCCGGTACTGTCACCGCCAATCCCGACCACAGCTGCGCGCGGTTCGCCGCTAAGCTCTGCAAGCAAGGCTTTCAGTGCCGCTTCAATTGAc                                              >  2:780488/1‑250 (MQ=255)
                                                                                                                                                                                                      tCAACATTGACCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGACCATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCACCGCCAATCCCGACCACAGATGCGCGCAGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCCATGCCGCTTCCCTTGAccc                                            >  2:645120/1‑249 (MQ=255)
                                                                                                                                                                                                      tCAACATTGACCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGCCCCCAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCACATCCATTGtctc                                            >  2:891180/1‑251 (MQ=255)
                                                                                                                                                                                                          cattGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGCGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATg                                        >  2:285799/3‑251 (MQ=255)
                                                                                                                                                                                                           attGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGt                                       <  1:397100/250‑1 (MQ=255)
                                                                                                                                                                                                                 cGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCAc                                 <  1:648851/251‑1 (MQ=255)
                                                                                                                                                                                                                 cGGCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCAc                                 <  1:422019/251‑1 (MQ=255)
                                                                                                                                                                                                                  ggCGCGTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCAcg                                >  1:788428/1‑251 (MQ=255)
                                                                                                                                                                                                                   cagcgTGGCACAAACGGGTAAGCTCTGTGCGTCCTTCAACCGCAGGGTGGTCTTTCCAAAATCACACCATCGCGTGCGTGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAAGCGGTGCGGGCGCCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGATTCCATTGCCTCAATGTAGTCAcgc                               <  2:908698/249‑1 (MQ=255)
                                                                                                                                                                                                                    cgcgTGGTACAAACGGGTAATCTATTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTAACGTCGGCATCAATCGTTGCGGGCGTCGCGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCCCGCCAGGTTTTAAGTGCAACTTCCATTGACTCAATGTAGTCAcgcg                              >  2:984034/1‑251 (MQ=255)
                                                                                                                                                                                                                    caggtgtcataaaCGGGTAACCTCTTCACTTCTTTCAACCGCAGTGTTGGTTGTCCACAATACGAACTCCGCGTGCGGGTTTTCGGCAAGCTCCGGGCGCAGCGCCAGCACGTTCCCGTCGGCATCAATAGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACACAATGTAGTCAcgcg                              <  2:8073/241‑1 (MQ=255)
                                                                                                                                                                                                                      cgTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCCTCCATTGACTCAATGTAGTCACGCgg                             >  2:667881/1‑250 (MQ=255)
                                                                                                                                                                                                                       gTGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGgat                           >  2:876233/1‑251 (MQ=255)
                                                                                                                                                                                                                       gTGGCACAAACGGGTAATCTCTTCACTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCTCGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGgat                           <  1:235364/251‑1 (MQ=255)
                                                                                                                                                                                                                        tGGCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGgatg                          >  1:490262/1‑251 (MQ=255)
                                                                                                                                                                                                                        aggCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCCATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCAGCAAGCACGGTTTTCATTGCCGCTTCCATTGACTCAATGTAGTAACGCggttg                          >  2:220480/2‑248 (MQ=255)
                                                                                                                                                                                                                         ggCACAAACGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGgatga                         <  1:399096/251‑1 (MQ=255)
                                                                                                                                                                                                                                aCGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATGGACTCAATGTAGTCACGCGGATGATGACGGa                  <  1:908508/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cGGGTAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGaa                 <  1:459766/251‑1 (MQ=255)
                                                                                                                                                                                                                                     tAATCTCTTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGAAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTgg             >  2:409896/1‑251 (MQ=255)
                                                                                                                                                                                                                                         ctctTCGCTTCTTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTAt         >  1:379193/1‑251 (MQ=255)
                                                                                                                                                                                                                                               tcGCTTCTTGCAACCGCAGTGTGGTCTTTCCCCAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTc       <  2:295173/248‑1 (MQ=255)
                                                                                                                                                                                                                                                   tccttcAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTAAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGc  <  1:598789/247‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        
ATTTATGCCCGGCGCTGCAACGTCCGCGACGAATATCCTGCGGGCGGGTGGTACCGGAAAGCAGAGCTGGCACCCAGTCGCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCTTCAACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGC  >  NC_000913/69427‑69916

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