Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 75,833 A→G A489A (GCT→GCC sgrR ← DNA‑binding transcriptional dual regulator SgrR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091375,8330AG96.4% 175.9 / ‑3.6 55A489A (GCT→GCCsgrRDNA‑binding transcriptional dual regulator SgrR
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (22/31);  total (22/33)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCGTGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTAGCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTCGATCTCTTTGATTTTCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTG  >  NC_000913/75588‑76072
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
gcgtGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTc                                                                                                                                                                                                                                             <  1:708454/251‑1 (MQ=255)
 cgtGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCg                                                                                                                                                                                                                                            <  1:195076/251‑1 (MQ=255)
    ggCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCGTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCt                                                                                                                                                                                                                                         <  1:678765/251‑1 (MQ=255)
            gTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAAGCAGGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGGATTCATGCGTAGGCCGCGCATACTGCTTTGCCCCTGAATGATCAGCCAGTGGTGGATAAGCTGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCGAGGTGCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTc                                                                                                                                                                                                                                 <  1:643490/251‑1 (MQ=255)
                agaACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTcgatc                                                                                                                                                                                                                             <  1:98900/251‑5 (MQ=255)
                  aaCGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTcgatcgg                                                                                                                                                                                                                           >  1:121339/1‑245 (MQ=255)
                              gtagtgATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGCTTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTTCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATTCAGtgttg                                                                                                                                                                                                               >  2:12260/1‑251 (MQ=255)
                                      ttttttATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAgc                                                                                                                                                                                                        <  2:576214/246‑2 (MQ=255)
                                                aCCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCa                                                                                                                                                                                             <  1:548215/251‑1 (MQ=255)
                                                aCCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCa                                                                                                                                                                                             <  1:50156/251‑1 (MQ=255)
                                                              agccGGTGGCGCAATCCCCGTGGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGaaca                                                                                                                                                                               <  2:158651/249‑3 (MQ=255)
                                                              aTCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGaac                                                                                                                                                                                >  1:206706/1‑249 (MQ=255)
                                                                       cgcAAACCACGCTGATTGAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGaacagcgaaaaat                                                                                                                                                                      <  1:489903/251‑12 (MQ=255)
                                                                             ccACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGaacagcgaaaaatccagcg                                                                                                                                                                >  1:449298/1‑234 (MQ=255)
                                                                                         aaaaTCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGCGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGt                                                                                                                                                   >  1:246722/1‑251 (MQ=255)
                                                                                             tCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTgac                                                                                                                                               >  1:526431/1‑249 (MQ=255)
                                                                                              cGAACCAGCCGTGGGGATTCATGCGTAGGCCGCGCCTACTGCGTTGCCCCTTAATGAGCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGgcg                                                                                                                                              <  2:653597/251‑1 (MQ=255)
                                                                                                 cccAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTg                                                                                                                                           <  2:595788/250‑1 (MQ=255)
                                                                                                 aCCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCCCCTCGCACAAGTGCGCGAACAGCGGAAAATCCAGCGGCAGCGTAAAGTTGGCGCTg                                                                                                                                           >  2:167223/1‑251 (MQ=255)
                                                                                                     gCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAAGACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGaacagcgaaaaatcc                                                                                                                                                                    <  1:76857/223‑14 (MQ=255)
                                                                                                           ggTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCGCCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGcc                                                                                                                                    >  1:205906/1‑248 (MQ=255)
                                                                                                                     cgtaGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGACGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGatagcgctctc                                                                                                                       >  2:536089/4‑243 (MQ=255)
                                                                                                                       cAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTAGCCAGATATCACTTTCGa                                                                                                                     <  2:431276/251‑1 (MQ=255)
                                                                                                                        aGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCa                                                                                                                                   >  1:143648/1‑236 (MQ=255)
                                                                                                                        aGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCa                                                                                                                                   <  2:143648/236‑1 (MQ=255)
                                                                                                                                 aTACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCCGATATCGCTCTCGATCcctccttc                                                                                                           >  2:551284/1‑249 (MQ=255)
                                                                                                                                           gCCCCTGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGa                                                                                                 <  1:162673/251‑1 (MQ=255)
                                                                                                                                           gCCCCGGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTg                                                                                                  <  1:709384/250‑1 (MQ=255)
                                                                                                                                                tGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCTCTGttaa                                                                                                                                       <  2:48133/208‑3 (MQ=255)
                                                                                                                                                tGAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGttaa                                                                                                                                       >  1:48133/1‑206 (MQ=255)
                                                                                                                                                 gAATGATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTTATTTCGCCAGTGCGCCAGCGCGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTa                                                                                            <  1:113786/250‑1 (MQ=255)
                                                                                                                                                     gATCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCCCTCCTTCATGCCACTGAGCGTAgcgg                                                                                        <  2:152535/250‑4 (MQ=255)
                                                                                                                                                      aTCAGCCAGTGGTGGATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGGGTTGGAGCAGCGGCCCCTCGCACACGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGGTGGCGCTGTTAAGCCAGAAATCGCTCTAGATCTCTCCTTCCTGCCACTGATCGTAgctgat                                                                                      >  2:519518/1‑246 (MQ=255)
                                                                                                                                                                  gtggataagcggCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGAttt                                                                           <  1:17773/241‑1 (MQ=255)
                                                                                                                                                                   tggataagcggCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCTTTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTCAGCCAGATATCGCTCTCGATCTCTCCTTAATGCCACTGATCGTAGCTGTTCTCTTTGAtttc                                                                          >  1:125343/10‑250 (MQ=255)
                                                                                                                                                                            ggCACCATTGCTTGGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGCTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGAGCGTAGCTGATCTCTTTGATTTCCAGCGTGAc                                                                 <  2:121576/251‑1 (MQ=255)
                                                                                                                                                                             gCACCAACGCTTTGCTGGCGACCAGTTGCTGGCCCCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGGGGCACTTCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTGAGCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTCGATCTCTTTGATTTTCAGCGTGAcc                                                                <  2:249012/251‑1 (MQ=255)
                                                                                                                                                                              cACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCt                                                               <  1:217257/251‑1 (MQ=255)
                                                                                                                                                                               aCCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGCTCGGAATCCAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGAtctc                                                                                                                 >  2:406102/1‑200 (MQ=255)
                                                                                                                                                                               aCCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCGCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGAtctc                                                                                                                 <  1:406102/200‑1 (MQ=255)
                                                                                                                                                                                  attgCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATAGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTg                                                           >  2:513627/4‑251 (MQ=255)
                                                                                                                                                                                           cTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGAGGCACCTCGCACAAGTGCGCGAACAGCGCAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATAACGCTCTCGATCCCTCCTTCATGCCACTGATCGTAGCTAATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGcccg                                                  >  2:541533/1‑249 (MQ=255)
                                                                                                                                                                                                   cAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATctgct                                          <  1:519518/251‑1 (MQ=255)
                                                                                                                                                                                                       tGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGCTTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCAt                                      >  2:560336/1‑251 (MQ=255)
                                                                                                                                                                                                       ggcgggCACCAGTTCGCCAGGTTCATTTCGCCAGTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCCCCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCAt                                      <  2:635438/247‑1 (MQ=255)
                                                                                                                                                                                                                    tcgccagggtcattgcGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCTCCTCGCACAAGCGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAAt                          <  2:246722/235‑1 (MQ=255)
                                                                                                                                                                                                                        cgggttcattccGCCATTGCGCCAGCGGGCGGCGGCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAAGCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCAccc                     <  2:277345/240‑1 (MQ=255)
                                                                                                                                                                                                                         aggttCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCg                    >  1:350258/4‑251 (MQ=255)
                                                                                                                                                                                                                               atttcGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTc              >  1:703806/4‑251 (MQ=255)
                                                                                                                                                                                                                                   cGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTccatg          >  2:143041/1‑247 (MQ=255)
                                                                                                                                                                                                                                   cGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACtg          >  1:58436/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aTTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGAt      <  1:317737/251‑1 (MQ=255)
                                                                                                                                                                                                                                         tGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATc     <  1:455182/250‑1 (MQ=255)
                                                                                                                                                                                                                                         ggcgcCAGCGGGCGGGGGCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCGGCAGCGTAAATTTGGCGCTGGTAAGCCAGATATCGCTCTCTATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCTCCCGATGCTCACTGTGATcc    <  2:13816/250‑1 (MQ=255)
                                                                                                                                                                                                                                           cgcCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAACAGCGAAAAATCCAGCTGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTg  <  1:210514/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
GCGTGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAAAACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTAGCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTCGATCTCTTTGATTTTCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCGATGCTCACTGTGATCCTG  >  NC_000913/75588‑76072

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