Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,369,184 T→C C315C (TGT→TGC ycjM → putative polysaccharide hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,369,1840TC98.3% 903.1 / 2.7 288C315C (TGT→TGCycjMputative polysaccharide hydrolase
Reads supporting (aligned to +/- strand):  ref base T (4/1);  new base C (146/137);  total (150/138)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-02

Only 100 of 310 total aligned reads displayed.

AGGCGATGACGAAGCACATATGGTGTACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGTGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGATGGCATCGGGCTAAACCCGC  >  NC_000913/1369094‑1369272
                                                                                          |                                                                                        
aGGCGATGACGAAGCACATATGGTGTACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTgcgc                                                                                        >  1:5576688/1‑93 (MQ=255)
aGGCGATGACGAAGCACATATGGTGTACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTgcgc                                                                                        >  1:5265336/1‑93 (MQ=255)
aGGCGATGACGAAGCACATATGGTGTACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTgcgc                                                                                        >  2:9159259/1‑93 (MQ=255)
aGGCGATGACGAAGCACATATGGTGTACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTgcgc                                                                                        >  1:21909/1‑93 (MQ=255)
 ggCGATGACGAAGCACATATGGTGTACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGcgcg                                                                                       <  2:8080304/93‑1 (MQ=255)
 ggCGATGACGAAGCACATATGGTGTACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGcgcg                                                                                       <  1:10011251/93‑1 (MQ=255)
    gATGACGAAGCACATATGGTGTACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTgg                                                                                    >  2:4609230/1‑93 (MQ=255)
    gATGACGAAGCACATATGGTGTACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTgg                                                                                    >  2:5285918/1‑93 (MQ=255)
    gATGACGAAGCACATATGGTGTACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTgg                                                                                    >  1:3824232/1‑93 (MQ=255)
     aTGACGAAGCACATCTGGTGGACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTggg                                                                                   <  1:5123211/93‑1 (MQ=255)
     aTGACGAAGCACATATGGTGTACCAGTTCTCGCTTCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTggg                                                                                   >  2:7090816/1‑93 (MQ=255)
     aTGACGAAGCACATATGGTGTACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTggg                                                                                   <  1:6784360/93‑1 (MQ=255)
       gACGAAGCACATATGGTGTACCCGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGAGCTTTGACCGTGGgca                                                                                 >  1:9860279/1‑92 (MQ=255)
        aCGAAGCACATATGGTGTACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGgcgc                                                                                >  1:7442404/1‑93 (MQ=255)
        aCGAAGCACATATGGTGTACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGgcgc                                                                                >  1:5749921/1‑93 (MQ=255)
         cGAAACACATATGGTGTACCAGTTCTCGCTCCCGCCACTGGTGCTGAATGCGGTGCAAAAACAGAACGTTGAGTCGCTTTGCGAGTGGGCGCa                                                                               <  2:5426232/93‑1 (MQ=255)
           aaGCACATATGGTGTACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAACAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCaaa                                                                             >  1:8289641/1‑93 (MQ=255)
            aGCACATATGGTGTACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCaaaa                                                                            <  1:9585327/93‑1 (MQ=255)
               acaTATGGTGTACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCt                                                                         <  2:1288401/93‑1 (MQ=255)
                caTATGGTGTACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTg                                                                        <  1:425557/93‑1 (MQ=255)
                caTATGGTGTACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTg                                                                        <  1:9176044/93‑1 (MQ=255)
                 atatGGTGTACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGa                                                                       <  1:6592472/93‑1 (MQ=255)
                  tatGGTGTACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAATAACAGAACGTTGAGGCGCTTTGCGCGAGGGCGCAAAACCTGac                                                                      >  2:4944071/1‑93 (MQ=255)
                  tatGGTGTACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGac                                                                      >  2:2595734/1‑93 (MQ=255)
                  tatGGTGTACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGac                                                                      >  1:2568906/1‑93 (MQ=255)
                  tatGGTGTACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGac                                                                      >  1:6312467/1‑93 (MQ=255)
                    tGGTGTACCAGTTTTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACTTGacac                                                                    >  2:8976517/1‑93 (MQ=255)
                     ggTGTACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGTGCGTGGGCGCAAAACCTGACACt                                                                   <  2:6959935/93‑1 (MQ=255)
                        gtACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTAcc                                                                >  2:8011390/1‑93 (MQ=255)
                        gtACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTAcc                                                                <  2:6398185/93‑1 (MQ=255)
                        gtACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGCGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACAGTAca                                                                >  1:4740153/1‑92 (MQ=255)
                         tACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCt                                                               <  1:2054199/93‑1 (MQ=255)
                           ccAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTc                                                             <  2:4299893/93‑1 (MQ=255)
                            cAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTcc                                                            >  1:5996136/1‑93 (MQ=255)
                            cAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTcc                                                            <  2:8609451/93‑1 (MQ=255)
                              gTTCTCGCTGACGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAg                                                          >  2:6547536/1‑93 (MQ=255)
                               ttCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGc                                                         >  2:7917972/1‑93 (MQ=255)
                               ttCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGc                                                         >  1:9632722/1‑93 (MQ=255)
                               ttCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGc                                                         >  2:8421651/1‑93 (MQ=255)
                               ttCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGc                                                         >  2:5565165/1‑93 (MQ=255)
                                 ctcGCTGCCGCCGCTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCaa                                                       <  2:3571682/93‑1 (MQ=255)
                                  tcGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCCac                                                      >  1:713056/1‑93 (MQ=255)
                                  tcGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAac                                                      >  1:4161904/1‑93 (MQ=255)
                                   cGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAaca                                                     <  1:618630/93‑1 (MQ=255)
                                   cGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAaca                                                     <  2:7808809/93‑1 (MQ=255)
                                     cTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAAcacc                                                   >  1:1707945/1‑93 (MQ=255)
                                     cTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAAcacc                                                   >  2:5243967/1‑93 (MQ=255)
                                       gccgccACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAAcaccac                                                 >  2:7452227/1‑93 (MQ=255)
                                        ccgccACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACaccacc                                                <  1:7907875/93‑1 (MQ=255)
                                        ccgccACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACaccacc                                                <  2:3486544/93‑1 (MQ=255)
                                        ccgccACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACaccacc                                                >  2:248632/1‑93 (MQ=255)
                                            cACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGt                                            >  2:3834960/1‑93 (MQ=255)
                                             aCTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGtt                                           >  1:8804236/1‑93 (MQ=255)
                                                ggTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTaa                                        >  1:8970776/1‑93 (MQ=255)
                                                ggTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTaa                                        <  1:6879766/93‑1 (MQ=255)
                                                     tgcATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTcc                                   <  1:9217504/93‑1 (MQ=255)
                                                     tgcATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTcc                                   <  1:9434694/93‑1 (MQ=255)
                                                     tgcATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTcc                                   >  2:3027666/1‑93 (MQ=255)
                                                      gcATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCt                                  >  2:9652475/1‑93 (MQ=255)
                                                      gcATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCt                                  <  2:6000407/93‑1 (MQ=255)
                                                      gcATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCt                                  <  1:1209931/93‑1 (MQ=255)
                                                      gcATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCt                                  >  1:779105/1‑93 (MQ=255)
                                                       gaTGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTc                                 <  1:8194710/92‑1 (MQ=255)
                                                       cATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGTGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTc                                 >  1:5669680/1‑93 (MQ=255)
                                                       cATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTc                                 >  2:1628524/1‑93 (MQ=255)
                                                         tGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGc                               <  2:4780669/93‑1 (MQ=255)
                                                          gCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGcc                              >  1:281514/1‑93 (MQ=255)
                                                          gCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGcc                              <  2:3378006/93‑1 (MQ=255)
                                                          gCGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGcc                              >  2:4918233/1‑93 (MQ=255)
                                                           cGGTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCct                             <  2:6177042/93‑1 (MQ=255)
                                                            ggTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCctc                            <  2:3316791/93‑1 (MQ=255)
                                                            ggTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCctc                            <  1:7300567/93‑1 (MQ=255)
                                                             gTGCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCctct                           <  2:1869000/93‑1 (MQ=255)
                                                               gCAAAAACAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCa                         <  2:7206626/93‑1 (MQ=255)
                                                                  aaaaCAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGa                      >  1:6455311/1‑93 (MQ=255)
                                                                   aaaCAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGAt                     >  1:7695873/1‑93 (MQ=255)
                                                                   aaaCAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGAt                     >  1:6154269/1‑93 (MQ=255)
                                                                   aaaCAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGAt                     <  2:2208063/93‑1 (MQ=255)
                                                                    aaCAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGATg                    <  2:8297030/93‑1 (MQ=255)
                                                                    aaCAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGATg                    >  1:4653351/1‑93 (MQ=255)
                                                                    aaCAGAACGTGGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGATg                    <  2:4039281/93‑1 (MQ=255)
                                                                     aCAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGATgg                   <  1:4447907/93‑1 (MQ=255)
                                                                     aCAGAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGATgg                   >  1:5435562/1‑93 (MQ=255)
                                                                        gAACGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGATGGCAt                <  1:1326143/93‑1 (MQ=255)
                                                                         aaCGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGATGGCATc               <  2:8820695/93‑1 (MQ=255)
                                                                         aaCGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGATGGCATc               <  2:1541993/93‑1 (MQ=255)
                                                                         aaCGTTGAGGCGCTTTGCGCGGGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGATGGCATc               <  1:8743848/93‑1 (MQ=255)
                                                                           cGTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGATGGCATCgg             <  1:8474823/93‑1 (MQ=255)
                                                                            gTTGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGATGGCATCggg            >  1:3630601/1‑93 (MQ=255)
                                                                             ttGAGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGATGGCATCGGGc           <  1:904481/93‑1 (MQ=255)
                                                                                ggtcgcTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGATGGCATCGGGCTaa        <  1:3943357/90‑1 (MQ=255)
                                                                                aGGCGCTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGATGGCATCGGGCTaa        <  2:9047691/93‑1 (MQ=255)
                                                                                 ggggcTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGATGGCATCGGGCTaaa       <  2:2448654/90‑1 (MQ=255)
                                                                                   cgcTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGATGGCATCGGGCTAAAcc     >  1:2891671/1‑93 (MQ=255)
                                                                                   cgcTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGATGGCATCGGGCTAAAcc     <  1:2387425/93‑1 (MQ=255)
                                                                                   cgcTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGATGGCATCGGGCTAAAcc     <  2:9829354/93‑1 (MQ=255)
                                                                                    gcTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGATGGCATCGGGCTAAAccc    <  1:798321/93‑1 (MQ=255)
                                                                                    ccTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGATGGCATCGGGCTAAAccc    <  2:9261026/92‑1 (MQ=255)
                                                                                     cTTTGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGATGGCATCGGGCTAAACCCg   >  1:3529215/1‑93 (MQ=255)
                                                                                      tttGCGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGATGGCATCGGGCTAAACCCGc  <  2:440836/93‑1 (MQ=255)
                                                                                          |                                                                                        
AGGCGATGACGAAGCACATATGGTGTACCAGTTCTCGCTGCCGCCACTGGTGCTGCATGCGGTGCAAAAACAGAACGTTGAGGCGCTTTGTGCGTGGGCGCAAAACCTGACACTACCTTCCAGCAACACCACCTGGTTTAACTTCCTCGCCTCTCACGATGGCATCGGGCTAAACCCGC  >  NC_000913/1369094‑1369272

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