Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 75,758 1 bp→GC coding (1542/1656 nt) 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,7571.G92.6% 167.7 / 3.0 54coding (1543/1656 nt)sgrRDNA‑binding transcriptional dual regulator SgrR
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base G (24/26);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
*NC_00091375,7580TC96.3% 161.0 / ‑4.2 54P514P (CCA→CCG) sgrRDNA‑binding transcriptional dual regulator SgrR
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (26/26);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-01

CGGCGTTATCCGGATCAGGTTCGACGGGTATTTTCTCAGCGCACGCGTACGCGTGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑CAATAA‑TGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAA‑AACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTAGCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTCGATCTCTTTGATTTTCAGCGTGAC  >  NC_000913/75538‑76009
                                                                                                                                                                                                                              ||                                                                                                                                                                                                                                                           
cGGCGTTATCCGGATCAGGTTCGACGGGTATTTTCTCAGCGCACGCGTACGCGTGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGc                                                                                                                                                                                                                                <  1:525507/251‑1 (MQ=255)
   cGTTATCCGGATCAGGTTCGACGGGTATTTTCTCAGCGCACGCGTACGCGTGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTgg                                                                                                                                                                                                                             >  1:635438/1‑251 (MQ=255)
   cGTTATCCGGATCAGGTTCGACGGGTATTTTCTCAGCGCACGCGTACGCGTGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTgg                                                                                                                                                                                                                             >  1:602795/1‑251 (MQ=255)
     ttATCCGGATCAGGTTCGACGGGTATTTTCTCAGCGCACGCGTACGCGTGGCACCCCGTTGAGCACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGCTTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACAATTGCTTTGCTGTCTGCCAGTTGCTGGc                                                                                                                                                                                                                            >  2:312596/1‑250 (MQ=255)
          cggctCAGGTTCGACGGGTATTTTCTCAGCGCACGAGTACGCGTGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCACCCAGCAATCATGGCTCCGGTGGCGCAAACCACGCTGATTTAAAATCGACCCAGCCGAGGGTGTTCATGCTTAGGCCGCGCATACTGGGGTGCCCCTGAATGATCAGCCAGGGGTG‑‑GTTAAGCGGCACCATTGCTTTGCGGGCGACGCGTTGCTgaccccagg                                                                                                                                                                                                                      >  2:455182/5‑243 (MQ=255)
          cGGATCCGGGTCGGCGGGTATTTTCTCAGCGCACGCGTACGCGTGGCACCCCGTTGAGAACGGCGTTAGTGTAGTTATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGt                                                                                                                                                                                                                      <  2:33579/251‑1 (MQ=255)
          cGGATCAGGTTCGACGGGTATTTTCTCAGCGCACGCGTACGCGAGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCCCGCTGATTTAAAATCGCACCAGCCGAGGGAATTCATGCGTAGGCCACGCATACTGCGTTGCCCCTGAATGTTCAGCCAGTGGGG‑GAGAAG‑CGGCACCATTGCTTTGCTGGGGACCAGTAGCGGGCCCCAGt                                                                                                                                                                                                                      >  2:633511/1‑251 (MQ=255)
            gATCAGGTTCGACGGGTATTTTCTCAGCGCACGCGTACGCGTGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTc                                                                                                                                                                                                                    >  1:667825/1‑251 (MQ=255)
            gATCAGGTTCGACGGGTATTTTCTCAGCGCACGCGTACGCGTGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTc                                                                                                                                                                                                                    <  1:223587/251‑1 (MQ=255)
               cAGGTTCGACGGGTATTTTCTCAGCGCACGCGTACGCGTGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGAG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGcc                                                                                                                                                                                                                 <  1:639120/251‑1 (MQ=255)
                         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)
                                                acgcgtGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCa                                                                                                                                                                                                                <  1:271217/219‑1 (MQ=255)
                                                acgcgtGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGTAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTAGCGACCAGTTGCTGGCACCAGTTCGCCa                                                                                                                                                                                                                >  2:271217/1‑219 (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                                                                    >  1:143648/1‑236 (MQ=255)
                                                                                                                                                                          aGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑‑GATAAGCGGCACCATTGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGGTTCATTTCGCCATTGCGCCAGCGGGCGGCGTCGGCTTGCCAGTCGATCGGAATACAGTGTTGGAGCAGCGGCACCTCGCACAAGTGCGCGAA‑CAGCGAAAAATCCAGCGGCAGCGTAAAGTTGGCGCTGTTAAGCCa                                                                    <  2:143648/236‑1 (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)
                                                                                                                                                                                                                              ||                                                                                                                                                                                                                                                           
CGGCGTTATCCGGATCAGGTTCGACGGGTATTTTCTCAGCGCACGCGTACGCGTGGCACCCCGTTGAGAACGGCGTTAGTGTAGTGATTTTGTTATCAACCAGCAATCATGGATCCGGTGGCGCAAACCACGCTGATTTAAAATCGAACCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTG‑CAATAA‑TGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGCGCCAGCGAGCAGCGTCGGCTTGCCAGTCAATGGGAATGCAATGTTGTAGCAGTGGCACTTCGCATAAATGTGCGAA‑AACAGAGAAGTCCAGCGGCAGGGTAAAGTTGGCGCTGTTTAGCCAGATATCACTTTCGATCTCTCCTGTATGCCACTGATCGTAGTCGATCTCTTTGATTTTCAGCGTGAC  >  NC_000913/75538‑76009

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