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.5% 179.0 / 0.2 66L418L (CTG→TTG) ugpBsn‑glycerol 3‑phosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (28/35);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

AAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAGTGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGG  >  NC_000913/3590829‑3591313
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
aaGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCct                                                                                                                                                                                                                                             <  1:217697/251‑1 (MQ=255)
 aGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCctg                                                                                                                                                                                                                                            <  1:439093/251‑1 (MQ=255)
          tgatgatgaGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGgtgt                                                                                                                                                                                                                                   >  1:639669/1‑251 (MQ=255)
                      gcggcgCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGt                                                                                                                                                                                                                       >  1:550593/1‑251 (MQ=255)
                        tgcgcGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGGGTCTTCTTACCGGTcc                                                                                                                                                                                                                     >  2:164845/2‑251 (MQ=255)
                          cgcgACAAGCAGATAAGGCAGCCAGCCCGCGCGGAACACCGGACGGGATGCAGAAATTACGGTTAATTCCTGAACAGTGCAGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGATCAACGTCAGTATCCAATGGCTGCTGTGGTGTCTTACTCCCGGTCcac                                                                                                                                                                                                                   >  2:91482/1‑251 (MQ=255)
                                       tAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGct                                                                                                                                                                                                     >  1:289638/1‑251 (MQ=255)
                                        aaGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCCCACGCTCTCC‑AGctc                                                                                                                                                                                                    >  1:271334/1‑251 (MQ=255)
                                         aGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑CCGCATCCGACCTTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGCCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGACCCCCCGCTCTCC‑AGctct                                                                                                                                                                                                   >  2:19168/1‑251 (MQ=255)
                                           gCAGTCAGCGTGAGCGGAACACCGGACGGGTTGATGACATTACGGTTAAGTCCTGTACTTTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCTACATTACACTGATTAAGACTTCTTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGGCCACACGCTCTCC‑AGCTCTTc                                                                                                                                                                                                 <  2:577589/251‑1 (MQ=255)
                                           gCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTc                                                                                                                                                                                                 <  1:163390/251‑1 (MQ=255)
                                             agccagcGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCAt                                                                                                                                                                                               >  1:540327/1‑251 (MQ=255)
                                             agccagcGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGAGTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCAt                                                                                                                                                                                               <  1:617892/251‑1 (MQ=255)
                                             agccagcGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAAATCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTAc                                                                                                                                                                                                                           >  1:320862/1‑224 (MQ=255)
                                             agccagcGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAAATCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTAc                                                                                                                                                                                                                           <  2:320862/224‑1 (MQ=255)
                                              gccagcGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCGGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATc                                                                                                                                                                                              <  1:454056/251‑1 (MQ=255)
                                               ccagcGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATcc                                                                                                                                                                                             >  1:163245/1‑251 (MQ=255)
                                               ccagcGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGGCCACACGCTCTCC‑AGCTCTTCATcc                                                                                                                                                                                             <  1:477640/251‑1 (MQ=255)
                                                cagcGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCACTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATcc                                                                                                                                                                                             >  1:404002/1‑250 (MQ=255)
                                                   cgcgAGCGGAACACCGGAGGGGATGATGACATTACGGTTAATTCCGGAACTGGGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCGCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACaa                                                                                                                                                                                         <  2:38450/251‑1 (MQ=255)
                                                     cgAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATc                                                                                                                                                                                       >  1:588614/1‑251 (MQ=255)
                                                      gAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCa                                                                                                                                                                                      <  2:639669/251‑1 (MQ=255)
                                                          ggAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCAcgcg                                                                                                                                                                                  >  1:553382/1‑251 (MQ=255)
                                                           gaccacCGGACTGGTTGATGAAATTACGGTTAAGTCCTGTACTGGGCCGGACTCGCTTCGCTTATCCGGCTTACACGCGATGCCAGGGCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCgg                                                                                                                                                                                 <  2:271334/248‑1 (MQ=255)
                                                                   gACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCGTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCgg                                                                                                                                                                         <  1:393228/251‑1 (MQ=255)
                                                                           gatgaCATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCGCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGc                                                                                                                                                                 <  1:214812/251‑1 (MQ=255)
                                                                           gatgaCATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATCACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAGTGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGGTCTGCGGCCTGTTGc                                                                                                                                                                 >  2:312664/1‑251 (MQ=255)
                                                                              gaCATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCa                                                                                                                                                              <  1:407956/251‑1 (MQ=255)
                                                                                 aTTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATCACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAGTGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGAc                                                                                                                                                           <  1:312664/251‑1 (MQ=255)
                                                                                              cctgaacgtggcctgATGCGCTTCTCTTATCCGGCCTACACGCGATGCCAGTGAA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATGAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCGTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTgg                                                                                                                                               <  2:649583/236‑1 (MQ=255)
                                                                                                     tgtgCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAGTGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGaacggc                                                                                                                                       >  1:103981/1‑251 (MQ=255)
                                                                                                           gactgcgctccctttATCCGGCCTACACGCGATGCCAGTGCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGTTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCATACGCAGCCCTTTGGTGaacggcaacggc                                                                                                                                 <  2:468750/238‑1 (MQ=255)
                                                                                                                cgcttcgcttATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCtt                                                                                                                            <  1:651505/251‑1 (MQ=255)
                                                                                                                 gcttcctttatccggccgaacctagATGCCGGGTCG‑GATAAGCGAA‑ACGCATCCGAAATTACACTGGTTAAGGCTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTa                                                                                                                           <  2:172558/227‑1 (MQ=255)
                                                                                                                  cttcgcttATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTAt                                                                                                                          >  1:278838/1‑251 (MQ=255)
                                                                                                                          agccGGCCTCCACGGGATGCAAGG‑CACGATAAGCGAA‑CCGCATCCTACATTACATTGATTAAGCCTTCGTCGATTTCTCAAATCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTACCCAATGCCTGCTGTGGTGACTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGAAACGGCGGCTTATTCAGCATc                                                                                                                  <  2:641473/249‑1 (MQ=255)
                                                                                                                          aTCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGctct                                                                                                                                                                                                   <  1:639702/170‑1 (MQ=255)
                                                                                                                          aTCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCCCTCC‑AGctct                                                                                                                                                                                                   >  2:639702/1‑170 (MQ=255)
                                                                                                                             cGGCCTACACGAGATGCCAGGTCA‑GATAAGCGAA‑ACGCATACGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGGACCCAATGACCACTGTGGTGTCTTCTTCCCTGTCCACACCCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGAGGCCTGATGCCCCGACGCAGCCCCTTGGTGAACGGCAACGACGGCTTaatatccacctgg                                                                                                               >  2:585795/1‑239 (MQ=255)
                                                                                                                                  tatgcgcaaTGCCAGGTCA‑GAAAAGCGAACAC‑AATCCAACATTACACTGATTAAGACTTCGTCGATTTCTCAAATCTGCGCAGCACCTGAATTCCCCGCTCAACGGCGTTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTATTCAGCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGAcgc                                                                                                            <  2:233989/241‑1 (MQ=255)
                                                                                                                                            gCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCATCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATcc                                                                                                <  1:621018/251‑1 (MQ=255)
                                                                                                                                                 gTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCCTCTGACGCCTCCCGGtggatcgga                                                                                            >  2:543457/1‑242 (MQ=255)
                                                                                                                                                 gTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGt                                                                                                    <  1:543457/242‑1 (MQ=255)
                                                                                                                                                   cA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCt                                                                                          <  1:157030/250‑1 (MQ=255)
                                                                                                                                                    gcagttaatcgaacagcaaccGACATTACACGGATTAATCCTTCGTCGATTTCTCAAAGCGGCGAAGCAACTGATTACCCCGCCCAACGGCGTTAGCCAATGCCTGCGTTGGTGTCTTCTTCCCGGTCCACACGCTCTCC‑AGCTCTTCATTCACAAGCCCGCGGATCTGCGGCATGTTGCCCAGCCGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCt                                                                                          <  2:271226/232‑1 (MQ=255)
                                                                                                                                                    a‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTTATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTgg                                                                                        <  1:34045/251‑1 (MQ=255)
                                                                                                                                                      agataagcgacactatcccGACATTACATAGATTAAGACTGCTGCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTgg                                                                                        <  2:473453/234‑1 (MQ=255)
                                                                                                                                                                 ccGCATCCGCCATTGCACTGATGAAGACTTCGTCGATTTCACAAAGCGGCGCAGCAACTGATTCCCACGCTCAACGGCGTTAGCCAAGGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCCCAATCCCGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCGCAt                                                                             <  2:80196/250‑1 (MQ=255)
                                                                                                                                                                   ccATCCGACATTACACTGATTAAGACTTCGGCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGTTCCACACGCTCTCC‑AGCTCTTCATCCACAATCTCGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATaa                                                                           <  1:679543/250‑1 (MQ=255)
                                                                                                                                                                    cATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGGATCCGCCCCTGGGTTTTTCTCATaa                                                                           >  1:501649/1‑250 (MQ=255)
                                                                                                                                                                          aCATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCt                                                                    >  1:327858/1‑251 (MQ=255)
                                                                                                                                                                                   tGATTAAGACTTCGTCGATTTCTCAAAGAGGCGCAGCAACTGATTTCAACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGGTTTTCTCATAAAAGCCCTGCTCACggg                                                           >  1:604485/1‑251 (MQ=255)
                                                                                                                                                                                               cgtcgATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACg                                               >  1:163242/1‑251 (MQ=255)
                                                                                                                                                                                                  tgttttCTCAAAGTGGCGCTGCAACTTTGTTCCACGCTAAAGGCCGGTATCCAATGCCTGTTGGGGTTCCTTCTTACCGGTCCACCCGCTCTCC‑AGCTCTTCAGCCACAACCACGCGGATCTGCGGCATGTTGCCCAGACTCAGCCCTGTGGTTAACGGCAACGGCGGCTTATTCAGCACCTTACTCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACgct                                             <  2:654885/247‑1 (MQ=255)
                                                                                                                                                                                                    aTTTCTCAAAGCGGAGCAGCAACTGATTTCCACGCTCAACGGCGGTCTCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACACTCACGCGGATCTGCGGCATATTGCCCAAACACAACCCTTTGGTGAACGGCAAAGGCGGCTTATTCAGCAACTGACGCGCCGCGGGCTCCGCCCCTGGGGTTTTTTACTACAAACCCTGCCCACGGGTAGGGGCATACgctcct                                          >  2:630968/1‑248 (MQ=255)
                                                                                                                                                                                                     tttCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCCCACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCCGCATCTGACGCGTcgcg                                                                                                      >  2:183689/1‑190 (MQ=255)
                                                                                                                                                                                                     tttCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTcgcg                                                                                                      <  1:183689/190‑1 (MQ=255)
                                                                                                                                                                                                          cAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTgg                                      <  1:80452/249‑1 (MQ=255)
                                                                                                                                                                                                            aaGCGGCGCAGCAACTGATTTCCACGCCCACCGCCGGTATCCAATGCCTGCTGTGGTTCCTTTTTACCGGCCCAC‑CCCCCTCCCACCTCTTCATCCACAATCCCCCGGATCTGCGCCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGGGCCTTATTCAGCATTTGACGCGCCGCGGTATCCCCCCCTGGGTTTTTTTCATAAAAACCCTCCTCACGGGTCGGGCCATCCGCTGCTTTggggtt                                  >  1:4815/1‑247 (MQ=255)
                                                                                                                                                                                                                  cgcagataCTGATTTCCCCGCTCAACGGCGGTAACCAATGCCTGCTGGGGTGTCTTCTGTCCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAg                            <  2:550593/244‑1 (MQ=255)
                                                                                                                                                                                                                  cgcAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCCTAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGc                              >  1:145169/1‑249 (MQ=255)
                                                                                                                                                                                                                   gcagcaACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTTTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGc                              >  1:404667/1‑248 (MQ=255)
                                                                                                                                                                                                                               ttCCACCCGCAACGGCGGTATCCAATGCCCGTTGTGGTGTCCTCTTACCGGTCCACACGCGCGCA‑AGGCCTTCAACCCAAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGtttt                <  1:495582/250‑1 (MQ=255)
                                                                                                                                                                                                                                      cTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGcc        >  1:94040/1‑251 (MQ=255)
                                                                                                                                                                                                                                            ggcggTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTTTCATACAAGCCCTGCTCACGGGTCCGGACATACGCCTCTTTGTTGTTTTGCAGATTACCGGTTTTCTTATgaccaccgg  >  2:297005/1‑243 (MQ=255)
                                                                                                                                                                                                                                            ggcggTAGCCAATGCCTGCTGGGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCgg  <  1:203254/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
AAGATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCA‑GATAAGCGAA‑ACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAGTGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCC‑AGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGG  >  NC_000913/3590829‑3591313

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