Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,657,601 G→C 16.7% intergenic (+34/‑765) 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,6010GC16.7% 457.1 / 82.8 175intergenic (+34/‑765)hdeD/gadEacid‑resistance membrane protein/gad regulon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (81/64);  new base C (14/15);  total (96/79)
Fisher's exact test for biased strand distribution p-value = 5.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.11e-01

Only 100 of 180 total aligned reads displayed.

ATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAAC  >  NC_000913/3657502‑3657700
                                                                                                   |                                                                                                   
 ttGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGa                                                                                                    >  1:152104/1‑100 (MQ=255)
  tGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTCCCCCGGTTGTCACCAGGATCATAGTCACTTGa                                                                                                    <  2:453387/99‑1 (MQ=255)
   gAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtg                                                                                                  <  2:2275073/100‑1 (MQ=255)
   gAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtg                                                                                                  >  2:3207463/1‑100 (MQ=255)
    aaCTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtgt                                                                                                 >  2:3356901/1‑100 (MQ=255)
    aaCTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGAGCATAGTCACTTGAtgt                                                                                                 >  1:3095215/1‑100 (MQ=255)
     aCTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtgtg                                                                                                <  1:2310399/100‑1 (MQ=255)
     aCTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtgtg                                                                                                <  1:1386052/100‑1 (MQ=255)
        gATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACt                                                                                             >  2:3227382/1‑100 (MQ=255)
         atatTTAGCGTCGTCAGCCTGTTCACCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTa                                                                                            <  2:1387476/100‑1 (MQ=38)
         atatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTa                                                                                            <  2:3711914/100‑1 (MQ=255)
          tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTAt                                                                                           <  2:3108958/100‑1 (MQ=255)
          tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTAt                                                                                           >  2:3015456/1‑100 (MQ=255)
           atTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATg                                                                                          <  1:2003172/100‑1 (MQ=39)
            tttAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGa                                                                                         <  1:768875/100‑1 (MQ=38)
            tttAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGa                                                                                         <  1:3792101/100‑1 (MQ=38)
            tttAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGa                                                                                         <  1:2363802/100‑1 (MQ=39)
              tAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATc                                                                                       >  1:471211/1‑100 (MQ=255)
              tAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATc                                                                                       <  2:2027080/100‑1 (MQ=255)
              tAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGATc                                                                                       <  1:176278/100‑1 (MQ=255)
                gcgcCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCg                                                                                     <  1:3327585/100‑1 (MQ=255)
                gcgcCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCg                                                                                     <  2:3873770/100‑1 (MQ=255)
                gcgcCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGATCCg                                                                                     <  1:3116063/100‑1 (MQ=38)
                 cgccgcCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGa                                                                                    >  2:2180831/1‑100 (MQ=255)
                 cgccgcCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGa                                                                                    >  1:3163157/1‑100 (MQ=255)
                   ccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGAtt                                                                                  >  1:3657434/1‑100 (MQ=255)
                   ccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGAtt                                                                                  <  2:1114720/100‑1 (MQ=255)
                    cgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGATCCgggtg                                                                                 <  1:1302863/100‑5 (MQ=37)
                     gccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTaa                                                                                >  1:1343663/1‑100 (MQ=255)
                     gccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTaa                                                                                >  2:2907229/1‑100 (MQ=255)
                      ccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAAt                                                                               <  2:3593855/100‑1 (MQ=255)
                      ccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAAt                                                                               <  2:1179745/100‑1 (MQ=255)
                         gCCTTTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGTAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATAct                                                                            <  2:62126/100‑1 (MQ=255)
                           cTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActct                                                                          >  1:2495137/1‑100 (MQ=255)
                            tGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActctc                                                                         >  2:769689/1‑100 (MQ=255)
                              ttCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACtctctc                                                                       >  2:1752467/1‑100 (MQ=255)
                               tCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTcc                                                                      >  1:1882714/1‑100 (MQ=255)
                                 aGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGc                                                                    >  1:408048/1‑100 (MQ=255)
                                    ttCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTAc                                                                 <  2:1469728/100‑1 (MQ=255)
                                      cGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGc                                                               <  2:857476/100‑1 (MQ=255)
                                       gCCAGTTTGTTCGTTAAGCAGCAATAATTACCCGGGTTGTCACCCGGATCATAGTCACTTGATGTGACTAAGATCCGATTAATACTCTCTCCGCTACGCa                                                              >  2:2368832/1‑100 (MQ=255)
                                       gCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCa                                                              >  1:197373/1‑100 (MQ=255)
                                       gCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCa                                                              >  2:3842457/1‑100 (MQ=255)
                                       gCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCa                                                              >  2:1514943/1‑100 (MQ=255)
                                        ccAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAg                                                             >  1:938407/1‑100 (MQ=255)
                                          aGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAgtg                                                           >  1:3208428/1‑100 (MQ=255)
                                           gTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAgtgt                                                          >  1:1215324/1‑100 (MQ=255)
                                              tGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGtgttgt                                                       <  1:2118650/100‑1 (MQ=255)
                                               gTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTa                                                      >  1:2231201/1‑100 (MQ=255)
                                               gTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTa                                                      >  2:3548021/1‑100 (MQ=255)
                                                  cGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGAt                                                   >  2:3273370/1‑100 (MQ=255)
                                                    ttAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTCGATCa                                                 >  2:3026029/1‑100 (MQ=255)
                                                        gcagcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTg                                             >  1:560374/1‑100 (MQ=255)
                                                        gcagcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTg                                             <  2:1890294/100‑1 (MQ=255)
                                                         cagcaATAATTACCCCGGTTGTCACCCGGAGCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGAGCAATTg                                             >  1:3346524/1‑99 (MQ=255)
                                                          agcaATACTTACCCCGGTTGTCACCCGGATCTTAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTgcg                                           <  2:3209010/100‑1 (MQ=255)
                                                           gcaATAATTACCCCGGTTGTCACCCGGCTCATCGTCAATTGATGTGACTATGATCCGATTACTACTCTCTCCGCTACGGAGTGTTGTAGATCACTTgcgc                                          >  2:541253/1‑100 (MQ=255)
                                                           gcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTgcgc                                          >  1:3579949/1‑100 (MQ=255)
                                                             aataatTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCAc                                        >  2:2231310/1‑100 (MQ=255)
                                                             aataatTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCAc                                        <  1:3849329/100‑1 (MQ=255)
                                                              ataatTACCCCGGTTGTCACCCGGAACATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTACGATACGCAGTTTTGTAGATCAATTGCGCa                                         >  2:3094756/1‑98 (MQ=255)
                                                                    aCCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCAtt                                 >  1:3236162/1‑100 (MQ=255)
                                                                     ccccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTg                                >  1:1069413/1‑100 (MQ=255)
                                                                          gTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAaata                           <  2:2734118/100‑1 (MQ=255)
                                                                           ttGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAaataa                          <  1:305471/100‑1 (MQ=255)
                                                                           ttGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAaataa                          >  2:1827542/1‑100 (MQ=255)
                                                                            tGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAaataat                         <  2:2445850/100‑1 (MQ=255)
                                                                              tCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTa                       >  1:3587921/1‑100 (MQ=255)
                                                                              tCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTa                       >  2:1326457/1‑100 (MQ=255)
                                                                              tCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTa                       <  1:1917592/100‑1 (MQ=255)
                                                                               taatcGGAGCATAGTCACATTACGTGACTATGATCCGATTAATACTCTCTCTGCTACGCAGTGTTGTCGATCAATTTCGCACTATCATTGAAATAATTAc                      >  2:397603/5‑100 (MQ=255)
                                                                               cACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAc                      >  2:720400/1‑100 (MQ=255)
                                                                               cACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAc                      >  2:969016/1‑100 (MQ=255)
                                                                               cACCCGGACCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAc                      >  1:2733550/1‑100 (MQ=255)
                                                                                aCCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAcc                     >  1:3424173/1‑100 (MQ=255)
                                                                                  tcGGATCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTg                   >  1:1407698/2‑100 (MQ=38)
                                                                                   cGGATCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCAGCTATGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTg                   >  2:3056926/1‑99 (MQ=37)
                                                                                    ggATCATAGTCACTTGATGTGGCTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCt                 <  1:986843/100‑1 (MQ=255)
                                                                                    ggATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCt                 <  2:2483260/100‑1 (MQ=255)
                                                                                    ggATCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCt                 >  2:3859809/1‑100 (MQ=38)
                                                                                     gATCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTa                >  1:1025319/1‑100 (MQ=38)
                                                                                     gATCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTa                >  2:1251838/1‑100 (MQ=255)
                                                                                      aTCATAGTCACGTGATGTGTCTATGATCCGATTAATTCTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAcc                     <  2:2457186/94‑1 (MQ=255)
                                                                                      aTCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAg               >  2:1972573/1‑100 (MQ=38)
                                                                                       tCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGt              >  1:2141163/1‑100 (MQ=38)
                                                                                        cATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAg               >  2:2388970/1‑98 (MQ=39)
                                                                                        cATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTg             >  2:1649743/1‑100 (MQ=39)
                                                                                        cATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTg             >  1:2534783/1‑100 (MQ=39)
                                                                                        cATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTg             >  2:1870591/1‑100 (MQ=39)
                                                                                         aTAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGa            >  2:2136893/1‑100 (MQ=39)
                                                                                          tAGTCACTTGATGTGACTATGATCCGATTAATACTCTTTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGat           <  2:2669417/100‑1 (MQ=255)
                                                                                          tAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGat           >  1:826755/1‑100 (MQ=255)
                                                                                           aGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGatt          >  1:3915977/1‑100 (MQ=255)
                                                                                            gTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGatta         >  1:20729/1‑100 (MQ=255)
                                                                                             tCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGattat        <  2:3657406/100‑1 (MQ=255)
                                                                                               aCTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTAttt      >  2:1919507/1‑100 (MQ=255)
                                                                                                nttGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTc     >  1:2591825/2‑100 (MQ=255)
                                                                                                cTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTc     <  2:1540850/100‑1 (MQ=255)
                                                                                                  tGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCaa   >  1:287945/1‑100 (MQ=255)
                                                                                                  tGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCaa   <  2:167957/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
ATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAAC  >  NC_000913/3657502‑3657700

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: