Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,657,599 T→A 17.4% 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,5990TA17.4% 413.9 / 78.2 173intergenic (+32/‑767)hdeD/gadEacid‑resistance membrane protein/gad regulon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (80/62);  new base A (15/15);  total (95/78)
Fisher's exact test for biased strand distribution p-value = 5.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

Only 100 of 179 total aligned reads displayed.

GAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCA  >  NC_000913/3657500‑3657698
                                                                                                   |                                                                                                   
gAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACt                                                                                                     >  2:3023025/1‑100 (MQ=255)
gAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACt                                                                                                     >  1:2761589/1‑100 (MQ=255)
gAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACt                                                                                                     >  1:167957/1‑100 (MQ=255)
  aTTGAACGGATATTTAGCCCCGCCAGCCGGTTCAGCTTCGCCAGTTTGTTCGTTAAGCATCAATAATTATCCCGGTTGTCACCCGTATCATAGTCACt                                                                                                     <  2:2432834/98‑1 (MQ=255)
   ttGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGa                                                                                                  >  1:152104/1‑100 (MQ=255)
    tGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTCCCCCGGTTGTCACCAGGATCATAGTCACTTGa                                                                                                  <  2:453387/99‑1 (MQ=255)
     gAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtg                                                                                                >  2:3207463/1‑100 (MQ=255)
     gAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtg                                                                                                <  2:2275073/100‑1 (MQ=255)
      aaCTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGAtgt                                                                                               >  2:3356901/1‑100 (MQ=255)
      aaCTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGAGCATAGTCACTTGAtgt                                                                                               >  1:3095215/1‑100 (MQ=255)
         tGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGAc                                                                                            <  1:3162667/100‑1 (MQ=255)
         tGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGAc                                                                                            >  1:201727/1‑100 (MQ=255)
          gATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACt                                                                                           >  1:3490382/1‑100 (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)
            tatTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGACCATAGTCACATCAAGTGACTAt                                                                                         <  1:2075537/100‑1 (MQ=255)
             atTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATg                                                                                        <  1:2764409/100‑1 (MQ=39)
              tttAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGa                                                                                       <  1:768875/100‑1 (MQ=38)
              tttAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGa                                                                                       <  1:3908202/100‑1 (MQ=38)
                tAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATc                                                                                     >  1:471211/1‑100 (MQ=255)
                tAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATc                                                                                     <  2:2027080/100‑1 (MQ=255)
                  gcgcCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCg                                                                                   <  2:3873770/100‑1 (MQ=255)
                  gcgcCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCg                                                                                   <  1:3327585/100‑1 (MQ=255)
                  gcgcCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGATCCg                                                                                   <  2:1896772/100‑1 (MQ=38)
                    gccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGAt                                                                                 >  1:2322125/1‑100 (MQ=255)
                     ccgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGAtt                                                                                <  2:1114720/100‑1 (MQ=255)
                      cgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGAt                                                                                 <  2:2162281/98‑1 (MQ=255)
                      cgccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACATCAAGTGACTATGATCCgggtg                                                                               <  1:1302863/100‑5 (MQ=37)
                       gccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTaa                                                                              >  1:1343663/1‑100 (MQ=255)
                        ccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAAt                                                                             <  2:1179745/100‑1 (MQ=255)
                        ccAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAAt                                                                             <  2:3593855/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:2415202/1‑100 (MQ=255)
                                   aGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGc                                                                  >  1:408048/1‑100 (MQ=255)
                                        cGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGc                                                             <  2:857476/100‑1 (MQ=255)
                                         gCCAGTTTGTTCGTTAAGCAGCAATAATTACCCGGGTTGTCACCCGGATCATAGTCACTTGATGTGACTAAGATCCGATTAATACTCTCTCCGCTACGCa                                                            >  2:2368832/1‑100 (MQ=255)
                                         gCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCa                                                            >  2:3842457/1‑100 (MQ=255)
                                            aGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAgtg                                                         >  1:3208428/1‑100 (MQ=255)
                                            aGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAgtg                                                         >  2:301373/1‑100 (MQ=255)
                                             gTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAgtgt                                                        >  1:1215324/1‑100 (MQ=255)
                                                tGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGtgttgt                                                     <  1:2118650/100‑1 (MQ=255)
                                                 gTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTa                                                    >  2:3548021/1‑100 (MQ=255)
                                                 gTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTa                                                    >  1:1363847/1‑100 (MQ=255)
                                                    cGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGAt                                                 >  2:3273370/1‑100 (MQ=255)
                                                    cGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGAt                                                 >  2:305471/1‑100 (MQ=255)
                                                        aaGCAGCAATAATTACCCCGGTTGTCACCCGGAACATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAAt                                             >  2:1512699/1‑100 (MQ=255)
                                                            agcaATACTTACCCCGGTTGTCACCCGGATCTTAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTgcg                                         <  2:3209010/100‑1 (MQ=255)
                                                             gcaATAATTACCCCGGTTGTCACCCGGCTCATCGTCAATTGATGTGACTATGATCCGATTACTACTCTCTCCGCTACGGAGTGTTGTAGATCACTTgcgc                                        >  2:541253/1‑100 (MQ=255)
                                                             gcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTgcgc                                        <  2:1215315/100‑1 (MQ=255)
                                                             gcaATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTgcgc                                        >  1:3579949/1‑100 (MQ=255)
                                                                ataatTACCCCGGTTGTCACCCGGAACATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTACGATACGCAGTTTTGTAGATCAATTGCGCa                                       >  2:3094756/1‑98 (MQ=255)
                                                                   atTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATc                                  <  1:64523/100‑1 (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)
                                                                       ccccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTg                              >  1:3385480/1‑100 (MQ=255)
                                                                       ccccGGTTGTCACCCGGATCATAGTCACTTGAGGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTg                              >  1:3741759/1‑100 (MQ=255)
                                                                         ccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGaa                            <  2:2085913/100‑1 (MQ=255)
                                                                         ccGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGaa                            <  1:985748/100‑1 (MQ=255)
                                                                             ttGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAaataa                        >  2:1827542/1‑100 (MQ=255)
                                                                             ttGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAaataa                        <  1:305471/100‑1 (MQ=255)
                                                                                tCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTa                     <  1:1917592/100‑1 (MQ=255)
                                                                                tCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTa                     >  1:3587921/1‑100 (MQ=255)
                                                                                tCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACACAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAAtt                      <  2:1882386/99‑1 (MQ=255)
                                                                                 taatcGGAGCATAGTCACATTACGTGACTATGATCCGATTAATACTCTCTCTGCTACGCAGTGTTGTCGATCAATTTCGCACTATCATTGAAATAATTAc                    >  2:397603/5‑100 (MQ=255)
                                                                                 cACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAc                    >  2:969016/1‑100 (MQ=255)
                                                                                 cACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAc                    >  2:720400/1‑100 (MQ=255)
                                                                                  aCCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTAcc                   <  1:3842486/100‑1 (MQ=255)
                                                                                   cccGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCt                  <  2:152104/100‑1 (MQ=255)
                                                                                   cccGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCt                  >  2:2181139/1‑100 (MQ=255)
                                                                                    tcGGATCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTg                 >  1:1407698/2‑100 (MQ=38)
                                                                                     cGGATCATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCAGCTATGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTg                 >  2:3056926/1‑99 (MQ=37)
                                                                                      ggATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCt               <  2:2483260/100‑1 (MQ=255)
                                                                                      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)
                                                                                        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           >  1:2534783/1‑100 (MQ=39)
                                                                                          cATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTg           >  2:1649743/1‑100 (MQ=39)
                                                                                          cATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTg           >  2:1870591/1‑100 (MQ=39)
                                                                                          cATAGTCACATCAAGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTg           >  1:1773506/1‑100 (MQ=39)
                                                                                           aTAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGa          <  1:3188185/100‑1 (MQ=255)
                                                                                            tAGTCACTTGATGTGACTATGATCCGATTAATACTCTTTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGat         <  2:2669417/100‑1 (MQ=255)
                                                                                            tAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGat         >  1:826755/1‑100 (MQ=255)
                                                                                             aGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGatt        >  1:3915977/1‑100 (MQ=255)
                                                                                               tCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGattat      <  2:3657406/100‑1 (MQ=255)
                                                                                                cACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGattatt     >  2:2867519/1‑100 (MQ=255)
                                                                                                  cTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTc   <  2:1540850/100‑1 (MQ=255)
                                                                                                   ttGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGAACAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCa  <  2:2686197/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GAATTGAACTGATATTTAGCGCCGCCAGCCTGTTCAGCTTCGCCAGTTTGTTCGTTAAGCAGCAATAATTACCCCGGTTGTCACCCGGATCATAGTCACTTGATGTGACTATGATCCGATTAATACTCTCTCCGCTACGCAGTGTTGTAGATCAATTGCGCACTATCATTGAAATAATTACCTGCTAGTGATTATTTCA  >  NC_000913/3657500‑3657698

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: