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→ACC) alsRDNA‑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

CGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAAGTTAGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCAT  >  NC_000913/4311970‑4312454
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
attttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑Gcaaggttggt                                                                                                                                                                                                                                             >  2:537251‑M2/88‑248 (MQ=255)
 ttttatgtggggtccttccctgatattacgggcactatttattcacaactctgccgaaaaacaggctgtcgtcagttttgacgtgaGGTAATGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGAG‑CCTGGTGAGAAGAAATCAACCATTGCGCCGGTTTTTTCCATATTAAGATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAggttggtt                                                                                                                                                                                                                                            >  2:166302‑M2/87‑251 (MQ=255)
             tcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTg                                                                                                                                                                                                                               >  2:171240‑M2/75‑251 (MQ=255)
              cattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGc                                                                                                                                                                                                                              <  1:89801‑M2/178‑1 (MQ=255)
                          ttacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACg                                                                                                                                                                                                                  >  2:111305‑M2/62‑251 (MQ=255)
                                     atttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACg                                                                                                                                                                                                       <  2:295170‑M2/201‑1 (MQ=255)
                                                 ctctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCgggg                                                                                                                                                                                          >  1:285417‑M2/39‑251 (MQ=255)
                                                     gacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTGCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTg                                                                                                                                                                                       <  1:430851‑M2/216‑1 (MQ=255)
                                                            aacacgctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATCTTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAgct                                                                                                                                                                                >  2:544250‑M2/28‑249 (MQ=255)
                                                                ggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAaataa                                                                                                                                                                           <  2:6429‑M2/228‑1 (MQ=255)
                                                                    gtcgtcagttttgacgtgaGGTAAGGGTACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGATGCATTACGACCTAATAACGG‑CGTTTCCGGTGCTGGTGAGCAaataatatg                                                                                                                                                                       <  2:549969‑M2/232‑2 (MQ=255)
                                                                        tcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCg                                                                                                                                                                   >  1:185317‑M2/16‑251 (MQ=255)
                                                                                            ggggACAG‑TAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAatat                                                                                                                                                                                                                                    >  2:90275/1‑167 (MQ=255)
                                                                                            ggggACAG‑TAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAatat                                                                                                                                                                                                                                    <  1:90275/167‑1 (MQ=255)
                                                                                             gggACAG‑TAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGa                                                                                                                                               >  1:381647/1‑250 (MQ=255)
                                                                                                    gTAACTCCGCGATTGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGACGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAg                                                                                                                                        <  2:438450/251‑1 (MQ=255)
                                                                                                    gTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGt                                                                                                                                                                                     >  2:311484/1‑206 (MQ=255)
                                                                                                    gTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGt                                                                                                                                                                                     <  1:311484/206‑1 (MQ=255)
                                                                                                            gcgATGGCGGAATTTCCTTTTTCTTATTGCAGCGCG‑CCAGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCAGGGCGACAGAGACGAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTGCGCTGTCGGTGAATGGTAGCTATGGGt                                                                                                                                <  2:204691/251‑1 (MQ=255)
                                                                                                                           ccTTTTTCTTATTTTAGCGCGGCCTG‑TGAGGAGGAATCCACAATTTCGCCGGTTTTTTCCATATTAATAATAGCTTGTTCAATGTTGAGCTGGGGGACAGGAACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTAGGGGTTATACAAATAATCt                                                                                                                  <  2:518494/250‑1 (MQ=255)
                                                                                                                                    tATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCTGGTTTTTTCCATATTACTATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAACAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGACGAGGCATTAAGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGCTCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGtt                                                                                                        >  2:83762/1‑251 (MQ=255)
                                                                                                                                          aGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCttttt                                                                                                  <  2:12206/251‑1 (MQ=255)
                                                                                                                                            cgcg‑cCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTAAGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTg                                                                                                >  2:52043/1‑251 (MQ=255)
                                                                                                                                                            gaaATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGAAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCg                                                                                 >  1:259513/1‑251 (MQ=255)
                                                                                                                                                                tCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCttt                                                                             >  2:388671/1‑251 (MQ=255)
                                                                                                                                                                  aacaaTTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCtttt                                                                            >  1:217614/1‑250 (MQ=255)
                                                                                                                                                                     aatggggcCGGTTTTTTCCATATAAATATTAGCTGGTTCAAGGTTGAGCTGGGCGAAAGAGACAAAAAAAGC‑TCCAAGAAAGGTTGGTT‑GTAATATCCTGGCGGAGGAATTCCGCCCTAA‑AACGGGCGTTTCCGGGGCTGGTGGGCAAAAAATATGACGGGCCGGTTCGGCTACCGTTGAATGGTAGCAATGGGT‑ATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATc                                                                       <  2:185317/245‑1 (MQ=255)
                                                                                                                                                                           gcCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTgg                                                                  >  1:293748/1‑251 (MQ=255)
                                                                                                                                                                                        atattaatattaGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATg                                                     >  2:143143/1‑251 (MQ=255)
                                                                                                                                                                                        atattaatattaGCTTGATCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCGATGGGTTATACAAATAATCGTTGCCCGGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATg                                                     >  1:567889/1‑251 (MQ=255)
                                                                                                                                                                                                attaGCTTGTTCAATGTTGAGCTGGGCGCCAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAACTAATATGATCAGCCAGTTTCGCTATCGGTGAATGGTAGCGATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGCATGGGTCACt                                              >  2:204700/1‑250 (MQ=255)
                                                                                                                                                                                                   aGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGTGTTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTAcc                                           <  1:204700/250‑1 (MQ=255)
                                                                                                                                                                                                           cAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCacaac                                   <  2:51311/250‑1 (MQ=255)
                                                                                                                                                                                                            aaTGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACAt                                 <  2:323640/251‑1 (MQ=255)
                                                                                                                                                                                                              tggtGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACAtct                               >  2:66980/4‑251 (MQ=255)
                                                                                                                                                                                                               gTTGAGCTGGGCGACAGAGACAAAAAACGCGTCCA‑GCACGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTCCGGCCGATTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACAtctc                              >  2:204102/1‑251 (MQ=255)
                                                                                                                                                                                                                  gAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGATGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCtccgt                           <  1:558061/251‑3 (MQ=255)
                                                                                                                                                                                                                  gAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAAt                                                                                                                                             <  1:318210/137‑1 (MQ=255)
                                                                                                                                                                                                                         gCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTac                    >  2:333511/1‑250 (MQ=255)
                                                                                                                                                                                                                          cgaaagagAAAAAAAAAGCGGCCA‑G‑TAAGTTTGTTTGTAATATTCTTGCCGAGGCATTACGACCTACTATCGG‑CGTTGCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTTAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTGTTACAGCTCTGGTTCGCCCGGAATGGGTCACTCCCAGCACACCATCTCCTTCCTGTaac                   <  2:573989/247‑1 (MQ=255)
                                                                                                                                                                                                                           gACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGGATGGGTCACTACCAGCACAACATCTCCTTCCTGTaac                   >  2:64444/1‑250 (MQ=255)
                                                                                                                                                                                                                                  aCAAAAAAAGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGc           >  1:53803/1‑251 (MQ=255)
                                                                                                                                                                                                                                      aaaaaGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCgg                                                                                                                                                  >  2:273923/1‑112 (MQ=255)
                                                                                                                                                                                                                                      aaaaaGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCgg                                                                                                                                                  <  1:273923/112‑1 (MQ=255)
                                                                                                                                                                                                                                          aGCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACACCGAAGCGGAcatc    >  2:72153/1‑250 (MQ=255)
                                                                                                                                                                                                                                           gCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGAcatcat  <  1:537251/251‑1 (MQ=255)
                                                                                                                                                                                                                                           gCGTCCA‑GCAAGGTTGGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGAcatcat  <  1:143143/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
CGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAAGTTAGTT‑GTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCAT  >  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: ‑