Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,926,894 C→T G37E (GGG→GAG)  asnC ← DNA‑binding transcriptional dual regulator AsnC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,926,8940CT97.6% 125.9 / ‑3.9 41G37E (GGG→GAG) asnCDNA‑binding transcriptional dual regulator AsnC
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (14/26);  total (14/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.54e-01

AGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGT‑AATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGAC  >  NC_000913/3926648‑3927139
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
agagCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGT‑AATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCgg                                                                                                                                                                                                                                                   >  1:123745/1‑251 (MQ=255)
  ggCGTCGAACGACCGGCACATCCTTTTTATAAAGATGCTGTAGTGGCCGGTTGTGT‑AATAGGCTTCAGTGATTTCAGCAAGTCTTTCCAGCTTTGCCAGCCCGGAAGGGTAGTCTTTGGCGCTCTTTAATAAAATGCCGATAAAGCCGCCTACGTCATAACCGGGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGAc                                                                                                                                                                                                                                                 <  2:121196/250‑1 (MQ=255)
   gCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGCAGTGGCCGGTTGTGT‑AATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGACGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGCGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGACTCGTCTCCGGACt                                                                                                                                                                                                                                                >  2:490787/1‑251 (MQ=255)
   gCGTCGATAGAACGGCAAATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGT‑AATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTGTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGACt                                                                                                                                                                                                                                                <  1:321517/251‑1 (MQ=255)
         atcgaACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGT‑AATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGACTGACGcc                                                                                                                                                                                                                                          >  1:537407/1‑251 (MQ=255)
                             tataAAGATGCTGTAGTGGCCGGTTGTGT‑AATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCg                                                                                                                                                                                                                       <  2:144241/250‑1 (MQ=255)
                                  aGATGCTGTAGGGGCCG‑TTTTGTGAATAGGCTTCAGTTACTTCATCAAGGTTTTCCAGCTTGGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTAAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAAGACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTaa                                                                                                                                                                                                                 <  2:336061/251‑1 (MQ=255)
                                    aTGCTGTAGTGGCCGGTTGGGT‑AATAGGCTTCAGTGACTTCATCAAGGTTTTCCAGCTTTGCTAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCAATAAAGCAGCTTACGCCATATCCGGGCTGCGCCGGGTTGACATCAATCCGCGCCCCGGGAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAGCGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGc                                                                                                                                                                                                               <  1:347220/251‑1 (MQ=255)
                                       cTGTAGTGGCCGGTTGTGT‑AATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGt                                                                                                                                                                                                            >  1:165305/1‑251 (MQ=255)
                                                      gtgt‑AATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATCCGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTTAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCa                                                                                                                                                                                              <  1:560880/250‑1 (MQ=255)
                                                              aGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGc                                                                                                                                                                                       <  1:78977/250‑1 (MQ=255)
                                                                   tCAGTGTCTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTCCTTCGGGATGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCa                                                                                                                                                                                  <  1:321148/250‑1 (MQ=255)
                                                                        gACTTCATCAAGGCTTCCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCAATAAAGCACCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGc                                                                                                                                                                            <  1:180702/251‑1 (MQ=255)
                                                                                  aGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTcc                                                                                                                                                                    >  1:448748/1‑249 (MQ=255)
                                                                                   ggCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTcc                                                                                                                                                                    >  1:187404/1‑248 (MQ=255)
                                                                                                 gCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGAt                                                                                                                                                    <  1:435371/250‑1 (MQ=255)
                                                                                                     gcgcGGAAGGGTAGTCTTGGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGTTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTGCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAAttt                                                                                                                                               <  2:238139/251‑1 (MQ=255)
                                                                                                                 agtgtgtggggcTCTTTAATATAATGCTGTAAAAGCAGCTTACGTCATACCCGAGCGGCTTCGGGTTGACAGCAATACGCGCCCCGTTAAGGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGAtt                                                                                                                                    <  2:471760/240‑1 (MQ=255)
                                                                                                                      ttGGCGCTCTTTAATATAATGCCGATAAAGCACCCTACGTAATATCCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCt                                                                                                                                <  2:165305/249‑1 (MQ=255)
                                                                                                                                ttAATATAAGGCCGATAAAGCAGCCTACTGCATAACCGGGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGAGTGTCGATCAGATAATTTTCCATAGGATGTTCTTATGCGGATtg                                                                                                                     <  1:490787/250‑1 (MQ=255)
                                                                                                                                        aTGCCGCTAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCg                                                                                                                                                                                                                               >  2:276041/1‑136 (MQ=255)
                                                                                                                                          gCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGattcattc                                                                                                          <  1:464049/251‑1 (MQ=255)
                                                                                                                                           ccGATAAAGCAGCCTACGTCATAACCTGGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCt                                                                                                         <  2:537407/251‑1 (MQ=255)
                                                                                                                                                aaaGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTAtttt                                                                                                    <  1:551334/251‑1 (MQ=255)
                                                                                                                                                   gcagcCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGc                                                                                                 >  1:57000/1‑251 (MQ=255)
                                                                                                                                                         tACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCcttctt                                                                                           >  1:16879/1‑251 (MQ=255)
                                                                                                                                                          ccGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCttt                                                                                          <  1:186598/250‑1 (MQ=255)
                                                                                                                                                                 aaCCGAGCTGCTTCGGGCTGACATCACTACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGCCTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAAtt                                                                                   >  2:386221/1‑250 (MQ=255)
                                                                                                                                                                  aCCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGa                                                                                  <  1:459918/251‑1 (MQ=255)
                                                                                                                                                                   ccGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGaa                                                                                 <  1:22396/251‑1 (MQ=255)
                                                                                                                                                                                                   accggcacggATCACCGCCTGCTTCCTTTTCTCTACTCGAACGTGACTCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTCAGCGGTGCGGGCATTGCCCATTACTGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGAGAGTTAGGCTTTTTATTGAATGat                                                 >  2:86342/10‑251 (MQ=255)
                                                                                                                                                                                                    ccGGTAGTGATCCCCGCCTGCTTCAGTTTCTCTACCCGAACGTGAATCGTCTGCGGACTGACGCCAAATTGTTTCGCCAGTTAGGCGTAAGCGGTGCGCTCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGTCTTCTTTTTTAATGAATCAAAAGTGAGGTAGGCGTTTTTTTGAATGatt                                                <  2:573548/251‑1 (MQ=255)
                                                                                                                                                                                                      ggTAATGATCCCCGCCTGCTTCAGTTTCTCTACTCGAACGTGAATCGTCTCCGGACTGACGCCAAATTGTTGCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAAGTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCGTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGattat                                              <  1:228652/251‑1 (MQ=255)
                                                                                                                                                                                                            gATCTCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTACAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCAt                                         <  1:459630/250‑1 (MQ=255)
                                                                                                                                                                                                            gATCCCCGCCTGCTTCATTTTCTCTACTCGCACGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCAtg                                        >  2:90451/1‑251 (MQ=255)
                                                                                                                                                                                                               ccccGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCAtgtgt                                     >  2:445081/1‑251 (MQ=255)
                                                                                                                                                                                                                              tttCTCTACTCGAACGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTttgttg                       <  1:93617/250‑1 (MQ=255)
                                                                                                                                                                                                                                   ctACTCGAACGGGAATCGTCTCCGGACTGACGCCAAATTGTTGCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTACAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCGTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCtt                    <  1:88767/248‑1 (MQ=255)
                                                                                                                                                                                                                                        cGAACGTGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATaa            <  2:187404/251‑1 (MQ=255)
                                                                                                                                                                                                                                              tGAATCGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAAc       >  1:372858/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   cGTCTCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGAc  <  2:448748/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
AGAGCGTCGATCGAACGGCACATCACTTTTATAAAGATGCTGTAGTGGCCGGTTGTGT‑AATAGGCTTCAGTGACTTCATCAAGGCTTTCCAGCTTTGCCAGCGCGGAAGGGTAGTCTTTGGCGCTCTTTAATATAATGCCGATAAAGCAGCCTACGTCATAACCGAGCTGCTTCGGGCTGACATCAATACGCGCCCCGGTAATGATCCCCGCCTGCTTCATTTTCTCTACTCGAACGTGAATCGTCCCCGGACTGACGCCAAATTGTTTCGCCAGTTCGGCGTAAGCGGTGCGCGCATTGCCCATTAATGCTTCCAGGATGCCACGGTCCAGATTGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGAC  >  NC_000913/3926648‑3927139

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