Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,088,091 G→C 17.8% intergenic (+34/+16) yiiG → / ← frvR DUF3829 family lipoprotein/putative frv operon regulator; contains a PTS EIIA domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,088,0910GC17.8% 363.1 / 77.2 160intergenic (+34/+16)yiiG/frvRDUF3829 family lipoprotein/putative frv operon regulator; contains a PTS EIIA domain
Reads supporting (aligned to +/- strand):  ref base G (56/74);  new base C (17/11);  total (75/85)
Fisher's exact test for biased strand distribution p-value = 9.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-05

Only 100 of 172 total aligned reads displayed.

CTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGC  >  NC_000913/4087994‑4088189
                                                                                                 |                                                                                                   
ctctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcggcg                                                                                                   <  1:2493515/100‑1 (MQ=255)
ctctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcggcg                                                                                                   <  2:3280095/100‑1 (MQ=255)
ctctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcggcg                                                                                                   <  1:1669570/100‑1 (MQ=255)
 tctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAGATGCAGTGTCGGATGCGgcgc                                                                                                  <  1:819626/100‑1 (MQ=255)
   tttCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCa                                                                                                 <  1:1638785/99‑1 (MQ=255)
        cGAGAGCATTATGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAATCGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGac                                                                                           <  2:390057/100‑1 (MQ=255)
        cGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGac                                                                                           <  1:647070/100‑1 (MQ=255)
        cGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGac                                                                                           <  2:3833717/100‑1 (MQ=255)
        cGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGac                                                                                           <  1:3454864/100‑1 (MQ=255)
        cGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGAGGCGGCGCATCCGac                                                                                           >  2:580166/1‑100 (MQ=255)
         gagaGTCTTGCGCGAGTACAGCGAAAGGGTTGATGACTATAATAGCCTGCGTTAACCGTTTTTGATGATAAATGCTGTGTTGGATGCGGCGCATCCGac                                                                                           <  2:2435975/99‑1 (MQ=255)
         gagaGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGaca                                                                                          <  2:3568440/100‑1 (MQ=255)
          agagCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGacac                                                                                         <  2:1312947/100‑1 (MQ=255)
          agagCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGacac                                                                                         <  1:2301564/100‑1 (MQ=255)
          agagCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGacac                                                                                         <  1:397298/100‑1 (MQ=255)
          agagCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGacac                                                                                         <  2:3047194/100‑1 (MQ=255)
             gCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGt                                                                                      >  2:1680967/1‑100 (MQ=255)
              cATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGt                                                                                      <  2:2031217/99‑1 (MQ=255)
                 tACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑c                                                                                 <  1:3781678/100‑1 (MQ=255)
                  aCGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑cc                                                                                <  1:1883979/100‑1 (MQ=255)
                  aCGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑cc                                                                                <  1:3277486/100‑1 (MQ=255)
                    gcgAGTACAACGAAATGGTTGATGACTATAATAGTCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑cccc                                                                              <  1:502729/100‑1 (MQ=255)
                       actACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTtg                                                                           >  1:2648040/3‑100 (MQ=255)
                          aCAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGa                                                                        <  2:3407680/100‑1 (MQ=255)
                           cAACGAAATGGTTGATGACTATAATAGGCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGaa                                                                       <  2:3260457/100‑1 (MQ=255)
                           cAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTTTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGaa                                                                       >  1:629813/1‑100 (MQ=255)
                           cAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGaa                                                                       <  1:481757/100‑1 (MQ=255)
                           cAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGAGTCGAAAGAGGCGCAACCGACACAGTTCA‑CCCCTTGTGaa                                                                       >  1:233504/1‑100 (MQ=255)
                            aaCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAg                                                                      <  2:934962/100‑1 (MQ=255)
                             aCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGt                                                                     <  1:2097703/100‑1 (MQ=255)
                             aCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGt                                                                     <  1:3850949/100‑1 (MQ=255)
                             aCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCACATCCGACACAGTTCA‑CCCCTTGTGAAGt                                                                     >  1:34374/1‑100 (MQ=255)
                              cGAAATGGTTTATAACTATAATAGCCTACGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCAGCCGACACAGTTCA‑CCCCTTGTGAAGTa                                                                    <  1:530106/100‑1 (MQ=255)
                               gAAATGATTGATGTATATAATAGCCTGCGTTAAACGTTTTAGATACTATATCCTGTGTCGGATGCGGCGCATCCAACACAGTTAT‑CCACTTTTGa                                                                        <  2:1998574/95‑1 (MQ=255)
                                    ggTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCa                                                              <  1:2028357/100‑1 (MQ=255)
                                    ggTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCa                                                              >  2:357089/1‑100 (MQ=255)
                                      ttaatgaCTATAATATCCTTCTTTAATCTTTTTTTTTTATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATa                                                            <  1:292874/97‑1 (MQ=255)
                                        gatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGTAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATaaa                                                          <  1:1105704/100‑1 (MQ=255)
                                        gatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATaaa                                                          >  1:35614/1‑100 (MQ=255)
                                         atgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAAt                                                         <  1:1707204/100‑1 (MQ=255)
                                          tgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATg                                                        <  1:357882/100‑1 (MQ=255)
                                           gaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGAATGCGACGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGc                                                       >  1:1026136/1‑100 (MQ=255)
                                            aCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGcc                                                      <  1:1068546/100‑1 (MQ=255)
                                            aCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGcc                                                      <  2:2804388/100‑1 (MQ=255)
                                             cTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCt                                                     >  2:1416904/1‑100 (MQ=255)
                                                  ataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCatat                                                <  2:1077322/100‑1 (MQ=255)
                                                   taGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATc                                               >  2:2866301/1‑100 (MQ=255)
                                                    aGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATcc                                              >  2:264746/1‑100 (MQ=255)
                                                     gCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCt                                             >  1:3037265/1‑100 (MQ=255)
                                                       cTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGt                                           >  2:3683950/1‑100 (MQ=255)
                                                        tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGtt                                          <  1:3788643/100‑1 (MQ=255)
                                                         gCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTa                                         >  1:3603016/1‑100 (MQ=255)
                                                              aaaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCg                                    <  2:539715/100‑1 (MQ=255)
                                                               aaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCgg                                   <  1:3104551/100‑1 (MQ=255)
                                                                aCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGc                                  >  2:1790912/1‑100 (MQ=255)
                                                                aCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCAGACACAGTTCA‑CCCCTTGTGAAGTAACTCGATAAATGTCTCTTATCCTGtt                                          >  2:2480584/1‑92 (MQ=255)
                                                                 cGTTTTTGATGATAAATGCAGTGTCGGATGTGGCGTATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCa                                 <  1:687381/100‑1 (MQ=255)
                                                                 cGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTTTGAAGTAACTCCATAGCTGCTTCATATCCTGTAAATCCGGc                                  <  2:1568578/99‑1 (MQ=255)
                                                                 cGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGc                                  <  1:3128309/99‑1 (MQ=255)
                                                                 cGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCa                                 >  1:2059576/1‑100 (MQ=255)
                                                                    ttttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaata                              <  2:576893/100‑1 (MQ=255)
                                                                     tttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGATATTGTTCA‑CCCCTTGTGAAGTAACTCCATACATGCCTCATATCCTGTTAACCCGGCaataa                             <  2:3173281/100‑1 (MQ=255)
                                                                     tttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataa                             <  1:3167790/100‑1 (MQ=255)
                                                                        gatgatAGATGCAGTGTCGGATGCGGCGCATAAGACACAGTTCA‑CCACTTGTGAAGTGACTCCCTAATTGCATCATATTCTGTTAAGCCGGCaataa                             >  2:2157849/1‑97 (MQ=255)
                                                                         atgatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtctc                         >  1:3103829/1‑100 (MQ=255)
                                                                           gatAAATGCAGTGTCGGATGCGGCGCATCCGACACCGTTCA‑CCCTTTGTGAAGTAACTCCATAAATGCCTCGTATACTGTTAACCCAGCAATAAtatc                         >  1:3720372/1‑98 (MQ=255)
                                                                           gatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccg                       >  2:1889750/1‑100 (MQ=255)
                                                                           gatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTcc                        >  2:1236922/1‑99 (MQ=255)
                                                                            atAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgc                      >  1:3721625/1‑100 (MQ=255)
                                                                               aaTGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCgg                   <  1:1565251/100‑1 (MQ=255)
                                                                               aaTGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCgg                   >  2:408007/1‑100 (MQ=255)
                                                                                 tGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTACGCCGGAt                 >  2:3784900/1‑100 (MQ=255)
                                                                                  gCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTTAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATg                <  1:117403/100‑1 (MQ=255)
                                                                                  gCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATg                >  2:3055044/1‑100 (MQ=255)
                                                                                    aGTGTCGGTTGCTGCGTATCCGACATAGTTCA‑CCTCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCt              <  1:2038419/100‑1 (MQ=255)
                                                                                    aGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGc               >  1:2618185/1‑99 (MQ=255)
                                                                                    aGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCt              <  2:470293/100‑1 (MQ=255)
                                                                                    aGTGTCGGATGCGGCGCATCCGACACAGTTAATCCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGc               >  2:730828/1‑100 (MQ=255)
                                                                                     gtgtCGGATGCGTCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTg             >  1:2572394/1‑100 (MQ=255)
                                                                                     gtgtCGGATGCGCCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAGAAGCCTCATCTCGTGTTAACCCTGCAATAATGTCGGCCGGATTCt              >  2:1505435/1‑99 (MQ=37)
                                                                                       gtCGGATGCGGCGCATCCTACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac           <  2:1501040/100‑1 (MQ=255)
                                                                                       gtCGGATGCGCCGGATCCGACACAGTTCA‑CCCCTTGTGCAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCCCCGCCGGATGCTGac           >  2:2416534/1‑100 (MQ=255)
                                                                                       gtCGGATGCGCCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac           >  1:995233/1‑100 (MQ=255)
                                                                                       gtCGGATGCGCCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac           >  1:1764556/1‑100 (MQ=255)
                                                                                       gtCGGATGCGCCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac           >  2:424594/1‑100 (MQ=255)
                                                                                       gtCGGATGCGCCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGACTCATATCATGTTAACCCAACAATAATCTCCACCGGATGCTGac           >  1:405747/1‑100 (MQ=255)
                                                                                       gtCGGATGCGCCGCATCCGACACAGTTCA‑CAACTTGTGAAGAAACTCCGTAAATGCCTCTCATCCCGTTAATCCGGCAATAATCTCCGCCGGATTCTGac           >  2:3310934/1‑100 (MQ=255)
                                                                                        tCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGaca          >  1:1389046/1‑100 (MQ=255)
                                                                                        tCGGATGCGCCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGaca          >  2:3894764/1‑100 (MQ=255)
                                                                                         cGGATGCTGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGaca          <  2:2568254/99‑1 (MQ=255)
                                                                                         cGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGacac         <  2:2618161/100‑1 (MQ=255)
                                                                                         cGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGacac         <  1:1182406/100‑1 (MQ=255)
                                                                                           gATGCGGCGCATCCGAAACAGTTCA‑CCCCTTGGGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGacac         <  2:3847356/98‑1 (MQ=255)
                                                                                             tGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCAGGCAATAATCTCCGCCGGATGCTGACACAata     <  2:1364406/100‑1 (MQ=255)
                                                                                              gcggcgCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAatat    >  1:1349395/1‑100 (MQ=255)
                                                                                              gcggcgCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAatat    <  2:1799021/100‑1 (MQ=255)
                                                                                               cggcgCCTCCGAAACAGTTGT‑CCCCTTGTGAAGTAACTCTATAAATGCCTCATATCCTGTTATTCCGGCATTAATTTCCGCCGGATGCTGACAAActat    <  2:670584/99‑1 (MQ=255)
                                                                                               cggcgCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATg   >  1:889394/1‑100 (MQ=255)
                                                                                                ggcgCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGc  >  2:1852550/1‑100 (MQ=255)
                                                                                                ggcgCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGc  >  2:779349/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
CTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCA‑CCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGC  >  NC_000913/4087994‑4088189

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑