Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,838,265 G→A P334S (CCT→TCT)  ascG ← DNA‑binding transcriptional repressor AscG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,838,2650GA97.3% 105.0 / NA 37P334S (CCT→TCT) ascGDNA‑binding transcriptional repressor AscG
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (25/11);  minor base C (1/0);  total (26/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGA‑ATTTTTACGCTGGAGAGCGCCGGAACGGTGTA‑GGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGCACCGCTACGCCGCGCTCGTGTAACGCTTTCATCGCACCTATCGC  >  NC_000913/2838042‑2838511
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                        
acgcTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAg                                                                                                                                                                                                                               <  1:138448/251‑1 (MQ=255)
          acaGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTCCTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTACGCTg                                                                                                                                                                                                                     >  1:375626/1‑251 (MQ=255)
                    cAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTcg                                                                                                                                                                                                           >  2:185376/1‑251 (MQ=255)
                                 tacgcgcgTCAGCAATGATGATACGGTTCATCAAATTCTCCAGCAATGACATGTAATGCGCCGATACATTCACAAATCATTCCAGTTTATAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCTGGAAAAACGTTGTCATCGACCCTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGAAGCAAGGAGTGAGTCGCGGCGGAGCAGTTTTCCGCTGAAGGTTTTCGTCGGTGATAAAt                                                                                                                                                                                               <  2:418662/246‑1 (MQ=255)
                                           aGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATcccc                                                                                                                                                                                           >  1:338005/1‑244 (MQ=255)
                                           aGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATcccc                                                                                                                                                                                           <  2:338005/244‑1 (MQ=255)
                                                 gatgaAAAGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGCTGACATGTGATGACAATGTTTATCGCGAAGAAGCACTGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGCTTTCGGCGGTGAGaaaatcaccacattcaccat                                                                                                                                                                               >  2:137354/1‑233 (MQ=255)
                                                                caAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAAGTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTcc                                                                                                                                                                <  2:74561/250‑1 (MQ=255)
                                                                 aaaTTCTCCAGCAAGGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTTTTTAAGGAAATTAAGTTCTGTAATGCAGAAAAAACGATGTCAGCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCaa                                                                                                                                                              <  2:399847/251‑1 (MQ=255)
                                                                                     ttAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCAtctc                                                                                                                                           >  1:146535/1‑250 (MQ=255)
                                                                                       aaTGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATACCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTTCTGAATCATATCAgc                                                                                                                                        >  2:9080/1‑250 (MQ=255)
                                                                                            gcCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCg                                                                                                                                   >  1:128314/1‑251 (MQ=255)
                                                                                              cGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGTGCATAAAAATCATACGTCCCATATTTTCCTGAAACATCTCAGTcccccg                                                                                                                                  >  2:230673/1‑245 (MQ=255)
                                                                                                  aCATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCCTCTCAGTCACCGGA‑Attt                                                                                                                            >  2:332499/1‑251 (MQ=255)
                                                                                                                tGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGGTTTCGGCGGTGAGAAATCCCCCCCATCGAGCATAAAAATCAGACGGCCAATAATTTCCTGAATCATCTCACACAACGGACATTTTTAC‑CCGGAGAgc                                                                                                              >  2:63586/1‑251 (MQ=255)
                                                                                                                       ttttAATTTACTGTTCTTTAAGGAAATTAATTTCTGTAATGCTGGAAAAACGATGTCATCGACACTAGTGACGATGACCTGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCAGCTGAAGGTTTTCGGAGGTGCGAAAACCCCGCCATCGAGGATAGAAAACCGCAGGCCAATAATTTCCTGCATCATCTCAATCCACGGA‑Attttttcgttagagtcggacagc                                                                                                        >  1:384722/1‑232 (MQ=255)
                                                                                                                         ttAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGA‑ATTTTTACGCTGGAGAGCGCCGGAAc                                                                                                      <  2:324094/250‑1 (MQ=255)
                                                                                                                                ccTGTTCTTTAAGGAAATTAATTTCTGTAATGCCGGAAAAACGATGTCATCGACACTAGTGACGATGCCATGTGATGACAATGTTTATCGCGACGCAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGGGAAATCCCCGCCATCGAGCATAAAAATCAGCCCTCCCATAATTTCCTGAATCATCCCCGTCACCAGA‑AATTTTACGCTGGAGAGCGCCGGGACGgtgttgg                                                                                              >  2:60783/2‑248 (MQ=255)
                                                                                                                                   gTTATTTAAGGAAATTAATTTCTTTAATGCAGGAAAAACGATGTCATCGACACTATGGACGATGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGGCCAATAATGTCCTGAATCATCTCAGTCACCTGA‑ATTTTTACGCTGGAGAGTGCCGGAACTGTGTA‑gggg                                                                                           <  2:273538/250‑1 (MQ=255)
                                                                                                                                    ttATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCCTCGACAATCGTGACGATGACATGTGATGACAATGTTTATCGCGCAGAAACAATGAGTGCGTCGCGGAGGATCAGTTTTCCGCTTAAGGTTTTCGGCGGTGAGAAATCCCCGACATCGAACATACAAATCACCCGCCCAATCATTTCCTGAAACAACTCAGTCAACGGA‑ATTTTTTCGCTTGAGATAGCCGGCACGAAGTA‑GGGGGc                                                                                         >  1:20564/1‑251 (MQ=255)
                                                                                                                                    ttATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGACGTCCAATAATTTCCTGAATCATCTCAGTCACCGAA‑ATTTTTACGCTGGAGAGCGCCGGAACGGTGTA‑GGGGGc                                                                                         >  2:31358/1‑251 (MQ=255)
                                                                                                                                       tttAAGGAAATTAATTTCTGTAATGCAGAAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGGTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATTAGCCGTCCAATAATTTTCTGAATAATATCAGACACCCGA‑ATTTTTACGCTGGAGATCGCCGGAACGGCGTA‑GGGGGCAAt                                                                                      >  2:417356/1‑251 (MQ=255)
                                                                                                                                                                             gTCATCGACACTAGTGACGCTGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGA‑ATTTTTACGCTGGGGAGCGCCGGAACGGTGTA‑GGGGGCAATAGCGATATCATCGAATCCGGTAACTGACACCTTCTCTg                                                >  1:328527/1‑251 (MQ=255)
                                                                                                                                                                               cATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTGTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGA‑ATTTTTACGCTGGAGAGCGCCGGAACGGTGTA‑GGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGc                                              <  1:294955/251‑1 (MQ=255)
                                                                                                                                                                                   gACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGA‑ATTTTTACGCTGGAGAGCGCCGGAACGGTGTA‑GGGGGCCATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGCAccc                                          >  1:201671/1‑250 (MQ=255)
                                                                                                                                                                                   gACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGA‑ATTTTTACGCTGGAGAGCGCCGGAACGGTGTA‑GGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGCACCg                                          >  1:27085/1‑251 (MQ=255)
                                                                                                                                                                                   gACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGAGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCCACGGA‑ATTTTTACGCTGGAGAGCGCCGGAACGGTGTA‑GGGGGCAATAGCGATATCATCGCATCCGCTAACTGACACCTGCTCTTGCACCg                                          >  1:191697/1‑251 (MQ=255)
                                                                                                                                                                                    acacTAGTGACGATGACATGTGATGACAATGTTTATCGCGCAGAAGCAATGAGTGAGTCGCGGCGGATCCGTTTTCCGCTGAAGGTTTTCGGCGGGGGGAAATCCCCGCCATCGAGCATAAAAATAAACCGTCCAATAATTTCCTGAATCATCTCAGTCCCCGGA‑ATTTTTAAGCTTGAGAGCGCCGGAACGGGGTA‑GGGGGCAAAAGCGATATCATCGAATCCGATAAATGAAACCTGCTCTTGCACCGc                                         >  2:26876/1‑251 (MQ=255)
                                                                                                                                                                                    acacTAGTGACGAGGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTGCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATGTCCTGAATCATCTCAGTCACCGGA‑ATTTTTACGCTGGAGAGCGCCGGAACGGTGTA‑TGGGGCAAGAGCGATATCATCGAATCCGATAACTGACACCGGCTCTGGCACCGc                                         <  2:390023/251‑1 (MQ=255)
                                                                                                                                                                                     cacTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGAAGAAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCAACGGGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGACACCGGA‑AATTTTACGCTGGAGCGCGCCGGAACGGTGTC‑AGGGGCAATAGCGGTATAATCGAATACGATAAATGACACCTGCTCTGGCACCGCt                                        >  1:277532/1‑251 (MQ=255)
                                                                                                                                                                                       cTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTACAATAATTTCCTGAATCATCTCAGTCACCGGG‑ATTTTTACGCTGGAGAGCGCCGGGACGGTGTA‑TGGGGCAATAGCGATATCATCGAATCCGAGAACTGACACCTGCTCTGGCACCGCTAc                                      >  2:33219/1‑251 (MQ=255)
                                                                                                                                                                                       cTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCAGCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGA‑ATTTTTACGCTGGAGAGCGCCGGAACGGTGTA‑gg                                                                                             <  1:33714/196‑1 (MQ=255)
                                                                                                                                                                                        tAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGGTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGA‑ATTTTTACGCTGGAGAGCGCCGGAACGGGGTA‑GGGGGCAATAGCGATATCATCGAATCAGAAAACTGACACCTGCTCTGGCACCGCTAcc                                     >  1:87567/1‑250 (MQ=255)
                                                                                                                                                                                               gATGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGA‑ATTTTTACGCTGGAGAGCGCCGGAACGGTGTA‑GGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTg                                                <  2:232239/233‑1 (MQ=255)
                                                                                                                                                                                               gATGACATGTGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGA‑ATTTTTACGCTGGAGAGCGCCGGAACGGTGTA‑GGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTg                                                >  1:232239/1‑233 (MQ=255)
                                                                                                                                                                                               gATGACATGAGATGACAATGTTTATCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGA‑ATTTTTACGCTGGAGAGCGCCGGAACGGTGTA‑GGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGCACCGCTACGCCgcgcc                              >  2:136274/1‑250 (MQ=255)
                                                                                                                                                                                                        tgatgaCAATGTTTAGCGCGAAGAAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGATGACGGGTTTAGGAGGAGAGAAATCCATGCCATCGAGTATAACAAAGGGCCGTCAAAGAAGTTCCTCAAGCCCCTCAATAGCCCGA‑ATTTTTGCGCTGGGAAACGACGGGACGG‑GTAAGGTGGCTAGAGGGATATCATCGAATCCAATAACTAGCACCTGCACTGGTACCACTAcgcccagccctcgttaa                      >  2:109926/1‑237 (MQ=255)
                                                                                                                                                                                                                           gaagaaGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGGGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGC‑ATTTTTACGCTGGAGCGCGCCCGAACGGTGTA‑GGGGGCAATAGCGATATCATCGAATCCGATTACTGACACCTGCTCTGGCACCGCTACGCCGCGCTCGTGTAACGCTTTCATCGCACCTATCGc  >  2:21664/1‑251 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                        
ACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGA‑ATTTTTACGCTGGAGAGCGCCGGAACGGTGTA‑GGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGCACCGCTACGCCGCGCTCGTGTAACGCTTTCATCGCACCTATCGC  >  NC_000913/2838042‑2838511

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