Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,657,603 T→A 17.6% 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,6030TA17.6% 403.5 / 73.7 167intergenic (+36/‑763)hdeD/gadEacid‑resistance membrane protein/gad regulon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (74/62);  new base A (14/15);  total (90/77)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

Only 100 of 178 total aligned reads displayed.

GAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAACCT  >  NC_000913/3657505‑3657702
                                                                                                  |                                                                                                   
gAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtg                                                                                                    >  2:3207463/1‑100 (MQ=255)
gAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtg                                                                                                    >  2:1850091/1‑100 (MQ=255)
 aaCTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtgt                                                                                                   >  2:3356901/1‑100 (MQ=255)
 aaCTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtgt                                                                                                   <  2:3922345/100‑1 (MQ=255)
  aCTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtgtg                                                                                                  >  1:1886684/1‑100 (MQ=255)
    tGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGAc                                                                                                >  1:201727/1‑100 (MQ=255)
     gATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACt                                                                                               >  2:3227382/1‑100 (MQ=255)
      atatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTa                                                                                              <  2:384627/100‑1 (MQ=255)
      atatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTa                                                                                              <  2:3711914/100‑1 (MQ=255)
       tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTaa                                                                                             >  1:3243137/1‑99 (MQ=255)
       tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTAt                                                                                             <  2:3108958/100‑1 (MQ=255)
       tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTAt                                                                                             >  2:3015456/1‑100 (MQ=255)
       tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGACCATAGTCACATCAAGTGACTAt                                                                                             <  1:2075537/100‑1 (MQ=255)
        atTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATg                                                                                            <  1:2764409/100‑1 (MQ=39)
         tttAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGa                                                                                           <  1:3792101/100‑1 (MQ=38)
         tttAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGa                                                                                           <  1:3908202/100‑1 (MQ=38)
           tAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATc                                                                                         <  1:2123447/100‑1 (MQ=255)
           tAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATc                                                                                         <  2:3224128/100‑1 (MQ=255)
             gcgcCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCg                                                                                       >  2:3794665/1‑100 (MQ=255)
             gcgcCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGATCCg                                                                                       <  2:1896772/100‑1 (MQ=38)
             gcgcCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGATCCg                                                                                       <  1:3116063/100‑1 (MQ=38)
              cgccgcCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGa                                                                                      >  2:2180831/1‑100 (MQ=255)
              cgccgcCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGa                                                                                      >  1:3163157/1‑100 (MQ=255)
               gccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGAt                                                                                     >  1:2322125/1‑100 (MQ=255)
               gccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGATCCggg                                                                                     <  2:1617418/100‑3 (MQ=38)
                ccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGAtt                                                                                    <  2:1114720/100‑1 (MQ=255)
                ccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGAtt                                                                                    >  1:3657434/1‑100 (MQ=255)
                 cgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGAt                                                                                     <  2:2162281/98‑1 (MQ=255)
                 cgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTa                                                                                   <  1:2958925/100‑1 (MQ=255)
                 cgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGATCCgggtg                                                                                   <  1:1302863/100‑5 (MQ=37)
                  gccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTaa                                                                                  >  2:2907229/1‑100 (MQ=255)
                  gccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTaa                                                                                  >  1:1343663/1‑100 (MQ=255)
                     aGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATAc                                                                               <  2:2277007/100‑1 (MQ=255)
                      gCCTTTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGTAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATAct                                                                              <  2:62126/100‑1 (MQ=255)
                      gCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATAct                                                                              <  2:3072086/100‑1 (MQ=255)
                       ccTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActc                                                                             >  2:363017/1‑100 (MQ=255)
                         tGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATActctc                                                                           >  2:769689/1‑100 (MQ=255)
                              aGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGc                                                                      >  1:2415202/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)
                                    gCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCa                                                                >  1:197373/1‑100 (MQ=255)
                                       aGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAgtg                                                             >  1:3208428/1‑100 (MQ=255)
                                        gTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAgtgt                                                            >  1:1215324/1‑100 (MQ=255)
                                        gTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAgtgt                                                            <  2:197683/100‑1 (MQ=255)
                                           tGTTCGTTAAGCAGCAATAATTACCCCTGTTGTCACCCGGCTCATAGTCACTTGATGTGACTATGATCCGCTTCATACTACCTCCGCTACGCAGtgttgt                                                         >  2:436934/1‑100 (MQ=255)
                                            gTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTa                                                        >  1:1363847/1‑100 (MQ=255)
                                            gTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTa                                                        >  2:3548021/1‑100 (MQ=255)
                                               cGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGAt                                                     >  2:3273370/1‑100 (MQ=255)
                                               cGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGAt                                                     >  2:305471/1‑100 (MQ=255)
                                               cGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGAt                                                     >  2:166826/1‑100 (MQ=255)
                                                 ttAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCa                                                   >  2:1787477/1‑100 (MQ=255)
                                                   aaGCAGCAATAATTACCCCGGTTGTCACCCGGAACATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAAt                                                 >  2:1512699/1‑100 (MQ=255)
                                                     gcagcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTg                                               >  1:560374/1‑100 (MQ=255)
                                                      cagcaATAATTACCCCGGTTGTCACCCGGAGCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGAGCAATTg                                               >  1:3346524/1‑99 (MQ=255)
                                                        gcaATAATTACCCCGGTTGTCACCCGGCTCATCGTCAATTGATGTGACTATGATCCGATTACTACTCTCTCCGCTACGGAGTGTTGTAGATCACTTgcgc                                            >  2:541253/1‑100 (MQ=255)
                                                        gcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTgcgc                                            <  2:1215315/100‑1 (MQ=255)
                                                        gcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTgcgc                                            >  1:3579949/1‑100 (MQ=255)
                                                         caATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGCTTAATACTCTCTCCGCTCCGCAGTGTTGTAGATAAATTgcg                                             >  2:2479147/1‑98 (MQ=255)
                                                          aataatTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCAc                                          <  1:3849329/100‑1 (MQ=255)
                                                          aataatTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCAc                                          >  2:2231310/1‑100 (MQ=255)
                                                           ataatTACCCCGGTTGTCACCCGGAACATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTACGATACGCAGTTTTGTAGATCAATTGCGCa                                           >  2:3094756/1‑98 (MQ=255)
                                                            taatTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTa                                        <  2:1118803/100‑1 (MQ=255)
                                                              atTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATc                                      <  1:64523/100‑1 (MQ=255)
                                                                 aCCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCAtt                                   >  1:2495138/1‑100 (MQ=255)
                                                                 aCCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCAtt                                   >  1:3236162/1‑100 (MQ=255)
                                                                 aCCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATAAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCAt                                    <  1:2629608/99‑1 (MQ=255)
                                                                  ccccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTg                                  >  1:1069413/1‑100 (MQ=255)
                                                                  ccccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTg                                  <  1:1980462/100‑1 (MQ=255)
                                                                    ccGGTTGTCACCCGTATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCTCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGaa                                <  2:2427360/100‑1 (MQ=255)
                                                                    ccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGaa                                <  2:2085913/100‑1 (MQ=255)
                                                                    ccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGaa                                <  1:985748/100‑1 (MQ=255)
                                                                       gTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAaata                             <  2:2734118/100‑1 (MQ=255)
                                                                        ttGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAaataa                            <  1:305471/100‑1 (MQ=255)
                                                                         tGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAaataat                           <  2:2445850/100‑1 (MQ=255)
                                                                          gTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAAtt                          >  1:3586602/1‑100 (MQ=255)
                                                                           tCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACACAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAAtt                          <  2:1882386/99‑1 (MQ=255)
                                                                            cACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAc                        >  2:969016/1‑100 (MQ=255)
                                                                            cACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAc                        >  2:720400/1‑100 (MQ=255)
                                                                            cACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAc                        >  2:2634912/1‑100 (MQ=255)
                                                                             aCCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAcc                       <  1:769693/100‑1 (MQ=255)
                                                                             aCCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAcc                       >  1:3424173/1‑100 (MQ=255)
                                                                              cccGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCt                      <  2:152104/100‑1 (MQ=255)
                                                                               tcGGATCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTg                     >  1:1407698/2‑100 (MQ=38)
                                                                                 ggATCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCt                   >  2:3859809/1‑100 (MQ=38)
                                                                                  gATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTa                  <  1:2231322/100‑1 (MQ=255)
                                                                                  gATCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTa                  >  1:1025319/1‑100 (MQ=38)
                                                                                   aTCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAg                 >  2:126418/1‑100 (MQ=255)
                                                                                    tCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGt                >  1:2141163/1‑100 (MQ=38)
                                                                                     cATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTg               >  1:1773506/1‑100 (MQ=39)
                                                                                     cATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTg               >  1:2534783/1‑100 (MQ=39)
                                                                                        aGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGatt            >  1:3915977/1‑100 (MQ=255)
                                                                                          tCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGattat          <  2:3657406/100‑1 (MQ=255)
                                                                                            aCTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTAttt        >  2:1919507/1‑100 (MQ=255)
                                                                                              ttGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGAACAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCa      <  2:2686197/100‑1 (MQ=255)
                                                                                               tGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCaa     <  2:167957/100‑1 (MQ=255)
                                                                                               tGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCaa     >  2:2381677/1‑100 (MQ=255)
                                                                                                gATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAAc    <  2:1346647/100‑1 (MQ=255)
                                                                                                gATGTGACTATGATCCGATTAATACTCTCTCAGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAAc    <  1:933151/100‑1 (MQ=255)
                                                                                                 aTGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCAAcc   <  1:466245/100‑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: