Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 419,294 A→G intergenic (+110/‑297) phoR → / → brnQ sensor histidine kinase PhoR/branched chain amino acid transporter BrnQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913419,2940AG97.4% 129.0 / ‑3.9 39intergenic (+110/‑297)phoR/brnQsensor histidine kinase PhoR/branched chain amino acid transporter BrnQ
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (16/22);  total (16/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

GGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTGCTTTAAAACGTTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGAAAAACTATTCTTCGCCGCGCCTGGTTGGGAGTATTTCCCGCTAAAATTGTTTAAATATACCGCTGTATCATCCCCAGG  >  NC_000913/419050‑419537
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
ggATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGAtt                                                                                                                                                                                                                                                >  1:107191/1‑250 (MQ=255)
 gATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGAtt                                                                                                                                                                                                                                                <  1:324618/249‑1 (MQ=255)
          tAGCCATCGTGAATCATGCAGTGAAATAGCACGAAAGTCGGCTGTATATTGAGAGTAAAGTAGGAAAAAGAACACGGTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTTAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCt                                                                                                                                                                                                                                      <  2:158947/250‑1 (MQ=255)
                    gAAACATGCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGaa                                                                                                                                                                                                                           <  1:530219/251‑1 (MQ=255)
                        cATGCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAg                                                                                                                                                                                                                       <  1:170547/251‑1 (MQ=255)
                              gtgAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGTGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGtttt                                                                                                                                                                                                                 <  1:515991/251‑1 (MQ=255)
                                    catcaCGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACtg                                                                                                                                                                                                           <  2:107191/251‑1 (MQ=255)
                                    catcaCGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACt                                                                                                                                                                                                            <  1:579623/250‑1 (MQ=255)
                                        aCGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACtggttg                                                                                                                                                                                                       <  1:396932/251‑1 (MQ=255)
                                          gAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACtggttggt                                                                                                                                                                                                     >  1:92993/1‑251 (MQ=255)
                                                                gTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACct                                                                                                                                                                               <  1:15709/251‑1 (MQ=255)
                                                                 tACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACctc                                                                                                                                                                              >  1:180905/1‑251 (MQ=255)
                                                                     gTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGt                                                                                                                                                                          <  1:95756/251‑1 (MQ=255)
                                                                           aaaGGAACAAGATTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGCTTTTCAAACGATTGGGGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTCTCGCTATATACCTCTGGttttt                                                                                                                                                                      >  2:85793/1‑249 (MQ=255)
                                                                               gAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCTATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTAACGGAAATAGTGTTTTATACTGGTTGTTGATTTCTTATCTCTATATACCTCTGGTTTTTATATcc                                                                                                                                                                >  2:158510/1‑251 (MQ=255)
                                                                                       ttCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTg                                                                                                                                                         <  2:215586/250‑1 (MQ=255)
                                                                                               tGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTGCTTTaaaa                                                                                                                                                 <  1:476302/250‑1 (MQ=255)
                                                                                                               ttAATTGCCAAAAACCGCGATTAATCCGCCTTTGTCATCTTTTATTTGCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTGCTTTAAAACGTTATAAGCGTTTaa                                                                                                                                 <  2:222436/250‑1 (MQ=255)
                                                                                                               ttAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGAGGCTGGCTTATTGTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGGGATTAGCCATGTCTTTTTCACGGAAATAGTGTGTTATACTGGTTGGTGATTTCTTATCGCTATATACCCCTGGTTTGTAGATCCCTCCTTGCTTTAAAACGTTATAAGCGTTTaa                                                                                                                                 <  2:534896/250‑1 (MQ=255)
                                                                                                                  aTTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTGCTTTAAAACGTTATAAGCGTTTAAATTg                                                                                                                             >  1:206786/1‑251 (MQ=255)
                                                                                                                       cAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTGCTTTAAAACGTTATAAGCGTTTAAATTGCGCtt                                                                                                                        >  2:429649/1‑251 (MQ=255)
                                                                                                                             cAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTTGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTGCTTTAAAACGTTATAAGCGTTTAAATTGCGCTTCAGGTg                                                                                                                  <  2:415250/251‑1 (MQ=255)
                                                                                                                                              ttggcATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTGCTTTAAAACGTTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATactgactg                                                                                                   <  2:328504/245‑1 (MQ=255)
                                                                                                                                                   atattttATTGCCATAAGCCAGTCGATGCGTGATTATTTTTTTTTAAGTCAAAATACGGGCGTTAGATTGTTCAACGAGTTGTTATTTTTTTTGCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACATCTGGGTGTTAGATCCCTCCTTGCTTTAAAACGTTATAAGCGTGTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCAtttac                                                                                            <  2:43674/248‑4 (MQ=255)
                                                                                                                                                             gCCATAAGCCAGTCGATGCTTGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGAGTTTACAACGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTGCTTTAAAACGTTATAAGCGGTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAAt                                                                                   <  2:446325/250‑1 (MQ=255)
                                                                                                                                                               cATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTGCTTTAAAACGTTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCg                                                                                 >  2:27095/1‑250 (MQ=255)
                                                                                                                                                                aTAAGCCCGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTGCTTTAAAACGTTATAAGCGTTGAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGaa                                                                               <  2:560300/251‑1 (MQ=255)
                                                                                                                                                                     ccAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTGCTTTAAAACGTTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGAAAAACt                                                                          >  1:120353/1‑251 (MQ=255)
                                                                                                                                                                       cgTCGATGCAGGCTTATTTTCTTTGAAGTCAAAATACGGGCATTCGATTTTACAACGATTGGAGATTTTTTGTTCGCGTGATTAGACATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTGCTTTAAACCGTTATAAGCGTTTACATTGCGCATCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGAACAACTAt                                                                        >  2:242658/2‑251 (MQ=255)
                                                                                                                                                                           gATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTGCTTTAAAACGTTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATGAACGCGGTAAATCGAAAAACTAttctt                                                                    >  2:389927/1‑251 (MQ=255)
                                                                                                                                                                                              tGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATTCCTCTGGTTTTTAGATCCCTCCTTGCTTTAAAACGTTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCAGTAACGCGGTAAATCGAAAAACTATTCTGCGCCGCGCCTGGTTGGGa                                                  <  2:69194/250‑1 (MQ=255)
                                                                                                                                                                                               gCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTGCTTTAAAACGTTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGAAAAACTATTCTTCGCCGCGCCTGGTTGGGAGt                                                <  1:514446/251‑1 (MQ=255)
                                                                                                                                                                                                      aaaTACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTGCTTTAAAACGTTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGAAAAACTATTCTTCGCCGCGCCTGGTTGGGAGTATTTccc                                         >  2:553941/1‑251 (MQ=255)
                                                                                                                                                                                                                            aaCGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTGCTTTAAAACGTTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGAAAAACTATTCTTCGCCGCGCCTGGTTGGGAGTATTTCCCGCTAAAATTGTTTAAATATAcc                   >  2:246391/1‑251 (MQ=255)
                                                                                                                                                                                                                            aaCGATTGGTGAGTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTGCTTTAAAACGTTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGAAAAACTATTCTTCGCCGCGCCTGGTTGGGAGTATTTCCCGCTAAAATTGTTTAAATATAcc                   <  1:76259/251‑1 (MQ=255)
                                                                                                                                                                                                                             aCGATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTGCTTTAAAACGTTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGAAAAACTATTCTTCGCCGCGCCTGGTTGGGAGTATTTCCCGCTAAAATTGTTTAAATATACCg                  >  2:463717/1‑251 (MQ=255)
                                                                                                                                                                                                                               gATTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTGCTTTAAAACGTTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGAAAAACTATTCTTCGCCGCGCCTGGTTGGGAGTATTTCCCGCTAAAATTGTTTAAATATACCGCt                >  1:420013/1‑251 (MQ=255)
                                                                                                                                                                                                                                aTTGGTGATTTTTTGTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTGCTTTAAAACGTTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGAAAAACTATTCTTCGCCGCGCCTGGTTGGGAGTATTTCCCGCTAAAATTGTTTAAATATACCGCTg               >  1:400036/1‑251 (MQ=255)
                                                                                                                                                                                                                                              gTTCGCGTGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTGCTTTAAAACGTTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGAAAAACTATTCTTCGCCGCGCCTGGTTGGGAGTATTTCCCGCTAAAATTGTTTAAATATACCGCTGTATCATCCCCAgg  <  1:30150/250‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
GGATTAGGGTTAGCGATCGTGAAACATGCTGTGAATCATCACGAAAGTCGCCTGAATATTGAGAGTACAGTAGGAAAAGGAACACGTTTCAGTTTTGTTATCCCGGAACGTTTAATTGCCAAAAACAGCGATTAATCCGCCTTTGTCATCTTTTATTGCCATAAGCCAGTCGATGCTGGCTTATTTTCTTTGCAGTCAAAATACGGGCGTTAGATTTTACAACGATTGGTGATTTTTTGTTCGCATGATTAGCCATGTCTTTTTCACGGAAATAGTGTTTTATACTGGTTGGTGATTTCTTATCGCTATATACCTCTGGTTTTTAGATCCCTCCTTGCTTTAAAACGTTATAAGCGTTTAAATTGCGCTTCAGGTGCTGTCATACTGACTGCATTAACGCGGTAAATCGAAAAACTATTCTTCGCCGCGCCTGGTTGGGAGTATTTCCCGCTAAAATTGTTTAAATATACCGCTGTATCATCCCCAGG  >  NC_000913/419050‑419537

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