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.

CTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGC  >  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)
                 tACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAc                                                                                 <  1:3781678/100‑1 (MQ=255)
                  aCGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAcc                                                                                <  1:1883979/100‑1 (MQ=255)
                  aCGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAcc                                                                                <  1:3277486/100‑1 (MQ=255)
                    gcgAGTACAACGAAATGGTTGATGACTATAATAGTCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAcccc                                                                              <  1:502729/100‑1 (MQ=255)
                       actACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTtg                                                                           >  1:2648040/3‑100 (MQ=255)
                          aCAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGa                                                                        <  2:3407680/100‑1 (MQ=255)
                           cAACGAAATGGTTGATGACTATAATAGGCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGaa                                                                       <  2:3260457/100‑1 (MQ=255)
                           cAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTTTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGaa                                                                       >  1:629813/1‑100 (MQ=255)
                           cAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGaa                                                                       <  1:481757/100‑1 (MQ=255)
                           cAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGAGTCGAAAGAGGCGCAACCGACACAGTTCACCCCTTGTGaa                                                                       >  1:233504/1‑100 (MQ=255)
                            aaCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAg                                                                      <  2:934962/100‑1 (MQ=255)
                             aCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGt                                                                     <  1:2097703/100‑1 (MQ=255)
                             aCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGt                                                                     <  1:3850949/100‑1 (MQ=255)
                             aCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCACATCCGACACAGTTCACCCCTTGTGAAGt                                                                     >  1:34374/1‑100 (MQ=255)
                              cGAAATGGTTTATAACTATAATAGCCTACGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCAGCCGACACAGTTCACCCCTTGTGAAGTa                                                                    <  1:530106/100‑1 (MQ=255)
                               gAAATGATTGATGTATATAATAGCCTGCGTTAAACGTTTTAGATACTATATCCTGTGTCGGATGCGGCGCATCCAACACAGTTATCCACTTTTGa                                                                        <  2:1998574/95‑1 (MQ=255)
                                    ggTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCa                                                              <  1:2028357/100‑1 (MQ=255)
                                    ggTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCa                                                              >  2:357089/1‑100 (MQ=255)
                                      ttaatgaCTATAATATCCTTCTTTAATCTTTTTTTTTTATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATa                                                            <  1:292874/97‑1 (MQ=255)
                                        gatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGTAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATaaa                                                          <  1:1105704/100‑1 (MQ=255)
                                        gatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATaaa                                                          >  1:35614/1‑100 (MQ=255)
                                         atgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAAt                                                         <  1:1707204/100‑1 (MQ=255)
                                          tgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATg                                                        <  1:357882/100‑1 (MQ=255)
                                           gaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGAATGCGACGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGc                                                       >  1:1026136/1‑100 (MQ=255)
                                            aCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGcc                                                      <  1:1068546/100‑1 (MQ=255)
                                            aCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGcc                                                      <  2:2804388/100‑1 (MQ=255)
                                             cTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCt                                                     >  2:1416904/1‑100 (MQ=255)
                                                  ataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCatat                                                <  2:1077322/100‑1 (MQ=255)
                                                   taGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATc                                               >  2:2866301/1‑100 (MQ=255)
                                                    aGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATcc                                              >  2:264746/1‑100 (MQ=255)
                                                     gCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCt                                             >  1:3037265/1‑100 (MQ=255)
                                                       cTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGt                                           >  2:3683950/1‑100 (MQ=255)
                                                        tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGtt                                          <  1:3788643/100‑1 (MQ=255)
                                                         gCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTa                                         >  1:3603016/1‑100 (MQ=255)
                                                              aaaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCg                                    <  2:539715/100‑1 (MQ=255)
                                                               aaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCgg                                   <  1:3104551/100‑1 (MQ=255)
                                                                aCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGc                                  >  2:1790912/1‑100 (MQ=255)
                                                                aCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCAGACACAGTTCACCCCTTGTGAAGTAACTCGATAAATGTCTCTTATCCTGtt                                          >  2:2480584/1‑92 (MQ=255)
                                                                 cGTTTTTGATGATAAATGCAGTGTCGGATGTGGCGTATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCa                                 <  1:687381/100‑1 (MQ=255)
                                                                 cGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTTTGAAGTAACTCCATAGCTGCTTCATATCCTGTAAATCCGGc                                  <  2:1568578/99‑1 (MQ=255)
                                                                 cGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGc                                  <  1:3128309/99‑1 (MQ=255)
                                                                 cGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCa                                 >  1:2059576/1‑100 (MQ=255)
                                                                    ttttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaata                              <  2:576893/100‑1 (MQ=255)
                                                                     tttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGATATTGTTCACCCCTTGTGAAGTAACTCCATACATGCCTCATATCCTGTTAACCCGGCaataa                             <  2:3173281/100‑1 (MQ=255)
                                                                     tttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataa                             <  1:3167790/100‑1 (MQ=255)
                                                                        gatgatAGATGCAGTGTCGGATGCGGCGCATAAGACACAGTTCACCACTTGTGAAGTGACTCCCTAATTGCATCATATTCTGTTAAGCCGGCaataa                             >  2:2157849/1‑97 (MQ=255)
                                                                         atgatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtctc                         >  1:3103829/1‑100 (MQ=255)
                                                                           gatAAATGCAGTGTCGGATGCGGCGCATCCGACACCGTTCACCCTTTGTGAAGTAACTCCATAAATGCCTCGTATACTGTTAACCCAGCAATAAtatc                         >  1:3720372/1‑98 (MQ=255)
                                                                           gatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccg                       >  2:1889750/1‑100 (MQ=255)
                                                                           gatAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTcc                        >  2:1236922/1‑99 (MQ=255)
                                                                            atAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgc                      >  1:3721625/1‑100 (MQ=255)
                                                                               aaTGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCgg                   <  1:1565251/100‑1 (MQ=255)
                                                                               aaTGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCgg                   >  2:408007/1‑100 (MQ=255)
                                                                                 tGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTACGCCGGAt                 >  2:3784900/1‑100 (MQ=255)
                                                                                  gCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTTAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATg                <  1:117403/100‑1 (MQ=255)
                                                                                  gCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATg                >  2:3055044/1‑100 (MQ=255)
                                                                                    aGTGTCGGTTGCTGCGTATCCGACATAGTTCACCTCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCt              <  1:2038419/100‑1 (MQ=255)
                                                                                    aGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGc               >  1:2618185/1‑99 (MQ=255)
                                                                                    aGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCt              <  2:470293/100‑1 (MQ=255)
                                                                                    aGTGTCGGATGCGGCGCATCCGACACAGTTAATCCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGc               >  2:730828/1‑100 (MQ=255)
                                                                                     gtgtCGGATGCGTCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTg             >  1:2572394/1‑100 (MQ=255)
                                                                                     gtgtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAGAAGCCTCATCTCGTGTTAACCCTGCAATAATGTCGGCCGGATTCt              >  2:1505435/1‑99 (MQ=37)
                                                                                       gtCGGATGCGGCGCATCCTACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac           <  2:1501040/100‑1 (MQ=255)
                                                                                       gtCGGATGCGCCGGATCCGACACAGTTCACCCCTTGTGCAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCCCCGCCGGATGCTGac           >  2:2416534/1‑100 (MQ=255)
                                                                                       gtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac           >  1:995233/1‑100 (MQ=255)
                                                                                       gtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac           >  1:1764556/1‑100 (MQ=255)
                                                                                       gtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGac           >  2:424594/1‑100 (MQ=255)
                                                                                       gtCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGACTCATATCATGTTAACCCAACAATAATCTCCACCGGATGCTGac           >  1:405747/1‑100 (MQ=255)
                                                                                       gtCGGATGCGCCGCATCCGACACAGTTCACAACTTGTGAAGAAACTCCGTAAATGCCTCTCATCCCGTTAATCCGGCAATAATCTCCGCCGGATTCTGac           >  2:3310934/1‑100 (MQ=255)
                                                                                        tCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGaca          >  1:1389046/1‑100 (MQ=255)
                                                                                        tCGGATGCGCCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGaca          >  2:3894764/1‑100 (MQ=255)
                                                                                         cGGATGCTGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGaca          <  2:2568254/99‑1 (MQ=255)
                                                                                         cGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGacac         <  2:2618161/100‑1 (MQ=255)
                                                                                         cGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGacac         <  1:1182406/100‑1 (MQ=255)
                                                                                           gATGCGGCGCATCCGAAACAGTTCACCCCTTGGGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGacac         <  2:3847356/98‑1 (MQ=255)
                                                                                             tGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCAGGCAATAATCTCCGCCGGATGCTGACACAata     <  2:1364406/100‑1 (MQ=255)
                                                                                              gcggcgCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAatat    >  1:1349395/1‑100 (MQ=255)
                                                                                              gcggcgCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAatat    <  2:1799021/100‑1 (MQ=255)
                                                                                               cggcgCCTCCGAAACAGTTGTCCCCTTGTGAAGTAACTCTATAAATGCCTCATATCCTGTTATTCCGGCATTAATTTCCGCCGGATGCTGACAAActat    <  2:670584/99‑1 (MQ=255)
                                                                                               cggcgCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATg   >  1:889394/1‑100 (MQ=255)
                                                                                                ggcgCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGc  >  2:1852550/1‑100 (MQ=255)
                                                                                                ggcgCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGc  >  2:779349/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
CTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGC  >  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: