Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 75,806 A→G N498N (AAT→AAC)  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,8060AG94.4% 156.0 / 0.7 54N498N (AAT→AAC) sgrRDNA‑binding transcriptional dual regulator SgrR
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base G (24/27);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

GGGTATTTTCTCAGCGCACGCGTACGCGTGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑CAATAA‑TGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAA‑AACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTAGCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTCGATCTCTTTGATTTTCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCAT  >  NC_000913/75563‑76036
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                       
gggTATTTTCTCAGCGCACGCGTACGCGTGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGATTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCCGGTTCAttt                                                                                                                                                                                                                                  >  2:541119/1‑249 (MQ=255)
 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:50156/251‑1 (MQ=255)
                                                                         aCCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCa                                                                                                                                                         <  1:548215/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)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                       
GGGTATTTTCTCAGCGCACGCGTACGCGTGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑CAATAA‑TGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAA‑AACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTAGCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTCGATCTCTTTGATTTTCAGCGTGACCTGGTGACTTGCCAGAATCTGCTGCAT  >  NC_000913/75563‑76036

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: ‑