Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,657,599 T→A 16.6% intergenic (+32/‑767) 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,5990TA16.6% 377.3 / 60.9 156intergenic (+32/‑767)hdeD/gadEacid‑resistance membrane protein/gad regulon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (74/56);  new base A (13/13);  total (87/69)
Fisher's exact test for biased strand distribution p-value = 5.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.96e-01

Only 100 of 157 total aligned reads displayed.

GAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTC  >  NC_000913/3657500‑3657697
                                                                                                   |                                                                                                  
gAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAATCACt                                                                                                    >  2:292527/1‑100 (MQ=255)
  aTTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGAACATAGTCACTTg                                                                                                  >  1:1389648/1‑100 (MQ=255)
   ttGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTg                                                                                                  <  1:3818140/99‑1 (MQ=255)
     gAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtg                                                                                               >  2:3555813/1‑100 (MQ=255)
     gAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtg                                                                                               >  1:1656550/1‑100 (MQ=255)
     gAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtg                                                                                               >  2:1274791/1‑100 (MQ=255)
       aCTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtgtg                                                                                             <  1:40143/100‑1 (MQ=255)
         tGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGAc                                                                                           <  2:2038110/100‑1 (MQ=255)
         tGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGAACATAGTCACTTGATGTGAc                                                                                           >  1:3161758/1‑100 (MQ=255)
          natatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGACCATAGTCACTTGATGTGACt                                                                                          >  2:1226043/2‑100 (MQ=255)
          gATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACt                                                                                          >  1:1127957/1‑100 (MQ=255)
            tatTTAGCGCCGCCAGTCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTa                                                                                         <  1:1496745/99‑1 (MQ=255)
            tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTAt                                                                                        <  2:3888707/100‑1 (MQ=255)
            tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTAt                                                                                        <  1:3692018/100‑1 (MQ=255)
            tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTAt                                                                                        >  2:1762369/1‑100 (MQ=255)
            tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTAt                                                                                        <  1:16152/100‑1 (MQ=39)
            tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCAACCGGAACATAGTCCCTTGATGTGACTa                                                                                         >  2:1591985/1‑99 (MQ=255)
               ttAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGa                                                                                      >  1:3277441/1‑99 (MQ=255)
               ttAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGAt                                                                                     <  1:887022/100‑1 (MQ=255)
               ttAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGAt                                                                                     <  2:2236201/100‑1 (MQ=38)
                 aGCGCTGCCAGCCTGTTCAGCTTTGCCAGTTTGTTCGTTAAGTAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGATcc                                                                                   <  1:100964/100‑1 (MQ=37)
                 aGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATcc                                                                                   >  2:1307229/1‑100 (MQ=255)
                  gcgcCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGATCCg                                                                                  <  2:2897238/100‑1 (MQ=38)
                   cgccgcCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGa                                                                                 <  2:2383996/100‑1 (MQ=255)
                    gccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGAt                                                                                <  2:1242579/100‑1 (MQ=255)
                    gccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGATCCggg                                                                                <  2:3166011/100‑3 (MQ=38)
                     ccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGAtt                                                                               <  1:619952/100‑1 (MQ=255)
                     ccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGAtt                                                                               <  1:548667/100‑1 (MQ=255)
                     ccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGAtt                                                                               <  1:1274805/100‑1 (MQ=255)
                     ccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGATCCgggt                                                                               <  2:363008/100‑4 (MQ=37)
                      cgccAGCCTGTTTGGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGATCCggg                                                                                <  2:455746/98‑3 (MQ=37)
                      cgccAGCCTGTTCAGCTTTGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTa                                                                              <  1:1365900/100‑1 (MQ=255)
                      cgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAGGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTa                                                                              <  1:209915/100‑1 (MQ=255)
                      cgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGAtt                                                                               <  1:1685106/99‑1 (MQ=255)
                       gccAGCCTTTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGTATCATAGTCACTTGATGTGACTATGATCCGAt                                                                                <  2:3547806/97‑1 (MQ=255)
                       gccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTaa                                                                             >  2:3254429/1‑100 (MQ=255)
                       gccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTaa                                                                             >  2:1631128/1‑100 (MQ=255)
                       gccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGAACATAGTCACTTGATGTGACTATGATCCGATTaa                                                                             >  1:3447228/1‑100 (MQ=255)
                        ccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGATCCgggt                                                                               <  1:576516/97‑4 (MQ=37)
                           gCCTGTTCAGTTTCGTTAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATAct                                                                         <  1:1099071/100‑1 (MQ=255)
                            ccTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActc                                                                        <  2:2898675/100‑1 (MQ=255)
                                 tCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTcc                                                                   <  2:1790362/100‑1 (MQ=255)
                                  cAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCg                                                                  <  2:1202908/100‑1 (MQ=255)
                                    gCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCt                                                                >  2:3359139/1‑100 (MQ=255)
                                       tCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACg                                                             >  2:2034593/1‑100 (MQ=255)
                                       tCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATAATCTCTACGATACg                                                             >  1:2086061/1‑100 (MQ=255)
                                         gCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGAACATAGTCACTTAATGTGACTATGATCCGATTAATACTCTCTCCGCTACGc                                                            >  2:2959632/1‑99 (MQ=255)
                                             gTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAgtgt                                                       <  1:1197051/100‑1 (MQ=255)
                                             gTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAgtg                                                        <  1:530657/99‑1 (MQ=255)
                                                tGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGtgttgt                                                    >  2:2609386/1‑100 (MQ=255)
                                                tGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGtgttgt                                                    <  1:2719831/100‑1 (MQ=255)
                                                   tCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGa                                                 >  1:1053253/1‑100 (MQ=255)
                                                    cGTTAAGCAGCATTAATTACCCCGGTTGTCACACGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGAt                                                <  2:1163471/100‑1 (MQ=255)
                                                    cGTTAAGCAGCAAGAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGAt                                                >  1:2662068/1‑100 (MQ=255)
                                                     gTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATc                                               <  2:3588354/100‑1 (MQ=255)
                                                      ttAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCa                                              >  2:2615160/1‑100 (MQ=255)
                                                        aaGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAAt                                            >  2:1199937/1‑100 (MQ=255)
                                                         cgcagcaATAATTACCCCGGTTGTCACCCGGATCATATTCACTTGATGTGACTATAATTCGAATAATACTCTCTCAGCTACGCAGAGTTGTAGATCAAtt                                           >  1:618041/2‑100 (MQ=255)
                                                         agcagcAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAAtt                                           >  1:2798877/1‑100 (MQ=255)
                                                           cagcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTgc                                         >  2:3563147/1‑100 (MQ=255)
                                                            agcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTgcg                                        >  2:1614831/1‑100 (MQ=255)
                                                             gcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTgcgc                                       >  1:3348642/1‑100 (MQ=255)
                                                             gcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTgcgc                                       >  1:3474758/1‑100 (MQ=255)
                                                             gcaATAATTACCCCGGTTATAACCCGGAAAATAGTAACTTGATGTGACTATGATCAGAATAATACTCTCTCCGCTACGCAGTATTGTAGATTAAt                                            >  2:383738/1‑95 (MQ=255)
                                                               aataatTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCattttagtgc                                     >  2:3436310/1‑91 (MQ=255)
                                                               aataatTACCCCGGTTGTCACCCGAATCATAGTCACTTGATGTGACTATAATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTgcgc                                       <  2:2532367/98‑1 (MQ=255)
                                                                 taatTACCCCGGTTGTCACCCGGAACATAGTCACTTGATGTGACTATGATCCGATTAATAGTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTa                                   >  2:2715852/1‑100 (MQ=255)
                                                                  aatTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTAt                                  >  2:1496962/1‑100 (MQ=255)
                                                                  aatTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTAt                                  <  2:2845186/100‑1 (MQ=255)
                                                                  aatTACCCCGGTTGTCACCCAGATCATAGTCACATGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTAt                                  <  1:801747/100‑1 (MQ=255)
                                                                   atTACCCCGGTTGTCACTCGGATCATAGTCACTTGATGTGACTATGATCAGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTAt                                  <  1:1031838/99‑1 (MQ=255)
                                                                   atTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATc                                 >  1:3809918/1‑100 (MQ=255)
                                                                    ttACCCCGGTTGTCACTCGGATCATAGTCATTTGATGTGACTATGATCTGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCa                                <  2:2158743/100‑1 (MQ=255)
                                                                       ccccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTg                             >  1:2920650/1‑100 (MQ=255)
                                                                       ccccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTg                             >  2:3200373/1‑100 (MQ=255)
                                                                       ccccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTg                             <  1:352640/100‑1 (MQ=255)
                                                                       ccccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTg                             >  2:3054137/1‑100 (MQ=255)
                                                                         ccGGTTGGCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGaa                           >  1:2055625/1‑100 (MQ=255)
                                                                          cGGTTGTCTCCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGaaa                          <  1:3889920/100‑1 (MQ=255)
                                                                          cGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGaaa                          <  2:1424250/100‑1 (MQ=255)
                                                                          cGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATACGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGaa                           <  1:1753215/99‑1 (MQ=255)
                                                                                tCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTa                    >  2:1075606/1‑100 (MQ=255)
                                                                                tCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTa                    <  2:1918538/100‑1 (MQ=255)
                                                                                 taatcGGATCATAGTCACATCAAGTGACTATGATCCGGTTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAc                   >  1:3572796/5‑100 (MQ=25)
                                                                                 taatcGGATCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAc                   >  2:564778/5‑100 (MQ=37)
                                                                                   cccTGATCATAGTCACTTGTTGTGTTTATGTTTCGATTAATACTCTTTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAcc                  <  1:1294355/99‑1 (MQ=255)
                                                                                   cccGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCCACGCAGTATTGTAGATCAATAACGCACTATCATTGAAATAATTACCt                 >  1:153784/1‑100 (MQ=255)
                                                                                     cGGATCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGc               >  2:677885/1‑100 (MQ=38)
                                                                                      ggATCATAGTCACTTGATGTGACTATGATTCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCt              <  1:1049186/100‑1 (MQ=255)
                                                                                      ggATCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCt              >  1:1786420/1‑100 (MQ=38)
                                                                                      ggATCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCt              >  1:994866/1‑100 (MQ=38)
                                                                                       gATCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTa             >  1:2825885/1‑100 (MQ=38)
                                                                                          cATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTg          >  1:2397039/1‑100 (MQ=39)
                                                                                          cATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGATACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTg          >  2:2961422/1‑100 (MQ=38)
                                                                                           aTAGTCACATGAAGTGACTATAATCCGATTAATACACTCTCCGCTACGCAGTGTTGTATCTCAAATACGCACTAACATTGAAATAAATACCTTCTaataa         >  2:3282162/1‑96 (MQ=255)
                                                                                            tAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTg          <  1:219319/98‑1 (MQ=255)
                                                                                            tAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGat        <  1:3422341/100‑1 (MQ=255)
                                                                                            tAGTCACTTGATGTGACTATAATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGat        >  2:637454/1‑100 (MQ=255)
                                                                                              gTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGatta      >  2:3760994/1‑100 (MQ=255)
                                                                                                  cTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTc  <  2:2441683/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
GAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTC  >  NC_000913/3657500‑3657697

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