Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,591,074 G→A L418L (CTG→TTG)  ugpB ← sn‑glycerol 3‑phosphate ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,591,0740GA95.7% 294.6 / NA 94L418L (CTG→TTG) ugpBsn‑glycerol 3‑phosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (35/55);  minor base C (4/0);  total (39/55)
Fisher's exact test for biased strand distribution p-value = 2.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAGTGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCG  >  NC_000913/3590827‑3591318
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
aaaaGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGc                                                                                                                                                                                                                                                   >  1:131469/1‑251 (MQ=255)
     aTAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCctgct                                                                                                                                                                                                                                              <  2:432214/251‑1 (MQ=255)
         cGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTgtgg                                                                                                                                                                                                                                          >  2:247091/1‑251 (MQ=255)
         cGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAAATGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATCAGCGAACCGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCACCTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTgtgg                                                                                                                                                                                                                                          >  1:37006/1‑251 (MQ=255)
             gatgatgaGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACCTTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGtc                                                                                                                                                                                                                                      >  1:98570/1‑251 (MQ=255)
              atgatgaGCTGCGTCGCGACCAGAAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTCCACGCGATGCCAGGTCAGATAAGCGAAACGCAACCGACATTCCACTGCTTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCACTGCCTGCTGTGGTGtct                                                                                                                                                                                                                                     >  2:612104/1‑251 (MQ=255)
                   gaGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTAACCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCAGCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTAAACGGCGGTATCCAATGCCTTCTGTGGTGTCTTCTTa                                                                                                                                                                                                                                <  2:130061/251‑1 (MQ=255)
                     gcggcggcgCGACCAGCAGAGAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGAGGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTAcc                                                                                                                                                                                                                              <  2:543265/248‑1 (MQ=255)
                        gcggcgCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCg                                                                                                                                                                                                                             <  1:19730/249‑1 (MQ=255)
                        gcggcgCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCg                                                                                                                                                                                                                             <  2:222292/249‑1 (MQ=255)
                          ggcgCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTGCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTcc                                                                                                                                                                                                                         <  1:251625/251‑1 (MQ=255)
                                aCCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCACTGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGc                                                                                                                                                                                                                   >  2:129171/1‑251 (MQ=255)
                                        aTAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGc                                                                                                                                                                                                           >  1:456098/1‑251 (MQ=255)
                                        aTAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTAGCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGc                                                                                                                                                                                                           <  1:36172/251‑1 (MQ=255)
                                          aaGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGctc                                                                                                                                                                                                         <  1:624926/251‑1 (MQ=255)
                                            ggCAGCCAGCGCGAGCGGACCACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTcc                                                                                                                                                                                                              <  1:226211/244‑1 (MQ=255)
                                            ggCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGAGGAGTCTTCTTACCGGGCCACACGCTCTcc                                                                                                                                                                                                              >  2:226211/1‑244 (MQ=255)
                                              cagccagcGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTTTCCCATGCCTGCTGTGGTGTCTTCTTTCCGGTCCACAAGCCCTCCCGCTCTTc                                                                                                                                                                                                      >  2:12933/1‑250 (MQ=255)
                                              cagccagcGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCa                                                                                                                                                                                                     <  1:183177/251‑1 (MQ=255)
                                                gccagcGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATc                                                                                                                                                                                                   <  2:539608/251‑1 (MQ=255)
                                                 ccagcGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGACCTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATACAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCGCCCGCTCTTCATcc                                                                                                                                                                                                  >  2:121263/1‑251 (MQ=255)
                                                     cgcgAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACaa                                                                                                                                                                                              <  1:197502/251‑1 (MQ=255)
                                                      gctaGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCTACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGTCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAAt                                                                                                                                                                                             <  1:249623/248‑1 (MQ=255)
                                                        gAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGTTGTCTTCTGACCGGTCCACACGCTCTCCATCTCTTCATCCACAATCa                                                                                                                                                                                           <  2:98570/251‑1 (MQ=255)
                                                              aaCACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGa                                                                                                                                                                                     <  2:453098/251‑1 (MQ=255)
                                                                  ccGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTg                                                                                                                                                                                 <  2:189864/251‑1 (MQ=255)
                                                                     gACGGTATGATGTCATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATGTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTTGGGTGGCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCTCGACCACGCGGATCTGCgg                                                                                                                                                                              <  2:268726/251‑1 (MQ=255)
                                                                         ggATGATTACATTCCGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCAACCGACATTACACTGATTAAGACATCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCCGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGCTCCGCGGcctc                                                                                                                                                                          >  2:223410/1‑248 (MQ=255)
                                                                         ggATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGCTGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCCGCTCTTCATCCCCAATCACGCGGATCTGCGGcctg                                                                                                                                                                          >  1:230084/1‑248 (MQ=255)
                                                                            tgatgaCATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTg                                                                                                                                                                       >  1:212362/1‑251 (MQ=255)
                                                                             gatgaCATTACGGTTAATTCCTGAACTGGGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAATCGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGc                                                                                                                                                                      <  2:475416/251‑1 (MQ=255)
                                                                                gaCATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGccc                                                                                                                                                                    <  1:361038/250‑1 (MQ=255)
                                                                                           tAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGccc                                                                                                                                                         <  1:341328/250‑1 (MQ=255)
                                                                                            aaTTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCtt                                                                                                                                                       <  1:172159/251‑1 (MQ=255)
                                                                                               tCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGCTTTCTAAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTgg                                                                                                                                                    >  1:615176/1‑251 (MQ=255)
                                                                                                  tGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTg                                                                                                                                                  >  1:525472/1‑250 (MQ=255)
                                                                                                   gaacggggCCGGATGCGCTTCGCTTATCCGGCCTCCACGCGATGCCAGGTCAGATAAGCGAAACCAATCCGATATTACACTGATTAAGACTTCGTCGATTTCTCAACGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTTCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCAAGCGGATCTGCGGAATGTTGCCCAGACGCAGCCCTTTGGTGaa                                                                                                                                                <  2:458320/244‑1 (MQ=255)
                                                                                                   gAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACTCGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGaa                                                                                                                                                <  1:420092/251‑1 (MQ=255)
                                                                                                   gAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACAATGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGa                                                                                                                                                 >  2:327527/1‑250 (MQ=255)
                                                                                                        gtgCCTGATGCGCTTCGCTAATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTGCCGGTCCACACGCGCTCCATCTCGTCATCCACAATCACGCGGATCTGCGGCATGTTTCCCAGACGCAGCCCTTTGGTGaacggc                                                                                                                                            <  2:58142/250‑1 (MQ=255)
                                                                                                         tgCCGGATGCGCTTCGCTTATCCGGACAAAACGCGATGCCCGGTCAGCTACGCGAAACGCATCCGACATTACCCTGATTAAGACTTCGTCGCTTTCTCAAAGAGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATACAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAAGCTCTCCAGCTCTTCATCCACAATCACGCGGCTCTGCGGCATGTTGCCCAGAAGCCGACCTTTGGTGaaaggcc                                                                                                                                           >  2:575608/1‑249 (MQ=255)
                                                                                                             ggATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGaacggcaacggc                                                                                                                                      <  1:113832/251‑1 (MQ=255)
                                                                                                             ggATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGaacggcaaaggc                                                                                                                                      >  2:525923/1‑251 (MQ=255)
                                                                                                                 gcgcTTCGCTTATCCGGCCTACACGCGATGCCAGGGCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCt                                                                                                                                  <  1:472581/251‑1 (MQ=255)
                                                                                                                  cgcttcgcttATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGCCATTAAACTGATTAAGACTTCGTCGATTTCTCAACGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCACTGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACCATCACGCGGATCTGCGGCATGTTGCCCAGACGCCGCCCTTTGGTGAACGGCAACGGCTGCtt                                                                                                                                 >  2:81344/1‑251 (MQ=255)
                                                                                                                  cgcttcgcttATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCtt                                                                                                                                 >  2:113439/1‑251 (MQ=255)
                                                                                                                           tacccGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGc                                                                                                                          <  2:131469/246‑1 (MQ=255)
                                                                                                                               cGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGaacggcaac                                                                                                                                         <  2:508859/230‑1 (MQ=255)
                                                                                                                               cGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGaacggcaac                                                                                                                                         >  1:508859/1‑230 (MQ=255)
                                                                                                                               cGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGa                                                                                                                    <  2:440258/251‑1 (MQ=255)
                                                                                                                                ggCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTACAGCTCTTCATCCACAATCACGCGGATCTGCGGCATTTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTTTTCAgactctgcc                                                                                                                   >  2:388295/1‑243 (MQ=255)
                                                                                                                                       acGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTcgc                                                                                                            <  1:517361/251‑1 (MQ=255)
                                                                                                                                       acGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTcgc                                                                                                            <  1:570619/251‑1 (MQ=255)
                                                                                                                                        cgcgATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCGGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTcgcg                                                                                                           <  1:81344/251‑1 (MQ=255)
                                                                                                                                        cgcgAGGCCAGGTCAGAGAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTcgct                                                                                                           <  2:332275/251‑2 (MQ=255)
                                                                                                                                            aTGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCAACTGACGCGTCGCGGTAt                                                                                                       >  2:194819/1‑251 (MQ=255)
                                                                                                                                              gCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATc                                                                                                      <  1:4258/250‑1 (MQ=255)
                                                                                                                                               ccAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGGGTCTTCTTACCGGTCCACCCGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGGGAACGGCAACGGCGGCTTATTCAACCTCTGACGCGTCGCGGTATccc                                                                                                    >  2:418194/1‑250 (MQ=255)
                                                                                                                                               ccAGGTCAGATAAGCGAAACGCAGCCGACTGTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTTTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCg                                                                                                    <  1:606311/251‑1 (MQ=255)
                                                                                                                                                cAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGc                                                                                                   <  1:150791/251‑1 (MQ=255)
                                                                                                                                                    caataTAAGCGAACCGCACCCGACTTTACACTGATTACGCCTCCGTCTATTTCTCAAAGCGGCGCAGCAACTGTGTCCCACGCTCAACGGCGGTATCCAATGCCGGCTGTGGGGTCTTCTTCCCGGTCCACACGCTCTCCAGCTCTTCATCCCCAATCACGCGGATCTGCGGCATGTGGCGCAGCCGCACCCCTTTGGGGAACGGCAACGGCGGCTTATTCAGCACCTGACGCGTCGCGGTGTCCGCCCCt                                                                                               <  2:303095/247‑1 (MQ=255)
                                                                                                                                                      aGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTgg                                                                                             <  1:369602/251‑1 (MQ=255)
                                                                                                                                                      aGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTgg                                                                                             <  1:425536/251‑1 (MQ=255)
                                                                                                                                                       gATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTggg                                                                                            >  2:337133/1‑251 (MQ=255)
                                                                                                                                                        aTAAGCGAAACGCATCCGACATTACACTGACTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGt                                                                                           <  1:245859/251‑1 (MQ=255)
                                                                                                                                                          aaGCGAAACGCATCCGACATTTCACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGtt                                                                                          <  2:525472/250‑1 (MQ=255)
                                                                                                                                                               aaaCGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTtctc                                                                                    <  1:247091/251‑1 (MQ=255)
                                                                                                                                                                        cGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGcc                                                                           <  1:247692/251‑1 (MQ=255)
                                                                                                                                                                           cATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCACTGCCTGCTGTGGGGTCTTCTTACCGGTCCACACGCTCTCCCGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGAGAACGGCAACGGCGGCTTATTCCGCCTCTGACGCGTCGCGGTAACCGCCCCTGGGGTTTTCTCATAAAAGCCCTg                                                                        >  2:629523/1‑251 (MQ=255)
                                                                                                                                                                            aTTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTTATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTTGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCAGAAAAGCCCTGc                                                                       <  1:546228/251‑1 (MQ=255)
                                                                                                                                                                                 acTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACg                                                                  >  1:596610/1‑251 (MQ=255)
                                                                                                                                                                                     aTTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAAGCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTGATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGggtc                                                              <  1:356833/251‑1 (MQ=255)
                                                                                                                                                                                         aGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGACTTCTTACCGGTCCACACGCTCTACAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCATTTGGTGCACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGAATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGgggcaggt                                                          >  1:566174/1‑251 (MQ=255)
                                                                                                                                                                                            cTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCAt                                                       >  1:584198/1‑251 (MQ=255)
                                                                                                                                                                                                gtcgATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCttt                                                                                                                                                      >  2:407964/1‑152 (MQ=255)
                                                                                                                                                                                                gtcgATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCCCGGGTCAGGTCATACgc                                                   >  2:107160/1‑251 (MQ=255)
                                                                                                                                                                                                         ccataGCGGCGCAGCAACTGATTGCCTCGCTCATCGGCGGTTTCCAATTCCAGCTGTGGTGTCGTCTTACCGTTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGAATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAGCGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGt                                          <  2:522443/247‑1 (MQ=255)
                                                                                                                                                                                                          cAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCGGCGGCATGTTGCCCAGACGCAGCCCTTTGGGGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTg                                         <  1:348391/251‑1 (MQ=255)
                                                                                                                                                                                                           aaaGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTg                                         >  1:413646/1‑250 (MQ=255)
                                                                                                                                                                                                           aaaGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGa                                        <  1:327527/251‑1 (MQ=255)
                                                                                                                                                                                                           aaaGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGa                                        <  1:629523/251‑1 (MQ=255)
                                                                                                                                                                                                              gcggcgCATCAACTGATTGCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTg                                     <  2:230602/251‑1 (MQ=255)
                                                                                                                                                                                                                            gATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTGCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGtt                       <  2:618872/251‑1 (MQ=255)
                                                                                                                                                                                                                             aTTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGttt                      <  1:280823/251‑1 (MQ=255)
                                                                                                                                                                                                                             aTTTCAACGCTCACCGGCGGTATCCAATGCCTGCTGCGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGGATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCCCGGGTCAGGGCATACGCTGCTTTTGTGATTGGCCGCTAACCGGttt                      >  2:459804/1‑251 (MQ=255)
                                                                                                                                                                                                                              tttCCACGCTCAACGGCGGGATCCCATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCCGACGCAGCCCTTTGGTGGACGGCAACGGCGGCTTATTCAGCATCTGGCGCGTCGCGGGATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCCCGGGTCAGGGCATCAGCTGCTTTCGGGATTGGCAGGTCACCGGtttt                     >  2:95685/1‑251 (MQ=255)
                                                                                                                                                                                                                               ttCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTc                    >  1:68589/1‑251 (MQ=255)
                                                                                                                                                                                                                               ttCCAAGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTc                    >  1:504138/1‑251 (MQ=255)
                                                                                                                                                                                                                                tCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGGTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTc                    <  2:279005/250‑1 (MQ=255)
                                                                                                                                                                                                                                 cctcgcgcAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCTCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGATTTTCTg                  <  2:263148/244‑1 (MQ=255)
                                                                                                                                                                                                                                     gCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGc              >  1:592460/1‑251 (MQ=255)
                                                                                                                                                                                                                                         aaCGGCGTTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCTGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACt          <  1:49491/251‑1 (MQ=255)
                                                                                                                                                                                                                                               ggTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGc      >  2:12973/1‑249 (MQ=255)
                                                                                                                                                                                                                                                 tATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACACTCTCCAGCTCTTCATCAACAATCACGAGGGTCTGCGGCATGTTGCCCAGAAGCAGCCCTTTGCTGAACGGCAACGGCGACTTATTAAGGATCTGACGCGTCGCTGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCGCACAGGTCAGGTCATACGCTGCTTTGGTGCTTGGCAGATATCCGGGTTTCTGATGCCACTCGGCAGCg  >  2:71753/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
AAAAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAGTGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCG  >  NC_000913/3590827‑3591318

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