Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,175,991 T→A K111* (AAG→TAG)  gatZ ← putative tagatose‑1,6‑bisphosphate aldolase 2 chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,175,9910TA95.1% 111.3 / ‑3.4 41K111* (AAG→TAG) gatZputative tagatose‑1,6‑bisphosphate aldolase 2 chaperone
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (17/22);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTTTACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCC  >  NC_000913/2175749‑2176235
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
ttCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCACTCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTCTGGGTTCCCCCGCGCAGGACATTGACGCATCAAGATGACTTTTACTGACGCCGGCACGAACATATTCCTATACCa                                                                                                                                                                                                                                                 >  2:460128/1‑248 (MQ=255)
ttCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGTAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCa                                                                                                                                                                                                                                                 <  1:460128/248‑1 (MQ=255)
      aggggTGATGGGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACgg                                                                                                                                                                                                                                        <  2:75035/249‑1 (MQ=255)
      aTGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGCTCCCCCGCGCAGGACAGTGACGCATCAAGATGAATTTTACTGTAGCCGGCCCGAACATTTTCCTTTCACAGCACGACgg                                                                                                                                                                                                                                        >  2:96713/1‑251 (MQ=255)
      aTGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACgg                                                                                                                                                                                                                                        >  1:414167/1‑251 (MQ=255)
                          gaATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTTGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGccgcagc                                                                                                                                                                                                                     <  1:37923/250‑3 (MQ=255)
                          gaATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGccgcatcc                                                                                                                                                                                                                    >  1:393539/1‑251 (MQ=255)
                                       atCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCACTCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACTTTGACGCATCAAGATGAATTTTACTGATGCCGGCACGAACATATTCCTATACCAGCTCGACGGTTTTTTACATCGCCGCATCCGCATTTTCTtgct                                                                                                                                                                                                       >  2:208812/2‑251 (MQ=255)
                                        tCACCGCCCGGACCCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTtgctg                                                                                                                                                                                                      >  1:117709/1‑251 (MQ=255)
                                        tCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTtgctg                                                                                                                                                                                                      >  1:10188/1‑251 (MQ=255)
                                                  gaaccggaacTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCACAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCg                                                                                                                                                                                                                                            >  1:410725/1‑203 (MQ=255)
                                                  gaaccggaacGTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCg                                                                                                                                                                                                                                            <  2:410725/203‑1 (MQ=255)
                                                   aaccggaacTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTg                                                                                                                                                                                           <  1:230699/251‑1 (MQ=255)
                                                             tCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAgat                                                                                                                                                                                  <  1:215814/250‑1 (MQ=255)
                                                                    cAATGACATAGCTCAGTTGCTCACGCTGGCAAGCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATcgccg                                                                                                                                                                          <  2:129095/251‑1 (MQ=255)
                                                                       tGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGccgccg                                                                                                                                                                       >  1:193528/1‑251 (MQ=255)
                                                                             agatCAGTTGCTCACACTGGCAATCTGACGCAACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCCGGACATTGACGCATCAAGATGAATTTTACTGCAGCCGGCACGAACATATTCCTATACAAGCTCGACGGATTTTTCTATCGCCCCATCCCCATTTTCTTGCTGCCAGCAGTGGGGCCCCAgctgaacaccaccgatact                                                                                                                                                                  >  2:313126/4‑236 (MQ=255)
                                                                             aGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGaat                                                                                                                                                                  <  1:265177/250‑1 (MQ=255)
                                                                                 cAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGTGCCCAGATGATCGCCGCCGAGAATAATg                                                                                                                                                              <  1:185177/250‑1 (MQ=255)
                                                                                      gCTCATGCTGGCAATCTGTCGCCACACTTTCCGCAGAAAAGCAAAGAACAGCAGCTCGTTCCGCAACCTTTTCTGGTGCTAACTGTATGGGATCCCCCGCGAAGGACATTGACGCAGCAAGATGAATTTTACTGAAGCCGGCACGAACATATCCCTTTACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTc                                                                                                                                                         <  2:166702/250‑1 (MQ=255)
                                                                                         cACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGt                                                                                                                                                     <  1:390264/251‑1 (MQ=255)
                                                                                                    tCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAAc                                                                                                                                          >  1:115690/1‑251 (MQ=255)
                                                                                                     ttGTCGCCATACTTTCCGCGGCAAAGCAAAGCACAGCATCTCGTTCCGCATCCGTTTCTGTTGTTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGAGGAATTTTACTGAAGCCTGCACGAACATATTCCGATACCAGCTCGACGGATTTTTCCATCGCCGCAGCCGCATTTTCTTGCTGCCAGCGGTTTGGCCCCATATGATCGCCGCCGATAATAATGCGTGCGCGTGCAAACCCAAc                                                                                                                                          <  2:177766/249‑1 (MQ=255)
                                                                                                         cGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTAt                                                                                                                                     <  1:299434/251‑1 (MQ=255)
                                                                                                                    ccGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTAGCGGCAATCGTa                                                                                                                          <  2:409349/251‑1 (MQ=255)
                                                                                                                       cagcaAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTaaa                                                                                                                        >  1:155738/1‑250 (MQ=255)
                                                                                                                          caaagcaaagCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAAc                                                                                                                     >  1:461053/1‑250 (MQ=255)
                                                                                                                              gcaaagccCAGCAGCTCGTTCCGCAACTGTTTCTGGTGCTAACGGTATTGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTGCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAAtt                                                                                                                <  2:320650/251‑1 (MQ=255)
                                                                                                                                              cgtgccGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGa                                                                                                           <  1:189437/236‑1 (MQ=255)
                                                                                                                                              cGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGAACTTGACGCATCAAGATGAATTTTACTGACGCCGGCCCGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGAATCCGCATTTTCTTGCTGCCAGCAATTTGGCCCCAGATGATCGCCGCCGAGCATAATGCGTTCCCGTGCAAACCCAACTTTATCGGCCATCGTAAAAACCAATTCGCGa                                                                                                           >  2:189437/1‑240 (MQ=255)
                                                                                                                                                 tcccccaCCGTTTCTGGTGCTAACGGTATGTGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTAGACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGgtg                                                                                              <  2:117709/244‑1 (MQ=255)
                                                                                                                                                   cGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTc                                                                                            >  1:88732/1‑250 (MQ=255)
                                                                                                                                                                 gTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCg                                                                             >  1:394459/1‑251 (MQ=255)
                                                                                                                                                                            gtGGGACCCCCCGCGCAGGACATTGACGCATCAAGAAGAATTTTACTGAAGCCGGCACGAACATATTCCAATCCCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGAAGTTTGGCCCCAGATGATCGCAGCCGAGAATAATGCGGTCGCGTGCAAACCCAACTTTACCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGCCATTCCGGGATAACCGCCAAattgatt                                                                  <  2:21286/250‑1 (MQ=255)
                                                                                                                                                                                gatcccccctccAGGACCTTGTCGCATCAAGAAGAATTTGACTGACGCCGGCACGAACATATTCCTATCCCAGCCCGATGGCTTTTTCCATCGCCGCATCCGTATTTTCTGGCTCCCAGCAGTTTGGCCCCAGATGATCTCCGCCGAGAATAATGCGTTCGCGTGCAAACCCACCTTTATCGGCAATCTTAAAAACAAATTAGCGAAAGTCGGCGGGTGTCATCCCGGTATAACCGCCAAATTGATTGAcc                                                              <  2:216768/240‑1 (MQ=255)
                                                                                                                                                                                                    gACGCATCAAGATGAATTTTATGTAAGCCGGCACGAACATATTCCTATACCAGCCGGACTGATTTTCCCATCGCCGCTTCCGCATCTTCTTGCTGCCAGCAGTTTGGCCCCAGATGACCCCCGCCGAGAATAATGCGTTCGCTGGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCTGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAAt                                          <  2:197301/251‑1 (MQ=255)
                                                                                                                                                                                                                          tGAAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCGGCACTTTGCGCGTGCTGTTg                    >  1:227022/1‑251 (MQ=255)
                                                                                                                                                                                                                           gAAGCCGGCACGAACATATTCCTATCCCAGCTCAACGGATTTTCCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGGTCGCGTGCAAACCCAACTTTACCGGCAATCGTAAAAACAAAGTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGtt                     <  2:122354/249‑1 (MQ=255)
                                                                                                                                                                                                                                      caaccATATCCCTATACCAGCCCGACGTATTTTTCCATCGCCGCAACAGCATTTTCGTGCGGCCAGCTGTTGGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATGCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGc         <  2:10188/247‑1 (MQ=255)
                                                                                                                                                                                                                                      cGAACAAATTCCTATACCAACTCGACGGATTTTCCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAAATGATCGCCGCCGAAAAAAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGc         <  2:461053/251‑1 (MQ=255)
                                                                                                                                                                                                                                             atTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAACCCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGAATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCTTGCTGTTGCGATCACATGCCAGcgcc  >  2:44084/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
TTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTGAAGCCGGCACGAACATATTCCTTTACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCC  >  NC_000913/2175749‑2176235

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: