Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,520,513 T→C 20.9% D64G (GAT→GGT)  hofQ ← DNA catabolic putative fimbrial transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,520,5130TC20.9% 371.7 / 95.8 195D64G (GAT→GGT) hofQDNA catabolic putative fimbrial transporter
Reads supporting (aligned to +/- strand):  ref base T (84/70);  new base C (13/28);  total (97/98)
Fisher's exact test for biased strand distribution p-value = 1.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 197 total aligned reads displayed.

TGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGG  >  NC_000913/3520414‑3520610
                                                                                                   |                                                                                                 
tGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACAt                                                                                                   >  1:600480/1‑100 (MQ=255)
   cAGGCAATGGACTGCACTGGGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACGACAGTTTGTAGTGCCTGCTTCCAGGGAACATca                                                                                                 <  1:2139196/99‑2 (MQ=255)
   cAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTg                                                                                                >  1:4046821/1‑100 (MQ=255)
     ggCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGtt                                                                                              >  2:3432195/1‑100 (MQ=255)
      gCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTa                                                                                             <  2:2503059/100‑1 (MQ=255)
       cAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAg                                                                                            <  1:1742604/100‑1 (MQ=255)
        aaTGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGa                                                                                           >  1:716640/1‑100 (MQ=255)
        aaTGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGa                                                                                           >  2:3582159/1‑100 (MQ=255)
           ggAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGt                                                                                        >  2:2198373/1‑100 (MQ=255)
           ggAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGt                                                                                        >  2:1928243/1‑100 (MQ=255)
            gAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTa                                                                                       <  1:2755707/100‑1 (MQ=255)
            gAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTa                                                                                       >  1:3542539/1‑100 (MQ=255)
            gAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTa                                                                                       <  1:1626994/100‑1 (MQ=255)
              aTGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAAc                                                                                     <  1:3524082/100‑1 (MQ=255)
              aTGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAAc                                                                                     <  2:2626869/100‑1 (MQ=255)
                gCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGa                                                                                   >  2:3649872/1‑100 (MQ=255)
                 cACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGac                                                                                  >  2:416316/1‑100 (MQ=255)
                 cACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGac                                                                                  >  1:2237496/1‑100 (MQ=255)
                 cACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGac                                                                                  >  1:3904801/1‑100 (MQ=255)
                 cACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGac                                                                                  >  2:1726696/1‑100 (MQ=255)
                  aCTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTCCCAGGGAACATCTGTTAGATGTAACGaca                                                                                 >  1:212986/1‑100 (MQ=255)
                   cTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGacac                                                                                >  1:3769023/1‑100 (MQ=255)
                     gagagaATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCg                                                                              >  1:3917972/1‑100 (MQ=255)
                     gagagaATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCg                                                                              >  1:78619/1‑100 (MQ=255)
                     gagagaAGGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCg                                                                              >  2:2313067/1‑100 (MQ=255)
                       gagaATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTa                                                                            >  1:3805862/1‑100 (MQ=255)
                        agaATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTAc                                                                           <  2:3267417/100‑1 (MQ=255)
                         gaATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTAcc                                                                          <  2:3856344/100‑1 (MQ=255)
                         aaaTGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTAcc                                                                          >  1:609494/2‑100 (MQ=255)
                          aaTGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACTACACCGTACCg                                                                         <  1:3601842/100‑1 (MQ=255)
                           aTGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGc                                                                        <  1:1861830/100‑1 (MQ=255)
                           aTGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGc                                                                        <  1:939630/100‑1 (MQ=255)
                            tGTAGCCATCCTGCCGCGTTATGAGCACGGCGCTTTTCACTACAGTTTCTATTGTCTGCTTCCCGGGAACATCGGTTAGATGTAACGACAGCTTACCGc                                                                        >  1:1440298/1‑99 (MQ=255)
                             gTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTg                                                                      >  2:580746/1‑100 (MQ=255)
                             gTTGCCTTCCAGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTg                                                                      <  1:2888155/100‑1 (MQ=255)
                              ttGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGa                                                                     <  1:4101679/100‑1 (MQ=255)
                               tGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGAc                                                                    <  1:2278996/100‑1 (MQ=255)
                               tGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGAc                                                                    >  2:1349696/1‑100 (MQ=255)
                               tGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGAc                                                                    >  2:750982/1‑100 (MQ=255)
                                gCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGAc                                                                    >  2:3038947/1‑99 (MQ=255)
                                  cTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTc                                                                 <  2:2834809/100‑1 (MQ=255)
                                    tCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTg                                                               <  2:3756858/100‑1 (MQ=255)
                                     ccTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTgg                                                              <  1:1399368/100‑1 (MQ=255)
                                       tGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCg                                                            <  2:3333946/100‑1 (MQ=255)
                                        gCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGa                                                           >  2:2137966/1‑100 (MQ=255)
                                        gCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGa                                                           >  2:905641/1‑100 (MQ=255)
                                        gCCGCGTTATCAGTCAGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGa                                                           >  2:319809/1‑100 (MQ=255)
                                           gcgTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGacac                                                        <  2:2031412/100‑1 (MQ=255)
                                           gcgTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGacac                                                        >  2:1511058/1‑100 (MQ=255)
                                              ttATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGAc                                                     >  2:590234/1‑100 (MQ=255)
                                               tATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGAcc                                                    >  2:617275/1‑100 (MQ=255)
                                                aTCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGACGTAACGTCACCGTACCGCTGACGTCTGGCGACACGACCa                                                   >  1:3091813/1‑100 (MQ=255)
                                                aTCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCa                                                   <  1:1852953/100‑1 (MQ=255)
                                                  cAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAgg                                                 <  2:1010243/100‑1 (MQ=255)
                                                   aGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAgg                                                 >  1:3881987/1‑99 (MQ=255)
                                                    gTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGtt                                               >  2:1896762/1‑100 (MQ=255)
                                                        ggCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAAc                                           <  2:4240253/100‑1 (MQ=255)
                                                          cgcTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCaagatcggaag                                         >  2:2703947/1‑90 (MQ=255)
                                                           gcTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTtc                                        <  2:3658077/100‑1 (MQ=255)
                                                              tttCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTcc                                     <  1:319808/100‑1 (MQ=255)
                                                              tttCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTcc                                     <  1:130245/100‑1 (MQ=255)
                                                                tCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCACTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTg                                   >  1:1680479/1‑100 (MQ=255)
                                                                tCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTg                                   <  2:1677170/100‑1 (MQ=255)
                                                                   ntACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTc                                >  2:854897/2‑100 (MQ=255)
                                                                   cTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTc                                >  2:1129329/1‑100 (MQ=255)
                                                                   cTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTc                                >  1:340663/1‑100 (MQ=255)
                                                                     aCAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTACCGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTcag                              >  2:4211741/1‑100 (MQ=255)
                                                                     aCAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTcag                              >  2:2962875/1‑100 (MQ=255)
                                                                       aGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTcagcc                            <  2:932625/100‑1 (MQ=255)
                                                                        gTTTGTAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTcagcca                           >  2:4355958/1‑100 (MQ=255)
                                                                         tttGTAGTGCCTGCTTGCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTcagccag                          <  2:1597870/100‑1 (MQ=255)
                                                                         tttGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTcagccag                          >  2:1190211/1‑100 (MQ=255)
                                                                          ttGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTcagccagc                         >  2:2838956/1‑100 (MQ=255)
                                                                           tGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAgcg                        <  2:1030384/100‑1 (MQ=255)
                                                                            gTAGTGCCTGCTTTCAGGGAACACCTTTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAgcgc                       <  1:1585739/100‑1 (MQ=255)
                                                                             tAGTGCCTGCTTCCAGGGAACATCTGTTAGCGGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGcc                      <  1:4351537/100‑1 (MQ=255)
                                                                             tAGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGcc                      <  2:2129380/100‑1 (MQ=255)
                                                                              aGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCt                     >  1:1989811/1‑100 (MQ=255)
                                                                              aGTGCCTGCTTCCAGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCt                     <  1:1726689/100‑1 (MQ=255)
                                                                               gTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTg                    >  1:3960739/1‑100 (MQ=255)
                                                                               gTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTg                    >  2:2423106/1‑100 (MQ=255)
                                                                                tGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGCCGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGc                   >  1:2020081/1‑100 (MQ=255)
                                                                                tGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGc                   >  2:1751703/1‑100 (MQ=255)
                                                                                tGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGc                   <  1:231135/100‑1 (MQ=255)
                                                                                tGCCTGCTTCCAGGGAACATCTCTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGc                   >  1:2553419/1‑100 (MQ=255)
                                                                                 gCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCa                  <  1:1087896/100‑1 (MQ=255)
                                                                                 gCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCa                  <  2:1397122/100‑1 (MQ=255)
                                                                                   cTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAac                >  1:3698823/1‑100 (MQ=255)
                                                                                     gCTTCCAGGGAACATCTGTTAGATGTAACGACACTGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAacac              <  1:3986715/100‑1 (MQ=255)
                                                                                         ccAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTg           >  1:2075111/1‑99 (MQ=255)
                                                                                         ccAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGa          <  2:4046851/100‑1 (MQ=255)
                                                                                          cAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAg         >  1:3123300/1‑100 (MQ=255)
                                                                                          cAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAg         <  1:1537129/100‑1 (MQ=255)
                                                                                          cAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAg         >  2:912630/1‑100 (MQ=255)
                                                                                           aGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGc        >  2:2845100/1‑100 (MQ=255)
                                                                                           aGGGAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGc        <  1:2892373/100‑1 (MQ=255)
                                                                                            gggAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCt       >  2:2620144/1‑100 (MQ=255)
                                                                                            gggAACACCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCt       >  2:2228806/1‑100 (MQ=255)
                                                                                                aCATCTGTGAGAGGTAACGACAGCGTACCGCTGAGGTCTGGCGACACGACCAGGTTCACCTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTa      <  1:1389545/97‑1 (MQ=255)
                                                                                                 cATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCgg  >  2:1830974/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
TGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGG  >  NC_000913/3520414‑3520610

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: