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

GAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCG‑CTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGA  >  NC_000913/3647725‑3648188
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
aaaaGCGGATCAGTTCCAGCAGCGCAATTAACGCACCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATTCGATGGAGTGAGAATGCATCCG‑CTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCag                                                                                                                                                                                                                       <  2:571267/250‑1 (MQ=255)
      ggATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCCTGCGATGGAGTGAGAATGCATCCG‑CTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTa                                                                                                                                                                                                                  >  1:48503/1‑250 (MQ=255)
              accagcagctacaTTACCGCCCCTCGACCGATTATTTCTTTATCAAGCGCTCGGTCTCTCCGTGCGATGGAGTGAGAATGAGTCCG‑CTTACTCATCCACTGCCTGTCACGGCGCATTTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCAAGAGCATCAGTTTACGTCACCATCCGTACTTCTGACCAGGACATTGAAATCGTTTATGTTCATTGCCAGAGTGAAAATGCCTCAGAGCTAGCTGGATa                                                                                                                                                                                                          <  2:17888/238‑1 (MQ=255)
                     gcgcAATTAACGCCACTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCG‑CTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACa                                                                                                                                                                                                  <  2:240393/251‑1 (MQ=255)
                      cgcAATTATCGTCCCTAGATCGATGATTTCTTTATCAATCACCCGTTTGCTCCATGCGATGGAGTGCGAATGCGTCCG‑CTTACTCTTCCACTGCCTGACACTGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGTCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATTTGCATTGCTATAGTGAAAATGCCTTAGAGCTAACTGGATAATCATACAg                                                                                                                                                                                                 <  2:427990/251‑1 (MQ=255)
                              aCGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCG‑CTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCa                                                                                                                                                                                         <  1:479137/251‑1 (MQ=255)
                                gCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGCGGGAGAATGCATCCG‑CTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGCGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGGGAAATTGCCTCAGAGCTAACTGGATAATAATACAATACATGAAgg                                                                                                                                                                                       >  2:358311/1‑251 (MQ=255)
                                gCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCG‑CTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAgg                                                                                                                                                                                       >  1:120911/1‑251 (MQ=255)
                                  cccTAGAACGATGCTTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCG‑CTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGtt                                                                                                                                                                                     >  2:578856/1‑251 (MQ=255)
                                  cccTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGCGTGAGAATGCATCCG‑CTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATAATACAGTACATGCAGGtt                                                                                                                                                                                     >  2:405127/1‑251 (MQ=255)
                                   ccTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCG‑CTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTta                                                                                                                                                                                    <  1:102878/251‑1 (MQ=255)
                                           gatgatTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCG‑CTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGCGCTAACTGGATAATCATACAGTACATGCAGGTTATaaaccc                                                                                                                                                                             >  2:254096/1‑247 (MQ=255)
                                              gatTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCG‑CTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAg                                                                                                                                                                                                 >  1:200853/1‑227 (MQ=255)
                                              gatTGCTTTATCAATCACCCGTTTTCTCCAGGCGATGGAGTGAGAATGCATCCG‑CTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACGCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCAGTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAg                                                                                                                                                                                                 <  2:200853/227‑1 (MQ=255)
                                               atTGTTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGTGAATTCAACCG‑CTTACTCATCCACTGCCTTGCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCCTCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGTATAATCATACAGTACATGCTGGTTATCAACCCAGCAc                                                                                                                                                                        <  2:110674/251‑1 (MQ=255)
                                                  gCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCG‑CTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTc                                                                                                                                                                     >  1:527198/1‑251 (MQ=255)
                                                     ttATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCG‑CTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCt                                                                                                                                                                   <  1:557937/250‑1 (MQ=255)
                                                          aaccccctggtgtctcCATGCGATGGATGGAGAATGATTCCG‑CATACTCATCGAATGCCTGTCACGGGGCATGACTCAGTGGTAGATAAGAGTTCAATCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTTACAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCAGAGAGTACAGGCAGGTTATAAAATCAGCATGTCCTTGCAAt                                                                                                                                                             <  2:444857/239‑1 (MQ=255)
                                                             caCCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCG‑CTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGTATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGt                                                                                                                                                          <  1:297781/251‑1 (MQ=255)
                                                              ccccGGTTTCTCAATGCGAAGGAGTGAGAATGCATCCG‑CTTACTCATCCACTGCCGGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTCCTTCTGACCAGGACTTTGAAAGCGTTTATGCGCATTGCTAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGtt                                                                                                                                                         <  2:147521/250‑1 (MQ=255)
                                                                      tctcCATTCGAGGGAGTGAGAATGCATCCGACTT‑CTCATCCTCGTCCTGTAACGTCGCATGTCGAATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACAGGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTAt                                                                                                                                                  <  2:128802/250‑1 (MQ=255)
                                                                              cGATGGAGTGAGAATGCATCCG‑CTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGc                                                                                                                                         >  1:362834/1‑251 (MQ=255)
                                                                                       gagaATGCATCCG‑CTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCg                                                                                                                                <  1:625245/251‑1 (MQ=255)
                                                                                         gaATGCATCCG‑CTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGtt                                                                                                                              <  1:661315/251‑1 (MQ=255)
                                                                                               aTCCG‑CTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGTGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGAGGCGTTAgatg                                                                                                                         >  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)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
GAAAGCGGATCAGTTCCAGCAGCGCAATTAACGCCCCTAGAACGATGATTGCTTTATCAATCACCCGTTTTCTCCATGCGATGGAGTGAGAATGCATCCG‑CTTACTCATCCACTGCCTGTCACGGCGCATGTCTCATTGTTAGATAAGAACTCTCTCACTCCGGCCAGAGCATCAGTTAACGGCACCACCCGTACTTCTGACCAGGACTTTGAAAGCGTTTATGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGA  >  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: ‑