Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,362,225 G→A intergenic (‑291/+326) cadC ← / ← pheU DNA‑binding transcriptional activator CadC/tRNA‑Phe

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,362,2250GA97.9% 147.7 / ‑4.6 47intergenic (‑291/+326)cadC/pheUDNA‑binding transcriptional activator CadC/tRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (25/21);  minor base T (0/1);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

CAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTTTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGCTTTATGGATGCGTTACGCTCCTGATGACGTCATTTGACGTTCAACAGCATCACGGGGCCGCACGACATTTCACGTCAGTTAGTGCTATAGCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTT  >  NC_000913/4361980‑4362457
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                        
aaaCCTATTTTTACCATCTATAGATGGTTTTATTTACATATTTTTGGTGAATGCAAGACGTTATTTTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCAAGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGc                                                                                                                                                                                                                                     <  1:128366/250‑1 (MQ=255)
  aCCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTTGTTAATTCAAGACGTTATTTTTACCAGCCATAAACTTCTTTATTACATAGTATTACGAAAGTATTTTTCTGAGAACCAGAATTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGTAGCACTTTCTATTTTAATGGATTGTTAGTCTTTTCATTAACAAGCGa                                                                                                                                                                                                                                   <  1:331344/251‑1 (MQ=255)
      aTTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTTTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACt                                                                                                                                                                                                                                <  2:346925/250‑1 (MQ=255)
                  tataGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTTTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTctct                                                                                                                                                                                                                   >  1:578243/1‑251 (MQ=255)
                   ataGATGGGTTTATTTACATATTATTGGTGAATGCAAGCCGTTATTTTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCCTCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCtt                                                                                                                                                                                                                  >  1:318648/1‑251 (MQ=255)
                     aGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTTTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTg                                                                                                                                                                                                                >  1:264706/1‑251 (MQ=255)
                             ttatttaCATATTATTGGCGAATGCAAGACGTTATTTTTACCAGCCATAAACTTCTTGATTGCATAGTATTACGAAAGGATTTTACTGAGAATCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCAATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATCGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTTTTCTCAtt                                                                                                                                                                                                        >  2:185437/1‑251 (MQ=255)
                                   aCATATTATTGGTGAATGCAAGACGTTATTTTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATCTTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGaa                                                                                                                                                                                                   >  2:561310/1‑250 (MQ=255)
                                                             tATTTTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCgg                                                                                                                                                                        >  2:661140/1‑251 (MQ=255)
                                                             tATTTTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCgg                                                                                                                                                                        >  2:403152/1‑251 (MQ=255)
                                                                   tACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGtgttgt                                                                                                                                                                  >  2:554840/1‑251 (MQ=255)
                                                                    aCCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGtgttgt                                                                                                                                                                  <  2:615309/250‑1 (MQ=255)
                                                                     ccAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGc                                                                                                                                                                >  1:544376/1‑251 (MQ=255)
                                                                        gCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGCttt                                                                                                                                                             <  2:653446/251‑1 (MQ=255)
                                                                              aacttataaaTTACATAGAATAACGAAAGGATTTTACTGAAAACAAGAACTAATATTACTTACCATCAGAATTCATCATATTTACCAAAAAAAAACGTTCTGAAAATCACCCATGTAAACTTGTCCGCAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTATGTGCTCATCTAGAAAATCTCATTAGTCTCCATTCAGAACGGGGTGGTGCTTTATGGAt                                                                                                                                                       >  2:42155/10‑251 (MQ=255)
                                                                                      gATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATGTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGCTTTATGGATGCGTTACg                                                                                                                                               <  1:90624/251‑1 (MQ=255)
                                                                                         gaCATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGCTTTATGGATGCGTTACGCt                                                                                                                                             <  2:652458/249‑1 (MQ=255)
                                                                                                      gAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGCTTTATGGATGCGTTACGCTCCTGATGACGTCAt                                                                                                                               >  1:465700/1‑251 (MQ=255)
                                                                                                                      gaACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGcc                                                                                                                                                                                   >  1:236680/1‑183 (MQ=255)
                                                                                                                      gaACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGcc                                                                                                                                                                                   <  2:236680/183‑1 (MQ=255)
                                                                                                                                    atTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGCTTTATGGATGCGTTACGCTCCTGATGACGTCATTTGACGTTCAACAGCATCACGGGGCCGCAc                                                                                                 <  1:401712/251‑1 (MQ=255)
                                                                                                                                              aTCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGTAAGTTGTCCGGAACACTTTCTAGTTGAATGGATTGTTAGTCTTTGCATGATCAAGCGAACTGATGCATTTCTCTTTGTTCTCATTCATAAAATCTCATCAGTCGCCATTCAGGCCGGTGGTGTGCTTTAATGATGCGTTCCGCTCCTGAGGACTTCATTTGACGTTCAACAGCGTAACGGGGGCGCGCGGCTTTTCAc                                                                                       <  2:67079/251‑1 (MQ=255)
                                                                                                                                               tCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGCTTTATGGATGCGTTACGCTCCTGATGACGTCATTTGACGTTCAACAGCATCACGGTGCCGCACGACATTTCACg                                                                                      >  1:553387/1‑251 (MQ=255)
                                                                                                                                                  aaaTTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGCTTTATGGATGCGTTACGCTCCTGATGACGTCATTTGACGTTCAACAGCATCACGGGGCCGCACGACATTTCACGTcc                                                                                   >  2:524305/1‑250 (MQ=255)
                                                                                                                                                               ttGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGCTTTATGGATGCGTTACGCTCCTGATGACGTCATTTGACGTTCAACAGCATCACGGGGCCGCACGACATTTCACGTCAGTTAGTGCtata                                                                       <  1:554840/250‑1 (MQ=255)
                                                                                                                                                                 gCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGCTTTATGGATGCGTTACGCTCCTGATGACGTCATTTGACGTTCAACAGCATCACGGGGCCGCACGACATTTCACGTCAGTTAGTGCTATAGc                                                                     <  1:524305/250‑1 (MQ=255)
                                                                                                                                                                    aacgAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGCTTTATGGATGCGTTACGCTCCTGATGCCGTCATTTGACGTTCAACAGAATCACGGGGCCGCACGACATTTCACGTCAGTTAGTGCTATAGCTc                                                                   >  2:510990/4‑249 (MQ=255)
                                                                                                                                                                    aaaGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGCTTTATGGATGCGTTACGCTCCTGATGACGTCATTTGACGTTCAACAGCATCACGGGGCCGCACGACATTTCACGTCAGTTAGTGCTATAGCTc                                                                   >  1:5858/1‑249 (MQ=255)
                                                                                                                                                                          aaaTGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGCTTTATGGATGCGTTACGCTCCTGATGACGTCATTTGACGTTCAACAGCATCACGGGGCCGCACGACATTTCACGTCAGTTAGTGCTATAGCTCAGGAACaa                                                           <  1:443795/251‑1 (MQ=255)
                                                                                                                                                                          aaaTGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGCTTTATGGATGCGTTACGCGCCTGATGACGTCATTTGACGTTCGACATCATCACGGGGCCGTACGCCATTTCACGTCAGTTAGTGCTATAGCTCATGAACaa                                                           <  2:631588/251‑1 (MQ=255)
                                                                                                                                                                                  aGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGCGTTATGGATGCGTTACGCTCCTGATGACGTCATTTGACGTTCAACAGCATCACGGGGCCGCACGACATTTCACGTCAGTTAGTGCTATAGCTCAGGAACAAATTTTccc                                                   <  1:661140/251‑1 (MQ=255)
                                                                                                                                                                                   gAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGCTTTATGGATGCGTTACGCTCCTGATGACGTCATTTGACGTTCAACAGCATCACGGGGCCGCACGACATTTCCCGTCAGTTAGTGCTATAGCTCAGGAACAAATTTTCCCg                                                  >  2:488128/1‑251 (MQ=255)
                                                                                                                                                                                    aaaaTAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGCTTTATGGATGCGTTACGCTCCTGATGACGTCATTTGACGTTCAACAGCATCACGGGGCCGCACGACATTTCACGTCAGTTAGTGCTATAGCTCAGGAACAAATTTTCCCGa                                                 <  1:480099/251‑1 (MQ=255)
                                                                                                                                                                                       ataataCATGGAAAATTGTCAGGAACACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAACACTCATCAGTCGCCGCTCAGACCGGTGTTGTGCTTTATGGATGCGTTACGCTCCTGATGACGTCATTTGACGTTCAACAGCATCACGGGGCAGCACGACATTTCACGTCAGTAAGTGATATAGCTAAGGAACAAATTTTCCCGattt                                              >  2:651309/1‑248 (MQ=255)
                                                                                                                                                                                         aatCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGCTTTATGGATGCGTTACGCTCCTGATGACGTCATTTGACGTTCAACAGCATCACGGGGCCGCACGACATTTCACGTCAGTTAGTGCTATAGCTCAGGAACAAATTTTCCCGAATTg                                             <  2:217430/250‑1 (MQ=255)
                                                                                                                                                                                         aatCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGGTCAGACCGGTGTTGTGCTTTATGGATGCGTTACGCTCCTGATGACGTCATTTGACGTTCAACAGCATCACGGGGCCGCACGGCATTTCACGTCAGGTAGTGCTATAGCTCAGGAACAAAGTTTCCCGAATTgg                                            <  2:264762/251‑1 (MQ=255)
                                                                                                                                                                                                  aaaaTTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGGCTTTGCATGAACAGGCGAGCTGATACATTTCTCTTTTATCTCATTCAGAAAATATCTCCAGTCGCCGGTCGGACCGGTGCTGTGCTTTATGGATGCGTTACGCTCCTGATCACGTCATTTGACGTTCATTAGCATCACGTGGCCGCAAGGCATTTCACGGCAGTAAATGCTATAGCTGAGAAACAAATTTTCCCGAATTGGGATATGccc                                   >  1:610969/1‑251 (MQ=255)
                                                                                                                                                                                                   aaaTTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGCTTTATGGATGCGTTACGCTCCTGATGACGTCATTTGACGTTCAACAGCATCACGGGGCCGCACGACATTTCACGTCAGTTAGTGCTATAGCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCg                                  <  1:229122/251‑1 (MQ=255)
                                                                                                                                                                                                       tGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGCTTTATGGATGCGTTACGCTCCTGATGACGTCATTTGACGTTCAACAGCATCACGGGGCCGCACGACATTTCACGTCAGTTAGTGCTATAGCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCaa                               >  2:279370/1‑250 (MQ=255)
                                                                                                                                                                                                         tCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTTATTCAGAAAATCTCATCAGTCGCCGTTCAGAc                                                                                                                                                                           <  2:546305/108‑1 (MQ=255)
                                                                                                                                                                                                         tCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGAc                                                                                                                                                                           >  1:546305/1‑108 (MQ=255)
                                                                                                                                                                                                             gAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGCTTTATGGATGCGTTACGCTCCTGATGACGTCATTTGACGTTCAACAGCATCACGGGGCCGCACGACATTTCACGTCAGTTAGTGCTATAGCTCAGGAACaaa                                                          <  2:564250/217‑1 (MQ=255)
                                                                                                                                                                                                             gAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGCTTTATGGATGCGTTACGCTCCTGATGACGTCATTTGACGTTCAACAGCATCACGGGGCCGCACGACATTTCACGTCAGTTAGTGCTATAGCTCAGGAACaaa                                                          >  1:564250/1‑217 (MQ=255)
                                                                                                                                                                                                                 tcTTACTATTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACCTTTCTCTTTGTTCTCTTTCAGAAAATCTCATCAGTCGCCGTTCAGATCGGTGTTGTGCTTTATGGATGCGTTACGCTCCTGATGACGTCATTTGACGTTCAACAGCATCACGGGGCCGCACGACATTTCACGTCAGTTAGTTCTATAGCTCAGGCACAAATTTTCTCGAATTGGGATATGCCCGCAACTTGCTGGTgt                    >  2:512539/2‑250 (MQ=255)
                                                                                                                                                                                                                        aTTTTAATGGATTGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGCTTTATGGATGCGTTACGCTCCTGATGACGTCATTTGACGTTCAACAGCATCACGGGGCCGCACGACATTTCACGTCAGTTAGTGCTATAGCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTggg              <  1:408670/250‑1 (MQ=255)
                                                                                                                                                                                                                          tttAATGGATTGTTCGTCTTTGAATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGCTTTATGGATGCGTTACGCTCCTGATGACGTCATTTGACGTTCAACTGCATCACGGGGCCGCACGACATTTCACGTCAGTTAGTGCTATAGCTCAGGACCAAATTTTCCCGATTTGTGATATGCCCGCAAATTGCTGGAGCTGTgggg             >  2:95135/1‑249 (MQ=255)
                                                                                                                                                                                                                                   ttGTTAGTCTTTGCATGAACAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGCTTTATGGATGCGTTACGCTCCTGATGACGTCATTTGACGTTCAACAGCATCACGGGGCCGCACGACATTTCACGTCAGTTAGTGCTATAGCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGtt  >  1:465911/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                        
CAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATGCAAGACGTTATTTTTACCAGCCATAAACTTCTTGATTACATAGTATTACGAAAGGATTTTACTGAGAACCAGAAGTAATATTCCTTACCATCAAAATTCATCATCTTTGCCAAAGAAAAATGTTCAGAAAATAATCCATGGAAAATTGTCCGGAGCACTTACTATTTTAATGGATTGTTAGTCTTTGCATGAGCAAGCGAACTGATACATTTCTCTTTGTTCTCATTCAGAAAATCTCATCAGTCGCCGTTCAGACCGGTGTTGTGCTTTATGGATGCGTTACGCTCCTGATGACGTCATTTGACGTTCAACAGCATCACGGGGCCGCACGACATTTCACGTCAGTTAGTGCTATAGCTCAGGAACAAATTTTCCCGAATTGGGATATGCCCGCAAATTGCTGGTGATGTGGGGAGAATCTGGTT  >  NC_000913/4361980‑4362457

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