Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 75,812 C→T E496E (GAG→GAA)  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,8120CT96.2% 166.6 / ‑3.4 53E496E (GAG→GAA) sgrRDNA‑binding transcriptional dual regulator SgrR
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (24/27);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

GGTATTTTCTCAGCGCACGCGTACGCGTGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑CAATAA‑TGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAA‑AACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTAGCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTCGATCTCTTTGATTTTCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCG  >  NC_000913/75564‑76054
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   
ggTATTTTCTCAGCGCACGCGTACGCGTGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑CAATAA‑TGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCAtctc                                                                                                                                                                                                                                                   <  1:357616/251‑1 (MQ=255)
 gTATTTTCTCAGCGCACGCGTACGCGTGTCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCg                                                                                                                                                                                                                                                  <  1:157712/251‑1 (MQ=255)
         tcAGCGCACGCGTACGCGTGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTgcg                                                                                                                                                                                                                                          >  2:548215/1‑251 (MQ=255)
             cgcACGAGTACGCGTGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGAACAGTTGCAGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGcc                                                                                                                                                                                                                                        >  1:607982/1‑249 (MQ=255)
                 cgcgtacgcgtGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCggg                                                                                                                                                                                                                                  <  1:421109/251‑3 (MQ=255)
                  gcgtacgcgtGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCgggc                                                                                                                                                                                                                                 >  1:277345/1‑248 (MQ=255)
                   cgtacgcgtGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCgggcg                                                                                                                                                                                                                                >  1:653597/1‑247 (MQ=255)
                    gtacgcgtGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGGATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGGCTA‑AG‑CGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCCCCAGTTCGCCCGGTTCATTTCGCCATTGCGCCAGCgggcgg                                                                                                                                                                                                                               >  2:50156/1‑246 (MQ=255)
                        gcgtGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTc                                                                                                                                                                                                                           <  1:708454/251‑1 (MQ=255)
                         cgtGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCg                                                                                                                                                                                                                          <  1:195076/251‑1 (MQ=255)
                            ggCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCGTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCt                                                                                                                                                                                                                       <  1:678765/251‑1 (MQ=255)
                                    gTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAAGCAGGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGGATTCATGCGTAGGCCGCGCATACTGCTTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCTGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCGAGGTGCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTc                                                                                                                                                                                                               <  1:643490/251‑1 (MQ=255)
                                        agaACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTcgatc                                                                                                                                                                                                           <  1:98900/251‑5 (MQ=255)
                                          aaCGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTcgatcgg                                                                                                                                                                                                         >  1:121339/1‑245 (MQ=255)
                                                      gtagtgATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGCTTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTTCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATTCAGtgttg                                                                                                                                                                                             >  2:12260/1‑251 (MQ=255)
                                                              ttttttATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAgc                                                                                                                                                                                      <  2:576214/246‑2 (MQ=255)
                                                                        aCCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCa                                                                                                                                                                           <  1:548215/251‑1 (MQ=255)
                                                                        aCCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCa                                                                                                                                                                           <  1:50156/251‑1 (MQ=255)
                                                                                      agccGGTGGCGCAATCCCCGTGGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGaaca                                                                                                                                                             <  2:158651/249‑3 (MQ=255)
                                                                                      aTCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGaac                                                                                                                                                              >  1:206706/1‑249 (MQ=255)
                                                                                               cgcAAACCACGCTGATTGAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGaacagcgaaaaat                                                                                                                                                    <  1:489903/251‑12 (MQ=255)
                                                                                                     ccACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGaacagcgaaaaatccagcg                                                                                                                                              >  1:449298/1‑234 (MQ=255)
                                                                                                                 aaaaTCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGCGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGt                                                                                                                                 >  1:246722/1‑251 (MQ=255)
                                                                                                                     tCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGTTgac                                                                                                                             >  1:526431/1‑249 (MQ=255)
                                                                                                                      cGAACCAGCCGTGGGGATTCATGCGTAGGCCGCGCCTACTGCGTTGCCCCTTAATGAGCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGTTGgcg                                                                                                                            <  2:653597/251‑1 (MQ=255)
                                                                                                                         cccAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTg                                                                                                                         <  2:595788/250‑1 (MQ=255)
                                                                                                                         aCCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCCCCTCGCACAAGTGCGCGAACAGCGGA‑AAATCCAGCGGCAGCGTAAAGTTGGCGCTg                                                                                                                         >  2:167223/1‑251 (MQ=255)
                                                                                                                             gCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAAGACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGaacagcgaaaaatcc                                                                                                                                                  <  1:76857/223‑14 (MQ=255)
                                                                                                                                   ggTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCGCCTCGCACAAGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGcc                                                                                                                  >  1:205906/1‑248 (MQ=255)
                                                                                                                                             cgtaGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGACGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGatagcgctctc                                                                                                     >  2:536089/4‑243 (MQ=255)
                                                                                                                                               cAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑CAATAA‑TGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAA‑AACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTAGCCAGATATCACTTTCGa                                                                                                   <  2:431276/251‑1 (MQ=255)
                                                                                                                                                aGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCa                                                                                                                 <  2:143648/236‑1 (MQ=255)
                                                                                                                                                aGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCa                                                                                                                 >  1:143648/1‑236 (MQ=255)
                                                                                                                                                         aTACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCCGATATCGCTCTCGATCcctccttc                                                                                         >  2:551284/1‑249 (MQ=255)
                                                                                                                                                                   gCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGa                                                                               <  1:162673/251‑1 (MQ=255)
                                                                                                                                                                   gCCCCGGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTg                                                                                <  1:709384/250‑1 (MQ=255)
                                                                                                                                                                        tGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCTCTGttaa                                                                                                                     <  2:48133/208‑3 (MQ=255)
                                                                                                                                                                        tGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGttaa                                                                                                                     >  1:48133/1‑206 (MQ=255)
                                                                                                                                                                         gAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTTATTTCGCCAGTGCGCCAGCGCGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTa                                                                          <  1:113786/250‑1 (MQ=255)
                                                                                                                                                                             gATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCCCTCCTTCATGCCACTGAGCGTAgcgg                                                                      <  2:152535/250‑4 (MQ=255)
                                                                                                                                                                              aTCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGGGTTGGAGCAGCGGCCCCTCGCACACGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGGTGGCGCTGTTAAGCCAGAAATCGCTCTAGATCTCTCCTTCCTGCCACTGATCGTAgctgat                                                                    >  2:519518/1‑246 (MQ=255)
                                                                                                                                                                                           gtggataagcggCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGAttt                                                         <  1:17773/241‑1 (MQ=255)
                                                                                                                                                                                            tggataagcggCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCTTTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTCAGCCAGATATCGCTCTCGATCTCTCCTTAATGCCACTGATCGTAGCTGTTCTCTTTGAtttc                                                        >  1:125343/10‑250 (MQ=255)
                                                                                                                                                                                                     ggCACCATTGCTTGGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGCTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGAGCGTAGCTGATCTCTTTGATTTCCAGCGTGAc                                               <  2:121576/251‑1 (MQ=255)
                                                                                                                                                                                                      gCACCAACGCTTTGCTGGCGACCAGTTGCTGGCCCCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGGGGCACTTCGCATAAATGTGCGAA‑AACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTGAGCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTCGATCTCTTTGATTTTCAGCGTGAcc                                              <  2:249012/251‑1 (MQ=255)
                                                                                                                                                                                                       cACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCt                                             <  1:217257/251‑1 (MQ=255)
                                                                                                                                                                                                        aCCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGCTCGGAATCCAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGAtctc                                                                                               >  2:406102/1‑200 (MQ=255)
                                                                                                                                                                                                        aCCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCGCGCACAAGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGAtctc                                                                                               <  1:406102/200‑1 (MQ=255)
                                                                                                                                                                                                           attgCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATAGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTg                                         >  2:513627/4‑251 (MQ=255)
                                                                                                                                                                                                                    cTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGAGGCACCTCGCACAAGTGCGCGAA‑CAGCGCAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATAACGCTCTCGATCCCTCCTTCATGCCACTGATCGTAGCTAATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGcccg                                >  2:541533/1‑249 (MQ=255)
                                                                                                                                                                                                                            cAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATctgct                        <  1:519518/251‑1 (MQ=255)
                                                                                                                                                                                                                                tGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGCTTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCAt                    >  2:560336/1‑251 (MQ=255)
                                                                                                                                                                                                                                ggcgggCACCAGTTCGCCAGGTTCATTTCGCCAGTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCCCCTCGCACAAGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCAt                    <  2:635438/247‑1 (MQ=255)
                                                                                                                                                                                                                                                  aggttCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCAGATATCGCTCTCGATCTCTCCTTCATGCCACTGATCGTAGCTGATCTCTTTGATTTCCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCg  >  1:350258/4‑251 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   
GGTATTTTCTCAGCGCACGCGTACGCGTGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑CAATAA‑TGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAA‑AACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTAGCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTCGATCTCTTTGATTTTCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCATGATCCCGGCAATCACCCG  >  NC_000913/75564‑76054

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