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

TTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTTTACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCC  >  NC_000913/2175749‑2176235
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
ttCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCACTCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTCTGGGTTCCCCCGCGCAGGACATTGACGCATCAAGATGACTTTTACTG‑ACGCCGGCACGAACATATTCCTATACCa                                                                                                                                                                                                                                                 >  2:460128/1‑248 (MQ=255)
ttCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGTAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTATACCa                                                                                                                                                                                                                                                 <  1:460128/248‑1 (MQ=255)
      aggggTGATGGGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTATACCAGCTCGACgg                                                                                                                                                                                                                                        <  2:75035/249‑1 (MQ=255)
      aTGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGCTCCCCCGCGCAGGACAGTGACGCATCAAGATGAATTTTACTG‑TAGCCGGCCCGAACATTTTCCTTTCACAGCACGACgg                                                                                                                                                                                                                                        >  2:96713/1‑251 (MQ=255)
      aTGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTATACCAGCTCGACgg                                                                                                                                                                                                                                        >  1:414167/1‑251 (MQ=255)
                          gaATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTTGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGccgcagc                                                                                                                                                                                                                     <  1:37923/250‑3 (MQ=255)
                          gaATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGccgcatcc                                                                                                                                                                                                                    >  1:393539/1‑251 (MQ=255)
                                       atCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCACTCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACTTTGACGCATCAAGATGAATTTTACTG‑ATGCCGGCACGAACATATTCCTATACCAGCTCGACGGTTTTTTACATCGCCGCATCCGCATTTTCTtgct                                                                                                                                                                                                       >  2:208812/2‑251 (MQ=255)
                                        tCACCGCCCGGACCCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTtgctg                                                                                                                                                                                                      >  1:117709/1‑251 (MQ=255)
                                        tCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTtgctg                                                                                                                                                                                                      >  1:10188/1‑251 (MQ=255)
                                                  gaaccggaacTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCACAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTATACCAGCTCg                                                                                                                                                                                                                                            >  1:410725/1‑203 (MQ=255)
                                                  gaaccggaacGTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTATACCAGCTCg                                                                                                                                                                                                                                            <  2:410725/203‑1 (MQ=255)
                                                   aaccggaacTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTg                                                                                                                                                                                           <  1:230699/251‑1 (MQ=255)
                                                             tCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAgat                                                                                                                                                                                  <  1:215814/250‑1 (MQ=255)
                                                                    cAATGACATAGCTCAGTTGCTCACGCTGGCAAGCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATcgccg                                                                                                                                                                          <  2:129095/251‑1 (MQ=255)
                                                                       tGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGccgccg                                                                                                                                                                       >  1:193528/1‑251 (MQ=255)
                                                                             agatCAGTTGCTCACACTGGCAATCTGACGCAACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCCGGACATTGACGCATCAAGATGAATTTTACTG‑CAGCCGGCACGAACATATTCCTATACAAGCTCGACGGATTTTTCTATCGCCCCATCCCCATTTTCTTGCTGCCAGCAGTGGGGCCCCAgctgaacaccaccgatact                                                                                                                                                                  >  2:313126/4‑236 (MQ=255)
                                                                             aGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGaat                                                                                                                                                                  <  1:265177/250‑1 (MQ=255)
                                                                                 cAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGTGCCCAGATGATCGCCGCCGAGAATAATg                                                                                                                                                              <  1:185177/250‑1 (MQ=255)
                                                                                      gCTCATGCTGGCAATCTGTCGCCACACTTTCCGCAGAAAAGCAAAGAACAGCAGCTCGTTCCGCAACCTTTTCTGGTGCTAACTGTATGGGATCCCCCGCGAAGGACATTGACGCAGCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATCCCTTTACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTc                                                                                                                                                         <  2:166702/250‑1 (MQ=255)
                                                                                         cACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGt                                                                                                                                                     <  1:390264/251‑1 (MQ=255)
                                                                                                    tCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAAc                                                                                                                                          >  1:115690/1‑251 (MQ=255)
                                                                                                     ttGTCGCCATACTTTCCGCGGCAAAGCAAAGCACAGCATCTCGTTCCGCATCCGTTTCTGTTGTTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGAGGAATTTTACTG‑AAGCCTGCACGAACATATTCCGATACCAGCTCGACGGATTTTTCCATCGCCGCAGCCGCATTTTCTTGCTGCCAGCGGTTTGGCCCCATATGATCGCCGCCGATAATAATGCGTGCGCGTGCAAACCCAAc                                                                                                                                          <  2:177766/249‑1 (MQ=255)
                                                                                                         cGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTAt                                                                                                                                     <  1:299434/251‑1 (MQ=255)
                                                                                                                    ccGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTAGCGGCAATCGTa                                                                                                                          <  2:409349/251‑1 (MQ=255)
                                                                                                                       cagcaAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTaaa                                                                                                                        >  1:155738/1‑250 (MQ=255)
                                                                                                                          caaagcaaagCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAAc                                                                                                                     >  1:461053/1‑250 (MQ=255)
                                                                                                                              gcaaagccCAGCAGCTCGTTCCGCAACTGTTTCTGGTGCTAACGGTATTGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTGCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAAtt                                                                                                                <  2:320650/251‑1 (MQ=255)
                                                                                                                                              cgtgccGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGa                                                                                                           <  1:189437/236‑1 (MQ=255)
                                                                                                                                              cGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGAACTTGACGCATCAAGATGAATTTTACTG‑ACGCCGGCCCGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGAATCCGCATTTTCTTGCTGCCAGCAATTTGGCCCCAGATGATCGCCGCCGAGCATAATGCGTTCCCGTGCAAACCCAACTTTATCGGCCATCGTAAAAACCAATTCGCGa                                                                                                           >  2:189437/1‑240 (MQ=255)
                                                                                                                                                 tcccccaCCGTTTCTGGTGCTAACGGTATGTGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTAGACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGgtg                                                                                              <  2:117709/244‑1 (MQ=255)
                                                                                                                                                   cGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTc                                                                                            >  1:88732/1‑250 (MQ=255)
                                                                                                                                                                 gTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCg                                                                             >  1:394459/1‑251 (MQ=255)
                                                                                                                                                                            gtGGGACCCCCCGCGCAGGACATTGACGCATCAAGAAGAATTTTACTG‑AAGCCGGCACGAACATATTCCAATCCCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGAAGTTTGGCCCCAGATGATCGCAGCCGAGAATAATGCGGTCGCGTGCAAACCCAACTTTACCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGCCATTCCGGGATAACCGCCAAattgatt                                                                  <  2:21286/250‑1 (MQ=255)
                                                                                                                                                                                gatcccccctccAGGACCTTGTCGCATCAAGAAGAATTTGACTG‑ACGCCGGCACGAACATATTCCTATCCCAGCCCGATGGCTTTTTCCATCGCCGCATCCGTATTTTCTGGCTCCCAGCAGTTTGGCCCCAGATGATCTCCGCCGAGAATAATGCGTTCGCGTGCAAACCCACCTTTATCGGCAATCTTAAAAACAAATTAGCGAAAGTCGGCGGGTGTCATCCCGGTATAACCGCCAAATTGATTGAcc                                                              <  2:216768/240‑1 (MQ=255)
                                                                                                                                                                                                    gACGCATCAAGATGAATTTTA‑TGTAAGCCGGCACGAACATATTCCTATACCAGCCGGACTGATTTTCCCATCGCCGCTTCCGCATCTTCTTGCTGCCAGCAGTTTGGCCCCAGATGACCCCCGCCGAGAATAATGCGTTCGCTGGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCTGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAAt                                          <  2:197301/251‑1 (MQ=255)
                                                                                                                                                                                                                          tG‑AAGCCGGCACGAACATATTCCTATACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCGGCACTTTGCGCGTGCTGTTg                    >  1:227022/1‑251 (MQ=255)
                                                                                                                                                                                                                           g‑AAGCCGGCACGAACATATTCCTATCCCAGCTCAACGGATTTTCCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGGTCGCGTGCAAACCCAACTTTACCGGCAATCGTAAAAACAAAGTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGtt                     <  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)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
TTCAACATGGGTGATGTGTACTGACTGAATGGCGCTGGCCTCACCGCCCGGAACCGGAACTTCGGTGCCAATGACATAGCTCAGTTGCTCACGCTGGCAATCTGTCGCCACACTTTCCGCAGCAAAGCAAAGCACAGCAGCTCGTTCCGCAACCGTTTCTGGTGCTAACGGTATGGGATCCCCCGCGCAGGACATTGACGCATCAAGATGAATTTTACTG‑AAGCCGGCACGAACATATTCCTTTACCAGCTCGACGGATTTTTCCATCGCCGCATCCGCATTTTCTTGCTGCCAGCAGTTTGGCCCCAGATGATCGCCGCCGAGAATAATGCGTTCGCGTGCAAACCCAACTTTATCGGCAATCGTAAAAACAAATTCGCGAAAGTCTGCCGGTGTCATTCCGGTATAACCGCCAAATTGATTGACCTGGTTTGACGTTGCTTCAATCAGCACTTTGCGCGTGCTGTTGCGATCAAATGCCAGCGCC  >  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: ‑