Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,312,215 A→G T259T (ACT→ACC alsR ← DNA‑binding transcriptional repressor AlsR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,312,2150AG97.8% 162.9 / ‑4.9 45T259T (ACT→ACCalsRDNA‑binding transcriptional repressor AlsR
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (24/20);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCAT  >  NC_000913/4311970‑4312454
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
attttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGcaaggttggt                                                                                                                                                                                                                                             >  2:537251‑M2/88‑248 (MQ=255)
 ttttatgtggggtccttccctgatattacgggcactatttattcacaactctgccgaaaaacaggctgtcgtcagttttgacgtgaGGTAATGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGAGCCTGGTGAGAAGAAATCAACCATTGCGCCGGTTTTTTCCATATTAAGATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAggttggtt                                                                                                                                                                                                                                            >  2:166302‑M2/87‑251 (MQ=255)
             tcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTg                                                                                                                                                                                                                               >  2:171240‑M2/75‑251 (MQ=255)
              cattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGc                                                                                                                                                                                                                              <  1:89801‑M2/178‑1 (MQ=255)
                          ttacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACg                                                                                                                                                                                                                  >  2:111305‑M2/62‑251 (MQ=255)
                                     atttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACg                                                                                                                                                                                                       <  2:295170‑M2/201‑1 (MQ=255)
                                                 ctctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCgggg                                                                                                                                                                                          >  1:285417‑M2/39‑251 (MQ=255)
                                                     gacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTGCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTg                                                                                                                                                                                       <  1:430851‑M2/216‑1 (MQ=255)
                                                            aacacgctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATCTTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAgct                                                                                                                                                                                >  2:544250‑M2/28‑249 (MQ=255)
                                                                ggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAaataa                                                                                                                                                                           <  2:6429‑M2/228‑1 (MQ=255)
                                                                    gtcgtcagttttgacgtgaGGTAAGGGTACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGATGCATTACGACCTAATAACGGCGTTTCCGGTGCTGGTGAGCAaataatatg                                                                                                                                                                       <  2:549969‑M2/232‑2 (MQ=255)
                                                                        tcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCg                                                                                                                                                                   >  1:185317‑M2/16‑251 (MQ=255)
                                                                                            ggggACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAatat                                                                                                                                                                                                                                    >  2:90275/1‑167 (MQ=255)
                                                                                            ggggACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAatat                                                                                                                                                                                                                                    <  1:90275/167‑1 (MQ=255)
                                                                                             gggACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGa                                                                                                                                               >  1:381647/1‑250 (MQ=255)
                                                                                                    gTAACTCCGCGATTGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGACGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAg                                                                                                                                        <  2:438450/251‑1 (MQ=255)
                                                                                                    gTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGt                                                                                                                                                                                     >  2:311484/1‑206 (MQ=255)
                                                                                                    gTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGt                                                                                                                                                                                     <  1:311484/206‑1 (MQ=255)
                                                                                                            gcgATGGCGGAATTTCCTTTTTCTTATTGCAGCGCGCCAGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCAGGGCGACAGAGACGAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTGCGCTGTCGGTGAATGGTAGCTATGGGt                                                                                                                                <  2:204691/251‑1 (MQ=255)
                                                                                                                           ccTTTTTCTTATTTTAGCGCGGCCTGTGAGGAGGAATCCACAATTTCGCCGGTTTTTTCCATATTAATAATAGCTTGTTCAATGTTGAGCTGGGGGACAGGAACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTAGGGGTTATACAAATAATCt                                                                                                                  <  2:518494/250‑1 (MQ=255)
                                                                                                                                    tATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCTGGTTTTTTCCATATTACTATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAACAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGACGAGGCATTAAGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGCTCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGtt                                                                                                        >  2:83762/1‑251 (MQ=255)
                                                                                                                                          aGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCttttt                                                                                                  <  2:12206/251‑1 (MQ=255)
                                                                                                                                            cgcgcCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTAAGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTg                                                                                                >  2:52043/1‑251 (MQ=255)
                                                                                                                                                            gaaATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGAAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCg                                                                                 >  1:259513/1‑251 (MQ=255)
                                                                                                                                                                tCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCttt                                                                             >  2:388671/1‑251 (MQ=255)
                                                                                                                                                                  aacaaTTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCtttt                                                                            >  1:217614/1‑250 (MQ=255)
                                                                                                                                                                     aatggggcCGGTTTTTTCCATATAAATATTAGCTGGTTCAAGGTTGAGCTGGGCGAAAGAGACAAAAAAAGCTCCAAGAAAGGTTGGTTGTAATATCCTGGCGGAGGAATTCCGCCCTAAAACGGGCGTTTCCGGGGCTGGTGGGCAAAAAATATGACGGGCCGGTTCGGCTACCGTTGAATGGTAGCAATGGGTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATc                                                                       <  2:185317/245‑1 (MQ=255)
                                                                                                                                                                           gcCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTgg                                                                  >  1:293748/1‑251 (MQ=255)
                                                                                                                                                                                        atattaatattaGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATg                                                     >  2:143143/1‑251 (MQ=255)
                                                                                                                                                                                        atattaatattaGCTTGATCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCGATGGGTTATACAAATAATCGTTGCCCGGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATg                                                     >  1:567889/1‑251 (MQ=255)
                                                                                                                                                                                                attaGCTTGTTCAATGTTGAGCTGGGCGCCAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAACTAATATGATCAGCCAGTTTCGCTATCGGTGAATGGTAGCGATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGCATGGGTCACt                                              >  2:204700/1‑250 (MQ=255)
                                                                                                                                                                                                   aGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGTGTTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTAcc                                           <  1:204700/250‑1 (MQ=255)
                                                                                                                                                                                                           cAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCacaac                                   <  2:51311/250‑1 (MQ=255)
                                                                                                                                                                                                            aaTGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACAt                                 <  2:323640/251‑1 (MQ=255)
                                                                                                                                                                                                              tggtGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACAtct                               >  2:66980/4‑251 (MQ=255)
                                                                                                                                                                                                               gTTGAGCTGGGCGACAGAGACAAAAAACGCGTCCAGCACGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTCCGGCCGATTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACAtctc                              >  2:204102/1‑251 (MQ=255)
                                                                                                                                                                                                                  gAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGATGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCtccgt                           <  1:558061/251‑3 (MQ=255)
                                                                                                                                                                                                                  gAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAAt                                                                                                                                             <  1:318210/137‑1 (MQ=255)
                                                                                                                                                                                                                         gCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTac                    >  2:333511/1‑250 (MQ=255)
                                                                                                                                                                                                                          cgaaagagAAAAAAAAAGCGGCCAGTAAGTTTGTTTGTAATATTCTTGCCGAGGCATTACGACCTACTATCGGCGTTGCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTTAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTGTTACAGCTCTGGTTCGCCCGGAATGGGTCACTCCCAGCACACCATCTCCTTCCTGTaac                   <  2:573989/247‑1 (MQ=255)
                                                                                                                                                                                                                           gACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGGATGGGTCACTACCAGCACAACATCTCCTTCCTGTaac                   >  2:64444/1‑250 (MQ=255)
                                                                                                                                                                                                                                  aCAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGc           >  1:53803/1‑251 (MQ=255)
                                                                                                                                                                                                                                      aaaaaGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCgg                                                                                                                                                  >  2:273923/1‑112 (MQ=255)
                                                                                                                                                                                                                                      aaaaaGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCgg                                                                                                                                                  <  1:273923/112‑1 (MQ=255)
                                                                                                                                                                                                                                          aGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACACCGAAGCGGAcatc    >  2:72153/1‑250 (MQ=255)
                                                                                                                                                                                                                                           gCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGAcatcat  <  1:537251/251‑1 (MQ=255)
                                                                                                                                                                                                                                           gCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGAcatcat  <  1:143143/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
CGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCAT  >  NC_000913/4311970‑4312454

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