Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,657,603 T→A 15.2% intergenic (+36/‑763) hdeD → / → gadE acid‑resistance membrane protein/gad regulon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,657,6030TA15.2% 483.1 / 61.0 170intergenic (+36/‑763)hdeD/gadEacid‑resistance membrane protein/gad regulon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (77/67);  new base A (11/15);  total (88/82)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 182 total aligned reads displayed.

GAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCT  >  NC_000913/3657505‑3657702
                                                                                                  |                                                                                                   
gAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtg                                                                                                    <  2:2012404/100‑1 (MQ=255)
gAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtg                                                                                                    <  1:3025454/100‑1 (MQ=255)
   cTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGa                                                                                                 >  1:2405141/1‑100 (MQ=255)
   cTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGa                                                                                                 >  2:1287653/1‑100 (MQ=255)
    tGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGAc                                                                                                <  2:3612310/100‑1 (MQ=255)
       tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTAt                                                                                             <  1:1132295/100‑1 (MQ=255)
       tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTAt                                                                                             >  2:295319/1‑100 (MQ=255)
       tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTAt                                                                                             <  1:561520/100‑1 (MQ=255)
       tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTAt                                                                                             <  1:3703276/100‑1 (MQ=39)
       tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTAt                                                                                             <  2:3051207/100‑1 (MQ=39)
        atTTATCGCTGCCAGCCTGTTCAACTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACGCCGGTTGTTACCCGGATCATAGTCACATCAAGTGACTAt                                                                                             <  1:2423608/99‑1 (MQ=37)
        atTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATg                                                                                            >  1:657037/1‑100 (MQ=255)
         tttAGCGCCGCCAGCCTGTTCAGCTTCGCAAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATTATCGTAACTTGATGTGACTATGa                                                                                           >  2:2587557/1‑100 (MQ=255)
            aGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATTATAGTCACTTGATGTGACTATGATcc                                                                                        <  2:2249335/100‑1 (MQ=255)
            aGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATcc                                                                                        >  2:2803923/1‑100 (MQ=255)
            aGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATcc                                                                                        >  2:2803840/1‑100 (MQ=255)
            aGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATcc                                                                                        >  2:1475021/1‑100 (MQ=255)
            aGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGATcc                                                                                        <  2:1514867/100‑1 (MQ=38)
             gcgcCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGATCCg                                                                                       <  1:3588592/100‑1 (MQ=255)
              cgccgcCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGa                                                                                      >  2:3260924/1‑100 (MQ=255)
               gccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGAt                                                                                     <  2:1737060/100‑1 (MQ=255)
               gccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGATCCggg                                                                                     <  1:1568851/100‑3 (MQ=38)
               gccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGATCCggg                                                                                     <  2:2334421/100‑3 (MQ=37)
                ccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGAtt                                                                                    >  1:1785300/1‑100 (MQ=255)
                  gccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAGGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTAAGATCCGATTaa                                                                                  <  1:1940472/100‑1 (MQ=255)
                       ccTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActc                                                                             <  2:3386642/100‑1 (MQ=255)
                        cTGTTCAGCTTCGCCAGTTTGTTCGTTACGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGCCTATGATCCGATTAATActc                                                                             <  1:1316044/99‑1 (MQ=255)
                        cTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActct                                                                            <  1:3273719/100‑1 (MQ=255)
                         tGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActctc                                                                           >  2:503103/1‑100 (MQ=255)
                            tCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTcc                                                                        >  2:984510/1‑100 (MQ=255)
                            tCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTcc                                                                        <  2:3277691/100‑1 (MQ=255)
                            tCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTcc                                                                        >  2:308406/1‑100 (MQ=255)
                            tCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTcc                                                                        >  1:3465340/1‑100 (MQ=255)
                               gCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCt                                                                     >  2:263642/1‑100 (MQ=255)
                                 ttCGCCAGTTTGTTCTTTAAGCAGCAATAATTACCCCGGTTGTCACCCTTATCATAGTCACTTCATGTGACTATGTTCCGATTAATACTCTCTCCGCTAc                                                                   <  1:497520/100‑1 (MQ=255)
                                 ttCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCTGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTAc                                                                   <  1:1641205/100‑1 (MQ=255)
                                 ttCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTAc                                                                   <  2:3536889/100‑1 (MQ=255)
                                 ttCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTAc                                                                   <  2:1930225/100‑1 (MQ=255)
                                  tCGCGAGTTTGTTCGTTAAGCAGCAAGAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTAc                                                                   <  1:3505985/99‑1 (MQ=255)
                                    gCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCa                                                                >  1:2789670/1‑100 (MQ=255)
                                     ccAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAg                                                               >  2:2047269/1‑100 (MQ=255)
                                       aGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAgtg                                                             >  2:2464743/1‑100 (MQ=255)
                                        gTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAgtgt                                                            <  2:456889/100‑1 (MQ=255)
                                        gTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAgtgt                                                            <  1:726328/100‑1 (MQ=255)
                                        gTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAgtgt                                                            >  2:381138/1‑100 (MQ=255)
                                          ttGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGtgttg                                                          <  1:1401455/100‑1 (MQ=255)
                                            gTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTa                                                        <  1:1899094/100‑1 (MQ=255)
                                             ttCGTTAAGCAGCAATAATTACTCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACGCTCTCCGCTACGCAGtgttg                                                          <  1:2251667/97‑1 (MQ=255)
                                             ttCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATTATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAg                                                       <  1:1985745/100‑1 (MQ=255)
                                              tCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGa                                                      <  1:897541/100‑1 (MQ=255)
                                                 ttAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCa                                                   <  1:137409/100‑1 (MQ=255)
                                                   aaGCAGCAATAATTACCCCGGTTGTCACCCGGAACATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAAt                                                 >  1:626986/1‑100 (MQ=255)
                                                    agcagcAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAAtt                                                >  1:2919975/1‑100 (MQ=255)
                                                    agcagcAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAAtt                                                <  2:2740765/100‑1 (MQ=255)
                                                    agcagcAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAAtt                                                >  2:3472466/1‑100 (MQ=255)
                                                      cagcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGGCTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTgc                                              >  2:111365/1‑100 (MQ=255)
                                                       agcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTgcg                                             <  1:252758/100‑1 (MQ=255)
                                                        gcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTgcgc                                            >  1:2814751/1‑100 (MQ=255)
                                                        gcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTgcgc                                            <  1:2891172/100‑1 (MQ=255)
                                                        gcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTgcgc                                            >  1:950651/1‑100 (MQ=255)
                                                        gcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTgcgc                                            <  1:3015806/100‑1 (MQ=255)
                                                         caATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCa                                           >  1:1990429/1‑100 (MQ=255)
                                                         caATAATTACCCCGGTTGTCACCCGGAGCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTGCGCTACGCAGTGTTGTGGATCAATTGCGCa                                           >  2:2523571/1‑100 (MQ=255)
                                                           ataatTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACt                                         >  1:929220/1‑100 (MQ=255)
                                                             aatTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTAt                                       <  2:1962711/100‑1 (MQ=255)
                                                              atTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATc                                      <  2:3516250/100‑1 (MQ=255)
                                                               ttACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCa                                     <  1:1475015/100‑1 (MQ=255)
                                                                  ccccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCAtt                                   >  2:1987538/1‑99 (MQ=255)
                                                                  ccccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTg                                  <  1:2401070/100‑1 (MQ=255)
                                                                  ccccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTg                                  <  2:1175671/100‑1 (MQ=255)
                                                                  ccccGGTTGTCACCCGGATAATAGTCAGTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTg                                  <  2:2062175/100‑1 (MQ=255)
                                                                   cccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGa                                 >  1:1752693/1‑100 (MQ=255)
                                                                    ccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTg                                  <  1:3157256/98‑1 (MQ=255)
                                                                    ccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGaa                                >  2:2074/1‑100 (MQ=255)
                                                                    ccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGaa                                >  1:2612457/1‑100 (MQ=255)
                                                                    ccGGTTGTCACCCGGAACATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCCCCGCTACGCAGGGTTGTAGATCAATTGCGCACTATCATGGaa                                >  1:218949/1‑100 (MQ=255)
                                                                     cGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGaaa                               <  2:2798784/100‑1 (MQ=255)
                                                                     cGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGaaa                               <  2:1835271/100‑1 (MQ=255)
                                                                           tCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTa                         >  2:1863201/1‑100 (MQ=255)
                                                                            cACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAc                        >  2:3657734/1‑100 (MQ=255)
                                                                            cACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAc                        >  2:872493/1‑100 (MQ=255)
                                                                             aCCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAcc                       <  2:657042/100‑1 (MQ=255)
                                                                              cccGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCt                      >  1:2388922/1‑100 (MQ=255)
                                                                              cccGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCt                      >  2:3285363/1‑100 (MQ=255)
                                                                              atcGGAACATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCACTACGAAGTGTTGTAGATCAATTGCGCACTATCATTGAaata                             >  1:2907132/3‑93 (MQ=25)
                                                                               ccGGATCATAGTCACTTGATGTGACTATGATCCGAGTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTg                     >  1:3052642/1‑100 (MQ=255)
                                                                                 ggATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCt                   <  2:1461721/100‑1 (MQ=255)
                                                                                    tCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGt                >  1:1218800/1‑100 (MQ=39)
                                                                                     cATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTg               >  1:1073343/1‑100 (MQ=39)
                                                                                     cATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTg               >  2:793396/1‑100 (MQ=255)
                                                                                       tAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGat             <  1:1049821/100‑1 (MQ=255)
                                                                                       tAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGat             >  2:1359327/1‑100 (MQ=255)
                                                                                             cTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTc       >  2:379074/1‑100 (MQ=255)
                                                                                              ttGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCa      <  2:3365661/100‑1 (MQ=255)
                                                                                               tGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCaa     >  2:2666474/1‑100 (MQ=255)
                                                                                               tGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCa      <  1:2952057/99‑1 (MQ=255)
                                                                                                gATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAAc    >  2:1260514/1‑100 (MQ=255)
                                                                                                 aTGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAAcc   <  1:1718466/100‑1 (MQ=255)
                                                                                                  tgtgACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCt  <  1:1870996/100‑1 (MQ=255)
                                                                                                  tgtgACTATGATCCGATTAATACTCTCTCCGCCACGCTGTGTTGTAGAGCAATTGCGGACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAAcc   <  1:2011536/99‑1 (MQ=255)
                                                                                                  |                                                                                                   
GAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCT  >  NC_000913/3657505‑3657702

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