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,0740GA97.8% 129.3 / ‑4.2 45L418L (CTG→TTG) ugpBsn‑glycerol 3‑phosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (21/23);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-01

GTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCG‑ACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAGTGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGC  >  NC_000913/3590838‑3591314
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 
gtgaggaggaGCAGCGGGGCGACCAGCAGATAAGGCAGCCAGCCCGAGAGGACCACCGGGCGGGATTCTGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGGTTCGCATACCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCACCCG‑ACATTACACTGATGAAGTCTGCGTCGATTTCTCAAAGCGGCGCAGCATCCGAGTTCCACGCTCAACGGCGGTACCCAATGCCTGCTgtggtg                                                                                                                                                                                                                                     <  2:69405/243‑1 (MQ=255)
     gatgaGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCG‑ACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCCACGGAGGTATCCAATGCCTGCTGTGGTGtctt                                                                                                                                                                                                                                 <  1:200186/250‑1 (MQ=255)
      atgaGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCG‑ACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGGGGGGtcttct                                                                                                                                                                                                                               >  2:416908/1‑251 (MQ=255)
           cTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCG‑ACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTACACGCTCAACGGCGGTATACAATGCCTGCTGTGGTGTCTTCTTACCg                                                                                                                                                                                                                          >  1:176172/1‑251 (MQ=255)
                               aaGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCG‑ACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGct                                                                                                                                                                                                      <  1:36441/250‑1 (MQ=255)
                                    agccagcGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCG‑ACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCAt                                                                                                                                                                                                >  1:28861/1‑251 (MQ=255)
                                      ccagcGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGTATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCG‑ACATTACACTGATTAAGACTTCGTCGATTGCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCATcc                                                                                                                                                                                              <  1:501731/251‑1 (MQ=255)
                                                 ggAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCG‑ACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTTCAGCTCTTCATCCACCATCAcgcg                                                                                                                                                                                   >  1:86246/1‑251 (MQ=255)
                                                 ggAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCG‑ACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCC‑CACCGCTCTCCAGCTCTTCCTCCACAATCAcgcg                                                                                                                                                                                   >  2:172296/1‑251 (MQ=255)
                                                      acCGGACGGGGTGATGACATTACGGATAATTCCTGAACTGTGCCGGATGCGCTGCGCTTATCCGGCCTACACGCGATGCCAGGACAGATAAGAGAAACGCATCCG‑ACATTAAACTGATTAAGACTTCGACGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTAAAAGACGGTATCCCATGACTGCTGGGGTGACTTCTTACCGGTCCACAA‑GTCCTTCTGCTTTTCATCCAAAAACACGCGGAttt                                                                                                                                                                              >  2:345771/1‑249 (MQ=255)
                                                         ggACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCG‑ACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCg                                                                                                                                                                           >  1:118346/1‑251 (MQ=255)
                                                            cggtattatgaCATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCG‑ACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCGGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCAGCCACAATCACGCGGATCTGCGGCa                                                                                                                                                                        <  2:473760/244‑1 (MQ=255)
                                                                         ttACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCG‑ACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACg                                                                                                                                                           >  1:43575/1‑251 (MQ=255)
                                                                              gTTAATTCCTGAACGCTGCCAGAGGCGCTTCGCTTAGCCGGCCTACATGCGATTCCAGGTCATACAAGCGAAACGCATCCG‑ACATTACCCCGATTAAGACTGCTTCGATTTCTCAAAGCGGCGCAGCAACGGATTTCCACGCCCAACGCCGCTAACCAATGCCTGCTGTGGTGTCTTCTTACCGGGCCACTC‑GCTCCCCAGCTCTTCACCCACAATCACGCGGATCTGCGGCATGTTGCCCATACGCAGcc                                                                                                                                                      <  2:364863/251‑1 (MQ=255)
                                                                                 aaTTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCG‑ACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCt                                                                                                                                                    <  2:449550/250‑1 (MQ=255)
                                                                                  aTTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCG‑ACATTACACTGATTAAGACTTCGTCGGTTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGTCTGCTGTGGTGTCTTCTTACCGGTACAAAC‑GCTCTCCACATCTGCATCCACAATCACGAGGATCTGCGGCATGTTGCCCAGACCCAGCCCttt                                                                                                                                                  >  1:387518/1‑251 (MQ=255)
                                                                                             gtgCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCG‑ACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGaacgg                                                                                                                                         <  1:479889/249‑1 (MQ=255)
                                                                                                  ggctGCGCTTCGCTTATCCGGCCTACCCGCGATGCCAGGTCAGATAAGCGAAACGCATCAG‑ACATTACACTTATTAAGACTTCGTCGATTGCTCAAAGCGGCGAAGCATCTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGaacggcaacggc                                                                                                                                  <  2:170636/248‑1 (MQ=255)
                                                                                                  ggATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCG‑ACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGaacggcaacggc                                                                                                                                  <  1:416908/251‑1 (MQ=255)
                                                                                                      gggcttcgcttACCCGGCCTACACGCGATGCCAGGTCAGATAAGGGAAACGCATCCG‑ACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATGTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACggcggc                                                                                                                               <  2:90683/248‑1 (MQ=255)
                                                                                                                tATCCGGCCTACCCGGGAAGCCAGGTCAGATAAGCGAAACGCATCCG‑ACATTACACTGATTAAGCCTGCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAGTGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGc                                                                                                                      <  2:158266/249‑1 (MQ=255)
                                                                                                                 aTCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCG‑ACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCATCCACAATCACGCGGCTCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATc                                                                                                                   <  1:40676/251‑1 (MQ=255)
                                                                                                                     ggCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCG‑ACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTACAGCTCTTCATCCACAATCAAGCGGATCTGCGGCATGTTGCCCAGAAGCAGCACTTTGGTGAACGGCAACGGCGGATTTTTCAGCATCTGAc                                                                                                               >  1:319107/1‑251 (MQ=255)
                                                                                                                           cacGCGATGCCAGGTCAGATAAGCGAAACGCATCCG‑ACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTcg                                                                                                         <  1:382905/251‑1 (MQ=255)
                                                                                                                               cgATGCCAGGTCAGATAAGCGAAACGCATCCG‑ACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGt                                                                                                     <  2:7789/251‑1 (MQ=255)
                                                                                                                                  ggCCAGGTCAGATAAGCGAAACGCATCCG‑ACATGACACTGAGTAAGCCTTCGGCGATTTCTCAAAGCGGCGCAGCACCGGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTGACCGGTCCCCAC‑GCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATc                                                                                                  <  2:67798/250‑1 (MQ=255)
                                                                                                                                    ccAGGTCAGATAAGCGAAACGCATCCG‑ACATTAAACTGATTAAGACTTCGTCGATTTCTCAAAGCAGCGCAGCAACTGATTTCCACGCTCAACGGCGGTAACCCATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCCGCCCTTCATCCAAAATCACGCGGATCTGCGGCATGTTGACCCGACGCAGCCCTTTGGGGGGCGGCAACGGCGCATTATTCaacagctgcactctcgcggtattcc                                                                                                >  2:183239/1‑227 (MQ=255)
                                                                                                                                     cAGGTCAGATAAGCGAAACGCATCCG‑ACATTACACTGATTAAGACTTCGTCGATTTCTCACAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAAACGCAGCCCTTTGGAGAACGGCCAAGGCGGCTTATTCAGCCTCTGCCGCGTCCCgggataccc                                                                                               >  2:489464/1‑244 (MQ=255)
                                                                                                                                           aGATAAGCGAAACGCA‑CCGTAAATTACCCTGAGTAAACCTCCGCCGATTGCTCAAAGCGGCGCTGCAACGGTTTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGGGGGGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTGCAGCCACAATCGCGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCGTGGGTTAACGGCAGCGGCGGCTTATTCAGCATCTGACGCGTCGCGGTAGCCGCGCCTgg                                                                                         <  2:319107/251‑1 (MQ=255)
                                                                                                                                            gATAAGCGAAACGCATCCG‑ACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGGATCCGCCCCTggg                                                                                        >  1:168472/1‑251 (MQ=255)
                                                                                                                                            aaTAAGCGAAACGCATCCG‑ACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTggg                                                                                        <  1:257151/250‑1 (MQ=255)
                                                                                                                                                       cGCATCCG‑ACAATACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGAAACTGATTTCAACGCTCAACGGCGGAATCCAATGCCTGCTGGGGGGTCTTCTTAACGGTCCACAC‑GCTCTCCAGCTCTTCATCCACCATCACGCGGATCTGCGGCCTGGTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGATTATTCACGACCTGACGCGTCCCGGTGTCCCACCCTGGGgtttttttct                                                                              >  2:494631/1‑245 (MQ=255)
                                                                                                                                                                         tGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGACCTGCTCACggg                                                            >  1:42770/1‑251 (MQ=255)
                                                                                                                                                                           aTTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGggtc                                                          <  1:79813/251‑1 (MQ=255)
                                                                                                                                                                            gtAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATACAAAGCCTGCTGGGGTGTGTTCTTACAGGTGCACAC‑GCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCTGCCCTTTGGTGGACGGAAACGGCGGCTTATTCAGCCACTGACGGGTCGCGGTAGTCTCCCCTGGGTTTTTCTCATAAAAgtactggtaacgtggcg                                                         >  2:99976/2‑235 (MQ=255)
                                                                                                                                                                                gccTTCGTCGGTTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCATCCCCAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCTCGGTATCCGCCCCTGGGATTTTCTCATTAAACCCCTGCTCACGgctcatgtc                                                     >  2:495945/3‑251 (MQ=255)
                                                                                                                                                                                         gacttcacaaggcggcgCAGCAACGAATTTCCACGCTCACCGGCGGTATCCCATGCCTGCTGTGGGCTCTACCTACCGGTCCACAC‑GCTCACCAGCTCTTTATCAACGATCACGCGGATCTGCGGCATGTTGCCCAGAGGCCGGCCATTGTGGAACGGCAACGGCGGGTTATTAAGCATCTGTCGCGGCGCTGTGTCCGTGTCTGGGTTTTTCGCATAAAAGCCCTGCGCAAGGGTCAGGTCATACgctg                                             <  2:278328/239‑1 (MQ=255)
                                                                                                                                                                                                cAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACCC‑GCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTgg                                       <  1:221583/249‑1 (MQ=255)
                                                                                                                                                                                                 aaaGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTgggg                                     >  1:244279/1‑248 (MQ=255)
                                                                                                                                                                                                 aaaGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGGTTTTCTCATAAAAGCCCTGCTCACGGGTCCGGTCATACGCTGCTTTggggg                                    >  1:106901/1‑248 (MQ=255)
                                                                                                                                                                                                     cagcgtcgcaACTGATTTCCACGCTCAACGGTGGTATCCAATGACTGCTGTGGGGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTGCACCCACAATCACGCGGATCTGCGGCAGGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCTGCGCGGTATCCGCCCCTGGGGTTTTCTCATAAAAGCCCTGCACACGGGTCAGGTCATACGCTGCTTTGGTGATTgg                                <  2:86246/244‑1 (MQ=255)
                                                                                                                                                                                                              aaCTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCCGGTCATACGCTGCTTTGGTGATTGGCAGATAAcc                       >  1:231545/1‑251 (MQ=255)
                                                                                                                                                                                                                        cACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACCC‑GCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGa             <  1:35584/251‑1 (MQ=255)
                                                                                                                                                                                                                         ccGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCATCCACAAGCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGAt            <  1:281679/250‑1 (MQ=255)
                                                                                                                                                                                                                                   gcggTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGc  >  1:87604/1‑251 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 
GTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCG‑ACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAGTGCCTGCTGTGGTGTCTTCTTACCGGTCCACAC‑GCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGC  >  NC_000913/3590838‑3591314

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