Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,658,780 A→G 80.0% T139A (ACG→GCG)  gadE → gad regulon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,658,7800AG80.0% 502.2 / 105.9 240T139A (ACG→GCG) gadEgad regulon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (24/24);  new base G (85/106);  total (109/131)
Fisher's exact test for biased strand distribution p-value = 5.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

Only 100 of 244 total aligned reads displayed.

CACTTGCCCCATAAGAATTCACAATTATGTTTTTCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGACGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCA  >  NC_000913/3658681‑3658879
                                                                                                   |                                                                                                   
     gCCCCATAAGAATTCACAATTATGTTTTTCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTc                                                                                                >  2:1304532/1‑100 (MQ=255)
     gCCCCATAAGAATTCACAATTATGTTTTTCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTc                                                                                                <  1:2589071/100‑1 (MQ=255)
      ctccATAAGAATTCACAATTATCTTTTTCACATAATTAGTTAAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTc                                                                                                <  2:3535509/97‑1 (MQ=255)
       cccATAATAATTCACAATTATGTTTTTCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGACGCTCaa                                                                                              >  1:2175835/1‑100 (MQ=255)
       cccATAAGAATTCACAATTATGTTTTTCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCaa                                                                                              >  2:250341/1‑100 (MQ=255)
         cATAAGAATTCACAATTATGTTTTTCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACATCGGCGCTCAat                                                                                             <  2:26483/99‑1 (MQ=255)
          aTAAGAATTCACAATTATGTTTGTCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAata                                                                                            >  2:2416751/1‑99 (MQ=255)
            aaGAATTCACAATTATGTTTTTCACATAATCAGTTCAAAATTATGCAACTTATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATAttt                                                                                         <  1:3607555/100‑1 (MQ=255)
             aGAATTCACAATTATGTTTTTCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTc                                                                                        >  2:3655128/1‑100 (MQ=255)
             aGAATTCACAATTATGTTTTACACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGACGCTCAATATTTc                                                                                        >  2:2553709/1‑100 (MQ=255)
              gAATTCACAATTATGTTTTTCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCg                                                                                       <  1:3637617/100‑1 (MQ=255)
               aaTTCACAATTATGTTTTTCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTctc                                                                                      <  1:2805948/100‑3 (MQ=255)
               aaTTCACAATTATGTTTTTCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGc                                                                                      <  1:102593/100‑1 (MQ=255)
               aaTTCACAATTATGTTTTTCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGACGCTCAATATTTCGc                                                                                      <  2:1650173/100‑1 (MQ=255)
                aTTCACAATTATGTTTTTCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGca                                                                                     >  2:3231048/1‑100 (MQ=255)
                aTTCACAATTATGTTTTTCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGca                                                                                     >  2:2474536/1‑100 (MQ=255)
                aTTCACAATTATGTTTTTCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGca                                                                                     <  1:945880/100‑1 (MQ=255)
                 ttCACAATTATGTTTTTCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTGGcaa                                                                                    >  2:757362/1‑100 (MQ=255)
                     caATTATGTTTTTCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGACGCTCAATATTTCGCaacaaa                                                                                <  1:478224/100‑1 (MQ=255)
                      aaTTATGTTTTTCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCaacaaa                                                                                <  2:41783/99‑1 (MQ=255)
                       aTTATGTTTTTCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCaacaaaca                                                                              <  1:3857710/100‑1 (MQ=255)
                       aTTATGTTTTTCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCaacaaaca                                                                              >  1:3172087/1‑100 (MQ=255)
                        ttATGTTNTTCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAAAAAACAt                                                                             >  2:183291/1‑100 (MQ=255)
                         tATGTTTTTCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACAtt                                                                            >  1:975375/1‑100 (MQ=255)
                           tGTTTTTCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTaa                                                                          >  1:2523213/1‑100 (MQ=255)
                            gTTTTTCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTaaa                                                                         >  1:374424/1‑100 (MQ=255)
                            gTTTTTCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGAAATATCACGTTGGTGATCAATATTTCGCAATAAACATTaa                                                                          >  2:3174107/1‑99 (MQ=255)
                             tttttCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGACGCTCAATATTTCGCAACAAACATTAAAg                                                                        <  2:1520315/100‑1 (MQ=255)
                               tttCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGACGCTCAATATTTCGCAACAAACATTAAAGAt                                                                      <  1:1125352/100‑1 (MQ=255)
                               tgtCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACACACACTaaa                                                                         <  2:2590663/95‑1 (MQ=255)
                                ttCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGAtt                                                                     >  1:2395173/1‑100 (MQ=255)
                                 tCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTc                                                                    <  1:2880147/100‑1 (MQ=255)
                                  cacaTAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCa                                                                   >  1:636422/1‑100 (MQ=255)
                                  cacaTAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCa                                                                   <  1:3947981/100‑1 (MQ=255)
                                   acaTAATCAGTTCAAAATTATGCAAATGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCa                                                                   <  2:1162585/99‑1 (MQ=255)
                                    caTAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGa                                                                 >  2:2163775/1‑100 (MQ=255)
                                      tAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGACGCTCAATATTTCGCAACAAACATTAAAGATTCAGaaa                                                               >  2:627378/1‑100 (MQ=255)
                                       aaTCAGTTAAAAATGATGAAACGGATTCTGAAAAATAAAAATGAAAGCAATATAACGTCGATGCGCAATATTTCGCAACAAACATTAAAGATTCAGaa                                                                <  2:3973244/98‑1 (MQ=255)
                                          cAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCa                                                           <  1:1957891/100‑1 (MQ=255)
                                           aGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCaa                                                          >  1:1668661/1‑100 (MQ=255)
                                            gTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAAc                                                         <  1:4119339/100‑1 (MQ=255)
                                            gTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAAc                                                         >  1:1877238/1‑100 (MQ=255)
                                            gTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAAc                                                         <  1:1002371/100‑1 (MQ=255)
                                            gTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAAc                                                         <  2:945402/100‑1 (MQ=255)
                                            gTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGACGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAAc                                                         <  1:2447460/100‑1 (MQ=255)
                                             ttCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACa                                                        <  1:3476453/100‑1 (MQ=255)
                                               cAAAATTATGCAACTGATTCTGAAAAATAAAAATTAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACAt                                                       <  2:350750/99‑1 (MQ=255)
                                               cAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACAtt                                                      >  1:1224014/1‑100 (MQ=255)
                                                aaaaTTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAAcn                                                        <  2:1582925/97‑2 (MQ=255)
                                                aaaaTTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTAGCAACAAACATTAAAGATTCAGAAATTCAACATTa                                                     >  2:88542/1‑100 (MQ=255)
                                                    ctATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTa                                                 <  1:4082991/99‑1 (MQ=255)
                                                      aTGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTaag                                               >  2:4082959/1‑98 (MQ=255)
                                                      aTGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAAAATTTCGCAACAAACATTACAGATTCAGACATTCAACATTATGTACa                                               >  2:2862866/1‑100 (MQ=255)
                                                      aTGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGACGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACa                                               >  2:1404530/1‑100 (MQ=255)
                                                        gCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAg                                             <  2:3645434/100‑1 (MQ=255)
                                                         cAACTGATTCTGAAAAATAAAAATGAAAGCAATATCAAGTCGACGCTCAATATTTAGCAAAAAACATTAAAGATTCAGAAATTCAACATTATGTACaa                                              >  2:2232801/1‑98 (MQ=255)
                                                         cAACTGATACTGAAAAATAAAAATGAAAGCAATATCACGCAGACGCTCAATATTTAGCAACAAACATTAAAGATTCGGAAGTTGAACTTTATGTACAAGc                                            >  2:196402/1‑100 (MQ=255)
                                                             tGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGaa                                        <  1:229858/100‑1 (MQ=255)
                                                             tGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGaa                                        >  1:1539431/1‑100 (MQ=255)
                                                             tGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGaa                                        <  1:4000924/100‑1 (MQ=255)
                                                               aTTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAAc                                      <  2:133852/100‑1 (MQ=255)
                                                               aTTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAAc                                      <  2:2252471/100‑1 (MQ=255)
                                                               aTTCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAAc                                      <  1:2766505/100‑1 (MQ=255)
                                                               aTTCTGAAAAATAAAAATGAAAGCAATATCACGTCGACGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAAc                                      >  1:1734836/1‑100 (MQ=255)
                                                                ttCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACt                                     <  1:4134551/100‑1 (MQ=255)
                                                                 tCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTa                                    >  1:3692675/1‑100 (MQ=255)
                                                                 tCTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTa                                    >  2:1615199/1‑100 (MQ=255)
                                                                  cTGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTaa                                   >  2:3025061/1‑100 (MQ=255)
                                                                   tGAAAAATAAAAATGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAg                                  >  1:2628993/1‑100 (MQ=255)
                                                                   tGAAAAATAAAAATGAAAGCAATATCACGTCGACGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAg                                  >  1:3491819/1‑100 (MQ=255)
                                                                    gAAAAATAAAAATGAAAGCAATATCACGTCGTCGCTCAATATTTCGCAACAAACATTAAAGATTCATAAATTCAACATTATGTACAAGCTGAAACTAAGa                                 <  2:1710968/100‑1 (MQ=255)
                                                                      aaaataaaaatGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATGATGTACAAGCTGAAACTAAGACg                               <  1:3532514/100‑1 (MQ=255)
                                                                        aataaaaatGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTa                             >  1:2758664/1‑100 (MQ=255)
                                                                         ataaaaatGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTAt                            >  1:2247307/1‑100 (MQ=255)
                                                                             aaatGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGc                        >  2:2029737/1‑100 (MQ=255)
                                                                             aaatGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGc                        <  1:1180844/100‑1 (MQ=255)
                                                                             aaatGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGc                        >  1:3888165/1‑100 (MQ=255)
                                                                              aatGAAAGCAATATCACGTCTACGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCg                       <  2:3799641/100‑1 (MQ=255)
                                                                              aatGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCg                       <  2:143197/100‑1 (MQ=255)
                                                                              aatGAAAGCAATATCACGTCGGCGCTCAATATTTCGAAACAAACATTAAAGATTCAGAAATTTAACATCATGTACAAGCTGAAACTAAGACGTATGAGCg                       >  2:2805234/1‑100 (MQ=255)
                                                                              aatGAAAGCAATATCACGTCGACGCTCAATATTTCGCAACAAACATTAAAGATTCAGACATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCg                       <  2:3373438/100‑1 (MQ=255)
                                                                               atGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGa                      <  2:1536067/100‑1 (MQ=255)
                                                                               atGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGa                      <  1:3319361/100‑1 (MQ=255)
                                                                               atGAAAGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGACTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGa                      <  2:3219289/100‑1 (MQ=255)
                                                                                   aaGCAATATCACGTCGACGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATc                  >  2:956163/1‑100 (MQ=255)
                                                                                    aGCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCg                 <  1:1359329/100‑1 (MQ=255)
                                                                                     gCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGt                >  1:3092539/1‑100 (MQ=255)
                                                                                     gCAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGt                >  1:3823880/1‑100 (MQ=255)
                                                                                      cAATATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTc               >  1:1294949/1‑100 (MQ=255)
                                                                                       aaTATCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCa              >  2:1776977/1‑100 (MQ=255)
                                                                                          atCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCAccc           <  2:3888136/100‑1 (MQ=255)
                                                                                          atCACGTCGACGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGCCATCGTCAccc           >  1:3935541/1‑100 (MQ=255)
                                                                                           tCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCt          >  2:3792288/1‑100 (MQ=255)
                                                                                           tCACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCt          >  1:4070856/1‑100 (MQ=255)
                                                                                            cACGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTg         <  2:3318391/100‑1 (MQ=255)
                                                                                              cGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTGACCCTggg       <  2:3288403/100‑1 (MQ=255)
                                                                                              cGTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTggg       <  1:3552824/100‑1 (MQ=255)
                                                                                              cGTCGACGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTggg       <  2:205776/100‑1 (MQ=255)
                                                                                               gTCGGCGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGt      <  1:2095398/100‑1 (MQ=255)
                                                                                               gTCGGCGATCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTgg        <  2:2772060/98‑1 (MQ=255)
                                                                                                   |                                                                                                   
CACTTGCCCCATAAGAATTCACAATTATGTTTTTCACATAATCAGTTCAAAATTATGCAACTGATTCTGAAAAATAAAAATGAAAGCAATATCACGTCGACGCTCAATATTTCGCAACAAACATTAAAGATTCAGAAATTCAACATTATGTACAAGCTGAAACTAAGACGTATGAGCGACATCGTCACCCTGGGTATCA  >  NC_000913/3658681‑3658879

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: