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

GTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAGTGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGC  >  NC_000913/3590838‑3591314
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 
gtgaggaggaGCAGCGGGGCGACCAGCAGATAAGGCAGCCAGCCCGAGAGGACCACCGGGCGGGATTCTGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGGTTCGCATACCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCACCCGACATTACACTGATGAAGTCTGCGTCGATTTCTCAAAGCGGCGCAGCATCCGAGTTCCACGCTCAACGGCGGTACCCAATGCCTGCTgtggtg                                                                                                                                                                                                                                     <  2:69405/243‑1 (MQ=255)
     gatgaGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCCACGGAGGTATCCAATGCCTGCTGTGGTGtctt                                                                                                                                                                                                                                 <  1:200186/250‑1 (MQ=255)
      atgaGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGGGGGGtcttct                                                                                                                                                                                                                               >  2:416908/1‑251 (MQ=255)
           cTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTACACGCTCAACGGCGGTATACAATGCCTGCTGTGGTGTCTTCTTACCg                                                                                                                                                                                                                          >  1:176172/1‑251 (MQ=255)
                               aaGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGct                                                                                                                                                                                                      <  1:36441/250‑1 (MQ=255)
                                    agccagcGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCAt                                                                                                                                                                                                >  1:28861/1‑251 (MQ=255)
                                      ccagcGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGTATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTGCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATcc                                                                                                                                                                                              <  1:501731/251‑1 (MQ=255)
                                                 ggAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTTCAGCTCTTCATCCACCATCAcgcg                                                                                                                                                                                   >  1:86246/1‑251 (MQ=255)
                                                 ggAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCCACCGCTCTCCAGCTCTTCCTCCACAATCAcgcg                                                                                                                                                                                   >  2:172296/1‑251 (MQ=255)
                                                      acCGGACGGGGTGATGACATTACGGATAATTCCTGAACTGTGCCGGATGCGCTGCGCTTATCCGGCCTACACGCGATGCCAGGACAGATAAGAGAAACGCATCCGACATTAAACTGATTAAGACTTCGACGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTAAAAGACGGTATCCCATGACTGCTGGGGTGACTTCTTACCGGTCCACAAGTCCTTCTGCTTTTCATCCAAAAACACGCGGAttt                                                                                                                                                                              >  2:345771/1‑249 (MQ=255)
                                                         ggACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCg                                                                                                                                                                           >  1:118346/1‑251 (MQ=255)
                                                            cggtattatgaCATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCGGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCAGCCACAATCACGCGGATCTGCGGCa                                                                                                                                                                        <  2:473760/244‑1 (MQ=255)
                                                                         ttACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACg                                                                                                                                                           >  1:43575/1‑251 (MQ=255)
                                                                              gTTAATTCCTGAACGCTGCCAGAGGCGCTTCGCTTAGCCGGCCTACATGCGATTCCAGGTCATACAAGCGAAACGCATCCGACATTACCCCGATTAAGACTGCTTCGATTTCTCAAAGCGGCGCAGCAACGGATTTCCACGCCCAACGCCGCTAACCAATGCCTGCTGTGGTGTCTTCTTACCGGGCCACTCGCTCCCCAGCTCTTCACCCACAATCACGCGGATCTGCGGCATGTTGCCCATACGCAGcc                                                                                                                                                      <  2:364863/251‑1 (MQ=255)
                                                                                 aaTTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCt                                                                                                                                                    <  2:449550/250‑1 (MQ=255)
                                                                                  aTTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGGTTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGTCTGCTGTGGTGTCTTCTTACCGGTACAAACGCTCTCCACATCTGCATCCACAATCACGAGGATCTGCGGCATGTTGCCCAGACCCAGCCCttt                                                                                                                                                  >  1:387518/1‑251 (MQ=255)
                                                                                             gtgCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGaacgg                                                                                                                                         <  1:479889/249‑1 (MQ=255)
                                                                                                  ggctGCGCTTCGCTTATCCGGCCTACCCGCGATGCCAGGTCAGATAAGCGAAACGCATCAGACATTACACTTATTAAGACTTCGTCGATTGCTCAAAGCGGCGAAGCATCTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGaacggcaacggc                                                                                                                                  <  2:170636/248‑1 (MQ=255)
                                                                                                  ggATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGaacggcaacggc                                                                                                                                  <  1:416908/251‑1 (MQ=255)
                                                                                                      gggcttcgcttACCCGGCCTACACGCGATGCCAGGTCAGATAAGGGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATGTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACggcggc                                                                                                                               <  2:90683/248‑1 (MQ=255)
                                                                                                                tATCCGGCCTACCCGGGAAGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGCCTGCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAGTGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGc                                                                                                                      <  2:158266/249‑1 (MQ=255)
                                                                                                                 aTCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGCTCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATc                                                                                                                   <  1:40676/251‑1 (MQ=255)
                                                                                                                     ggCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTACAGCTCTTCATCCACAATCAAGCGGATCTGCGGCATGTTGCCCAGAAGCAGCACTTTGGTGAACGGCAACGGCGGATTTTTCAGCATCTGAc                                                                                                               >  1:319107/1‑251 (MQ=255)
                                                                                                                           cacGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTcg                                                                                                         <  1:382905/251‑1 (MQ=255)
                                                                                                                               cgATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGt                                                                                                     <  2:7789/251‑1 (MQ=255)
                                                                                                                                  ggCCAGGTCAGATAAGCGAAACGCATCCGACATGACACTGAGTAAGCCTTCGGCGATTTCTCAAAGCGGCGCAGCACCGGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTGACCGGTCCCCACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATc                                                                                                  <  2:67798/250‑1 (MQ=255)
                                                                                                                                    ccAGGTCAGATAAGCGAAACGCATCCGACATTAAACTGATTAAGACTTCGTCGATTTCTCAAAGCAGCGCAGCAACTGATTTCCACGCTCAACGGCGGTAACCCATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCCGCCCTTCATCCAAAATCACGCGGATCTGCGGCATGTTGACCCGACGCAGCCCTTTGGGGGGCGGCAACGGCGCATTATTCaacagctgcactctcgcggtattcc                                                                                                >  2:183239/1‑227 (MQ=255)
                                                                                                                                     cAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCACAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAAACGCAGCCCTTTGGAGAACGGCCAAGGCGGCTTATTCAGCCTCTGCCGCGTCCCgggataccc                                                                                               >  2:489464/1‑244 (MQ=255)
                                                                                                                                           aGATAAGCGAAACGCACCGTAAATTACCCTGAGTAAACCTCCGCCGATTGCTCAAAGCGGCGCTGCAACGGTTTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGGGGGGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTGCAGCCACAATCGCGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCGTGGGTTAACGGCAGCGGCGGCTTATTCAGCATCTGACGCGTCGCGGTAGCCGCGCCTgg                                                                                         <  2:319107/251‑1 (MQ=255)
                                                                                                                                            gATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGGATCCGCCCCTggg                                                                                        >  1:168472/1‑251 (MQ=255)
                                                                                                                                            aaTAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTggg                                                                                        <  1:257151/250‑1 (MQ=255)
                                                                                                                                                       cGCATCCGACAATACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGAAACTGATTTCAACGCTCAACGGCGGAATCCAATGCCTGCTGGGGGGTCTTCTTAACGGTCCACACGCTCTCCAGCTCTTCATCCACCATCACGCGGATCTGCGGCCTGGTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGATTATTCACGACCTGACGCGTCCCGGTGTCCCACCCTGGGgtttttttct                                                                              >  2:494631/1‑245 (MQ=255)
                                                                                                                                                                         tGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGACCTGCTCACggg                                                            >  1:42770/1‑251 (MQ=255)
                                                                                                                                                                           aTTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGggtc                                                          <  1:79813/251‑1 (MQ=255)
                                                                                                                                                                            gtAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATACAAAGCCTGCTGGGGTGTGTTCTTACAGGTGCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCTGCCCTTTGGTGGACGGAAACGGCGGCTTATTCAGCCACTGACGGGTCGCGGTAGTCTCCCCTGGGTTTTTCTCATAAAAgtactggtaacgtggcg                                                         >  2:99976/2‑235 (MQ=255)
                                                                                                                                                                                gccTTCGTCGGTTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCCCAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCTCGGTATCCGCCCCTGGGATTTTCTCATTAAACCCCTGCTCACGgctcatgtc                                                     >  2:495945/3‑251 (MQ=255)
                                                                                                                                                                                         gacttcacaaggcggcgCAGCAACGAATTTCCACGCTCACCGGCGGTATCCCATGCCTGCTGTGGGCTCTACCTACCGGTCCACACGCTCACCAGCTCTTTATCAACGATCACGCGGATCTGCGGCATGTTGCCCAGAGGCCGGCCATTGTGGAACGGCAACGGCGGGTTATTAAGCATCTGTCGCGGCGCTGTGTCCGTGTCTGGGTTTTTCGCATAAAAGCCCTGCGCAAGGGTCAGGTCATACgctg                                             <  2:278328/239‑1 (MQ=255)
                                                                                                                                                                                                cAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACCCGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTgg                                       <  1:221583/249‑1 (MQ=255)
                                                                                                                                                                                                 aaaGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTgggg                                     >  1:244279/1‑248 (MQ=255)
                                                                                                                                                                                                 aaaGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGGTTTTCTCATAAAAGCCCTGCTCACGGGTCCGGTCATACGCTGCTTTggggg                                    >  1:106901/1‑248 (MQ=255)
                                                                                                                                                                                                     cagcgtcgcaACTGATTTCCACGCTCAACGGTGGTATCCAATGACTGCTGTGGGGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTGCACCCACAATCACGCGGATCTGCGGCAGGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCTGCGCGGTATCCGCCCCTGGGGTTTTCTCATAAAAGCCCTGCACACGGGTCAGGTCATACGCTGCTTTGGTGATTgg                                <  2:86246/244‑1 (MQ=255)
                                                                                                                                                                                                              aaCTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCCGGTCATACGCTGCTTTGGTGATTGGCAGATAAcc                       >  1:231545/1‑251 (MQ=255)
                                                                                                                                                                                                                        cACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACCCGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGa             <  1:35584/251‑1 (MQ=255)
                                                                                                                                                                                                                         ccGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAAGCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGAt            <  1:281679/250‑1 (MQ=255)
                                                                                                                                                                                                                                   gcggTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGc  >  1:87604/1‑251 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 
GTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAGTGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGC  >  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: