Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,647,949 C→T intergenic (‑161/‑114) dinQ ← / → agrA membrane toxin DinQ/small RNA AgrA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,647,9490CT97.6% 121.9 / ‑4.5 42intergenic (‑161/‑114)dinQ/agrAmembrane toxin DinQ/small RNA AgrA
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (17/23);  total (17/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

GAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGA  >  NC_000913/3647725‑3648188
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
aaaaGCGGATCAGTTCCAGCAGCGCAATTAACGCACCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATTCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCag                                                                                                                                                                                                                       <  2:571267/250‑1 (MQ=255)
      ggATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCCTGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTa                                                                                                                                                                                                                  >  1:48503/1‑250 (MQ=255)
              accagcagctacaTTACCGCCCCTCGACCGATTATTTCTTTATCAAGCGCTCGGTCTCTCCGTGCGATGGAGTGAGAATGAGTCCGCTTACTCATCCACTGCCTGTCACGGCGCATTTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCAAGAGCATCAGTTTACGTCACCATCCGTACTTCTGACCAGGACATTGAAATCGTTTATGTTCATTGCCAGAGTGAAAATGCCTCAGAGCTAGCTGGATa                                                                                                                                                                                                          <  2:17888/238‑1 (MQ=255)
                     gcgcAATTAACGCCACTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACa                                                                                                                                                                                                  <  2:240393/251‑1 (MQ=255)
                      cgcAATTATCGTCCCTAGATCGATGATTTCTTTATCAATCACCCGTTTGCTCCATGCGATGGAGTGCGAATGCGTCCGCTTACTCTTCCACTGCCTGACACTGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGTCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATTTGCATTGCTATAGTGAAAATGCCTTAGAGCTAACTGGATAATCATACAg                                                                                                                                                                                                 <  2:427990/251‑1 (MQ=255)
                              aCGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCa                                                                                                                                                                                         <  1:479137/251‑1 (MQ=255)
                                gCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGCGGGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGCGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGGGAAATTGCCTCAGAGCTAACTGGATAATAATACAATACATGAAgg                                                                                                                                                                                       >  2:358311/1‑251 (MQ=255)
                                gCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAgg                                                                                                                                                                                       >  1:120911/1‑251 (MQ=255)
                                  cccTAGAACGATGCTTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGtt                                                                                                                                                                                     >  2:578856/1‑251 (MQ=255)
                                  cccTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGCGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATAATACAGTACATGCAGGtt                                                                                                                                                                                     >  2:405127/1‑251 (MQ=255)
                                   ccTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTta                                                                                                                                                                                    <  1:102878/251‑1 (MQ=255)
                                           gatgatTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGCGCTAACTGGATAATCATACAGTACATGCAGGTTATaaaccc                                                                                                                                                                             >  2:254096/1‑247 (MQ=255)
                                              gatTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAg                                                                                                                                                                                                 >  1:200853/1‑227 (MQ=255)
                                              gatTGCTTTATCAATCACCCGTTTTCTCCAGGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACGCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCAGTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAg                                                                                                                                                                                                 <  2:200853/227‑1 (MQ=255)
                                               atTGTTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGTGAATTCAACCGCTTACTCATCCACTGCCTTGCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCCTCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGTATAATCATACAGTACATGCTGGTTATCAACCCAGCAc                                                                                                                                                                        <  2:110674/251‑1 (MQ=255)
                                                  gCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTc                                                                                                                                                                     >  1:527198/1‑251 (MQ=255)
                                                     ttATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCt                                                                                                                                                                   <  1:557937/250‑1 (MQ=255)
                                                          aaccccctggtgtctcCATGCGATGGATGGAGAATGATTCCGCATACTCATCGAATGCCTGTCACGGGGCATGACTCAGTGGTAGATAAGAGTTCAATCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTTACAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCAGAGAGTACAGGCAGGTTATAAAATCAGCATGTCCTTGCAAt                                                                                                                                                             <  2:444857/239‑1 (MQ=255)
                                                             caCCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGTATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGt                                                                                                                                                          <  1:297781/251‑1 (MQ=255)
                                                              ccccGGTTTCTCAATGCGAAGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCGGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTCCTTCTGACCAGGACTTTGAAAGCGTTTATGCGCATTGCTAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGtt                                                                                                                                                         <  2:147521/250‑1 (MQ=255)
                                                                      tctcCATTCGAGGGAGTGAGAATGCATCCGACTTCTCATCCTCGTCCTGTAACGTCGCATGTCGAATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACAGGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTAt                                                                                                                                                  <  2:128802/250‑1 (MQ=255)
                                                                              cGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGc                                                                                                                                         >  1:362834/1‑251 (MQ=255)
                                                                                       gagaATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCg                                                                                                                                <  1:625245/251‑1 (MQ=255)
                                                                                         gaATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGtt                                                                                                                              <  1:661315/251‑1 (MQ=255)
                                                                                               aTCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGAGGCGTTAgatg                                                                                                                         >  2:430800/1‑250 (MQ=255)
                                                                                                             tccacttctggTCACGGCGCATTTCGCATTGTTAGATAAGAACTCTATCACTCCGGCCTTAGCATCAGTGAACGGCACCACACGTACTTCTTACCAGGACATTGAAATCGTTCATGGGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATATATGACATGTATGTTATAAAACCAGAACTTCCTTGCAATAGTTTCAGTATGTTATTAGCATTGATGCGTTAGATGATGGCTAtctc                                                                                                              <  2:416015/238‑1 (MQ=255)
                                                                                                                            ggCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGccacc                                                                                            >  1:237211/1‑251 (MQ=255)
                                                                                                                                        aTTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTg                                                                                >  1:214262/1‑251 (MQ=255)
                                                                                                                                                           ctcACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGAc                                                               <  1:31207/249‑1 (MQ=255)
                                                                                                                                                                 ccGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCaga                                                                                                  >  2:474858/1‑208 (MQ=255)
                                                                                                                                                                 ccGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCaga                                                                                                  <  1:474858/208‑1 (MQ=255)
                                                                                                                                                                          gcctcatttAACGGCCTCACCCGTCCCTTTGACCAGGACTTTGAAAGCGTTTATGTGCATTGGCCGAGTGAATATGCCTCAGAGCTATCTGGTTAAGCATACAGTACATGCAGGTTTTAAAACCAGCCCGGCCTTGCAATAGTTTCAGTATGGGATTAGCATTGATGCTTTAGATGATTGCTATCTCTCTCCAGTCAGAGTGACCAACGCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTTGGtttttta                                              <  2:368930/244‑2 (MQ=255)
                                                                                                                                                                                    aCGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGAt                                    <  1:405127/251‑1 (MQ=255)
                                                                                                                                                                                        caccacCCGTCCTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACATTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCgg                                <  1:80659/251‑1 (MQ=255)
                                                                                                                                                                                               cGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCATAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGCTGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGATGGAAAGTAAACAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGc                         >  2:323052/1‑251 (MQ=255)
                                                                                                                                                                                                gTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCt                        >  2:524062/1‑251 (MQ=255)
                                                                                                                                                                                                     tCTGACCAGGACTTTGAAAGCGGTTATGTGCAGTGCCAGAGGTAAAATGCCTCAGATCTAACTGGATAATCATACAGTACATGCATGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGACAGTAAAAAACCGACGCAAAGCCGGTTTTTTTACATCCGGATTCGGACAATGCtgaa                     <  2:224629/249‑4 (MQ=255)
                                                                                                                                                                                                       tGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATg                 >  2:317613/1‑251 (MQ=255)
                                                                                                                                                                                                        gACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAata                   <  1:578856/248‑1 (MQ=255)
                                                                                                                                                                                                                      aaGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGa  <  1:3135/251‑1 (MQ=255)
                                                                                                                                                                                                                             tATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttt                                                 >  2:340732/1‑197 (MQ=255)
                                                                                                                                                                                                                             tATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttt                                                 <  1:340732/197‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
GAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCGCTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGA  >  NC_000913/3647725‑3648188

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: