Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,041,102 A→G 7.6% N207S (AAC→AGC)  yedZ → inner membrane heme subunit for periplasmic YedYZ reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,041,1020AG 7.6% 163.8 / 2.9 52N207S (AAC→AGC) yedZinner membrane heme subunit for periplasmic YedYZ reductase
Reads supporting (aligned to +/- strand):  ref base A (20/28);  new base G (2/2);  total (22/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

GCTTTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATT  >  NC_000913/2040870‑2041339
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                             
gcttTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAAtt                                                                                                                                                                                                                                         <  1:267740/239‑1 (MQ=255)
gcttTAGCGTTCAATTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAAtt                                                                                                                                                                                                                                         <  1:256401/239‑1 (MQ=255)
                        gCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAATCAGGTTAACAAAAAATAATCTGATTACATTCATATTTccac                                                                                                                                                                                                                 >  2:868459/1‑237 (MQ=255)
                        gCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCg                                                                                                                                                                                                                  >  1:624621/1‑238 (MQ=255)
                        gCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGa                                                                                                                                                                                                                 >  1:45174/1‑239 (MQ=255)
                         gtattccggcaaaaCCGGGCAAGCACTGGCAACATTTTGATAATTTCGTCTATCTTGTCGCTATTCTGGCCCCAATTCATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGa                                                                                                                                                                                                                 <  2:192119/228‑1 (MQ=255)
                         gggtggggaaaaaaCGCGGAAAGCACTGGCAACATTTGGATATTTTCGTCTAGTTGGGCGCGATTCGGGCCCCAATACATTATGTGTGGTCCGTGAAGATTATTGCACCGCAGCTGCTCATATACGCTGGGCTGGGTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGa                                                                                                                                                                                                                 <  2:78190/228‑1 (MQ=255)
                            tGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACaaa                                                                                                                                                                                                             >  2:198713/1‑239 (MQ=255)
                               ctcaaaaaCCCGGCCAGGCCGGGCAAACGTTTGATTATTTCGTCTTACTTGTCGCGAATCTTGCCCCAAATACATATCTGTGGTCCGTGAAGATTATTTCCCCCCAGGCGCTCATCTACGCTGGGGTTGGTGGACTGGCTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTa                                                                                                                                                                                                          <  2:34136/235‑1 (MQ=255)
                                    aaaCTCGGCAAGCCCTTGCAACCGTGTTATAATTTCGTCTATCTTGTCCCGATTTCGGCCCCAATTCATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCAACTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATATCTTCATATTGCCGACAAAGTAcgcg                                                                                                                                                                                                      <  2:589467/238‑1 (MQ=255)
                                     aaCTCGGCAAGCCCTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTAcgc                                                                                                                                                                                                       <  1:651860/236‑1 (MQ=255)
                                               gCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGATTAACTTCATATTCCCGACAAAGTACGCGTTAACCCTg                                                                                                                                                                                             >  2:651860/1‑236 (MQ=255)
                                                 aCTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGGTCACAACCAATTATCTGTATAACTTCTTATTTCCGACAAAATACGCCGTAACCCTCCCGGa                                                                                                                                                                                        >  2:85831/1‑239 (MQ=255)
                                                                ttatttttgttttttttGTCGGGTTTTTGGCCCCAATTCATTATGTGTGGACCGTGAAGGTTATTTCACCGCCGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGcc                                                                                                                                                                         <  2:217795/224‑1 (MQ=255)
                                                                    tttcctttttttttttcggtttCTGGCCCCAATTCCTTTCCTGTGGTCCGGGAAGATTATTTCCCCGCAGCCGTTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTAt                                                                                                                                                                     <  2:131998/219‑1 (MQ=255)
                                                                    tttCGTTTATCTTGTCGCGATTCTGGCCCCAATACATTAACTGGGGGCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTAt                                                                                                                                                                     <  1:248800/239‑1 (MQ=255)
                                                                    tttCGTCTTTCTTGTCGCGATTCAGGCCCCAATACATTTTCTGGGGTCCGTGAAGATTTTTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTAt                                                                                                                                                                     <  2:689370/239‑1 (MQ=255)
                                                                         tCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGc                                                                                                                                                                >  1:841931/1‑239 (MQ=255)
                                                                         tCTATCTTGTCGCGATTCAGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGc                                                                                                                                                                >  1:689370/1‑239 (MQ=255)
                                                                                       gttttGGCCCCAATACATTATCTGGGGTCCGTGAAATTTATTTCCCCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCa                                                                                                                                                  <  2:624621/235‑1 (MQ=255)
                                                                                       aTTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTTCCGACAAAGt                                                                                                                                                                                                           >  2:365789/1‑182 (MQ=255)
                                                                                       aTTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGt                                                                                                                                                                                                           <  1:365789/182‑1 (MQ=255)
                                                                                        ttCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCa                                                                                                                                                  <  1:301161/238‑1 (MQ=255)
                                                                                                           aTCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGa                                                                                                                                                                                                                 >  1:743409/1‑157 (MQ=255)
                                                                                                           aTCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGa                                                                                                                                                                                                                 <  2:743409/157‑1 (MQ=255)
                                                                                                             cTGTGGTCCGTGAAGATTATTTCAGCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAAtt                                                                                                                                                                                                                                         >  1:869139/1‑131 (MQ=255)
                                                                                                             cTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGTTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAAtt                                                                                                                                                                                                                                         <  2:869139/131‑1 (MQ=255)
                                                                                                                tgGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAGCAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGtcatc                                                                                                                          >  1:305387/1‑239 (MQ=255)
                                                                                                                tgGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGCCAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGtcatc                                                                                                                          >  1:734191/1‑239 (MQ=255)
                                                                                                                 ggTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATaca                                                                                                                                                           <  1:129831/205‑1 (MQ=255)
                                                                                                                 ggTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATaca                                                                                                                                                           >  2:129831/1‑205 (MQ=255)
                                                                                                                 ggTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGATATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGAGGCGCGAGCGCCCTATCCGGCATACCCAGATCGCATCAGGTTATGGTTCGTTTGtcaaaa                                                                                                                         >  2:354248/1‑236 (MQ=255)
                                                                                                                        gAAGATTATTTCCCCGAAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCaa                                                                                                                  <  1:354248/239‑1 (MQ=255)
                                                                                                                          agtttattTCACCGCAGCCGCTCAACTACGCTGGGCTGGTTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGTTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAAt                                                                                                                <  1:868459/236‑1 (MQ=255)
                                                                                                                                      ccctcccGTTCTTCTACGGTGGGCTGGATTTACTTCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCGCAAATCACCAAaataat                                                                                                    <  2:734191/234‑1 (MQ=255)
                                                                                                                                                         tgggtttggttttcTGCTTTTTGCCTTACGTTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCTTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAaga                                                                                 <  2:675788/226‑1 (MQ=255)
                                                                                                                                                           ggctggctGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAACGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGCCATACACAGATCCCACTAGGTTATGTTTCGTTTGTCATCAACCTCAAATACCCAAAATAAATTATTTTATTTCAATAagaaa                                                                               >  2:99144/1‑237 (MQ=255)
                                                                                                                                                              ttgttttactttttttAGCCTTACGGTATAAGAAGTTGCTTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAGCAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAagatagat                                                                            <  2:305387/227‑1 (MQ=255)
                                                                                                                                                                   ttACTGTTTTTAGCCTTACGGTATAAAAATTTTCGTTCCCTTTTTAACCGGTTCCCCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGt                                                                        <  2:841931/237‑1 (MQ=255)
                                                                                                                                                                              aGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAAtttattttat                                                                                           >  2:433989/1‑208 (MQ=255)
                                                                                                                                                                              aGCCTTACGGTATAAGAAGTTGCGTTCCCTTTGTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAAtttattttat                                                                                           <  1:433989/208‑1 (MQ=255)
                                                                                                                                                                                   tACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAGCAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGAAAAGTAGAACTGAGaaaaacag                                                       >  2:784154/1‑234 (MQ=255)
                                                                                                                                                                                                  tGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGa                                                                                                                                                                              <  2:104108/105‑1 (MQ=255)
                                                                                                                                                                                                  tGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGa                                                                                                                                                                              >  1:104108/1‑105 (MQ=255)
                                                                                                                                                                                                  tGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTAt                                                                                                                                                                     <  2:574571/114‑1 (MQ=255)
                                                                                                                                                                                                  tGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTAt                                                                                                                                                                     >  1:574571/1‑114 (MQ=255)
                                                                                                                                                                                                  tGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCCACATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCtca                                        >  1:35731/1‑237 (MQ=255)
                                                                                                                                                                                                  tGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAAAACTGAAAAAGCCATTCTCTCGTTTCCTaa                                        >  2:86148/1‑239 (MQ=255)
                                                                                                                                                                                                                    ggTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCt                      <  1:198713/239‑1 (MQ=255)
                                                                                                                                                                                                                            aaacatgTTCCCAACAAATTTACGTTATAACTTAAATTTGCCGAAAAAGTACGCGTTATCCCTGCCGGATGGGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTa              <  2:35731/233‑1 (MQ=255)
                                                                                                                                                                                                                            aaGGAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGGACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTa              <  1:85831/239‑1 (MQ=255)
                                                                                                                                                                                                                                   ttCACAGCAAATTATCTGTATACCTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAata         <  1:784154/237‑1 (MQ=255)
                                                                                                                                                                                                                                        aacaaATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTTTTGCATTTTGCTATATGTTCAAATAAAACAtt  >  2:155624/1‑239 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                             
GCTTTAGCGTTCACTTCAACCCAGGCGATGCAGCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATT  >  NC_000913/2040870‑2041339

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