Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,744,794 A→G 9.6% G156G (GGA→GGG)  yiaM → 2,3‑diketo‑L‑gulonate TRAP transporter small permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,744,7940AG 9.6% 218.8 / 3.2 97G156G (GGA→GGG) yiaM2,3‑diketo‑L‑gulonate TRAP transporter small permease protein
Reads supporting (aligned to +/- strand):  ref base A (44/40);  new base G (3/8);  total (49/48)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAA  >  NC_000913/3744646‑3744938
                                                                                                                                                    |                                                                                                                                                
ggAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGa                                                                                                                                                  >  2:675690/1‑149 (MQ=255)
  aGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCCGGTGTTGTTTCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCACGGAGc                                                                                                                                                >  1:21642/1‑149 (MQ=255)
   gCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCt                                                                                                                                               >  1:563384/1‑149 (MQ=255)
     gATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCCCTTGTCTGCCCACCCGCCTCGTTATCGCATTCTTTGAACTGCGTCCTTTATGTCAACTCACCCCGCGTAGCCATTCGCTACCGTCGCCACCGCCCGGAGCTTa                                                                                                                                             >  2:104840/1‑149 (MQ=255)
       ttattCACCGTGCCTCGGTTTACTCATCGGCCGGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGCCCATTATATCAACTCACCACGCGTAGCAATTCTCTACCGTCGCCACCGCAAGGAGCTTAAt                                                                                                                                           <  1:213719/149‑1 (MQ=255)
       ttattCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGGCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCAACCGCAAGGAGCTTAAt                                                                                                                                           >  1:513842/1‑149 (MQ=255)
         attCACCGATCCTCGGTTTCCCCATCGGCCCGATGTATGTCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCCACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCa                                                                                                                                         <  2:172433/149‑1 (MQ=255)
         attCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTCTATTAACTCTTCATGATGCGCAAACGTCTAACATCCCCCCCGCAAGGAGCTTAATCa                                                                                                                                         >  2:640576/1‑149 (MQ=255)
         attCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCGCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTCACGTCGCCACCGCAAGGACCTTAATCa                                                                                                                                         >  2:190289/1‑149 (MQ=255)
           tCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATg                                                                                                                                       >  1:402662/1‑149 (MQ=255)
           tCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATg                                                                                                                                       >  1:672108/1‑149 (MQ=255)
            cACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATgg                                                                                                                                      >  2:707123/1‑149 (MQ=255)
             aCCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCAGCTGTTGCGTTTTGGCAGTCTGTGGACTGCGTCAGTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATAATGGc                                                                                                                                     >  2:648059/1‑149 (MQ=255)
             aCCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGc                                                                                                                                     >  1:685015/1‑149 (MQ=255)
             aCCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTAGCCACCGCAAGGAGCTTAATCATGGc                                                                                                                                     >  1:92666/1‑149 (MQ=255)
             aCCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATAGCAACTCATCGCGCGTGGCAATTCCCTAACGTCGCCACCGCAAGGGGCTTAATCATGGc                                                                                                                                     >  2:293271/1‑149 (MQ=255)
               cGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCACCTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCtg                                                                                                                                   <  1:635851/149‑1 (MQ=255)
               cGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCCCAAGGAGCTTAATCATGGCtg                                                                                                                                   >  1:477424/1‑149 (MQ=255)
                  tCCTCGTTTTTCCGATCGAACTGATGTATGCCTTTTGTATTCCCACCATCCTCGTTATCGAATTCTTTGAACTGCGTCATTTATTTCAACTCATCACGAGTACCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGc                                                                                                                                <  2:238588/149‑1 (MQ=255)
                        gTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCAGGGCTGTGCTGAttt                                                                                                                          >  1:213949/1‑149 (MQ=255)
                            aCCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCt                                                                                                                      <  1:694803/149‑1 (MQ=255)
                             cccTTCGCCGTGCTGGATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGTCTTCTTTGAACTGCGTCATTTATATCACCTCCTCACGCGTAGTAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCCTGGCTGTGCTGATTTTTCTg                                                                                                                     <  1:111433/149‑1 (MQ=255)
                             cccATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTg                                                                                                                     <  1:106759/149‑1 (MQ=255)
                             cccATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATACCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTg                                                                                                                     <  1:76967/149‑1 (MQ=255)
                              ccATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCCTCCATCGGAGCAATTCTCTAACGTCGCCACAGCAAGGCGCTTAATCATGGCTGTGCTGATTGTTCTgg                                                                                                                    >  1:360752/1‑149 (MQ=255)
                              ccATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTgg                                                                                                                    >  1:130386/1‑149 (MQ=255)
                              ccATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGAGTTTTCCgg                                                                                                                    >  2:214542/1‑149 (MQ=255)
                               cATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCCCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTggg                                                                                                                   <  1:105048/149‑1 (MQ=255)
                                 tCGGCCTGATGTATGCCGCTTGTCTGCCCGCCAGCCTCGTTGTCGCATTCTTTGGACTGCGTCATGTTTATTATTTCATCACGCGTTGCGATTCTTTTACGTGACGACACACAGGGGCTTAATCATGGCTGTGGCGAgtttttggggag                                                                                                                 >  2:189154/1‑142 (MQ=255)
                                   ggCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGctgt                                                                                                               >  1:573842/1‑149 (MQ=255)
                                   ggCCTGAGGTATGCCGCTTGTCTGCCCACCCGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGctgt                                                                                                               <  2:661246/149‑1 (MQ=255)
                                    gCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCGCGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGctgtc                                                                                                              <  2:627769/149‑1 (MQ=255)
                                     ccTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGctgtct                                                                                                             >  1:508188/1‑149 (MQ=255)
                                       tGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTGTCTGGGcggtctgt                                                                                                           >  2:615116/1‑149 (MQ=255)
                                       tGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATGTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGctgtccgt                                                                                                           >  2:565822/1‑149 (MQ=255)
                                          tGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTTTCTGTTgg                                                                                                        <  2:203077/149‑1 (MQ=255)
                                            tATGCCGCTAGTCTGCACACCAGCCTCGATATCGCAGTCTTTGGACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCGGTGCTGATTTTTCGGGGCAGTCCGTTGggt                                                                                                      >  2:68752/1‑149 (MQ=255)
                                              tGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTCGggtgg                                                                                                    >  1:14686/1‑149 (MQ=255)
                                              tGCCGCTTGTCTGCCCAACAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCGACTCATCACGCGTAGCAATTCTCTAACGTCGCCACAGCAAGGAGCTTAATCATGGCTGTCCGGATTTTTCTGGGCTGTCTGTTGggtgg                                                                                                    >  2:555258/1‑149 (MQ=255)
                                               gCCGCTTGTCTGCCGACCAGCCTCGAAATGTACAACGCGGCACTGCGTCCTTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGCGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGgggggt                                                                                                   >  1:555279/1‑149 (MQ=255)
                                                ccGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCATCTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTa                                                                                                  <  2:696648/149‑1 (MQ=255)
                                                  gCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCGGTGCTGATTTTTCTGGGCTGGCTGTTGGGGGGTAtt                                                                                                >  2:320293/1‑149 (MQ=255)
                                                   tctggcTGCCCCCCAGCCTCGTTATCTCATTCTTCGAACTGCGTCATTTATATCATCTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTg                                                                                               <  2:736442/144‑1 (MQ=255)
                                                   cctGTCGGCCCACCAGCCTCCTTATCCCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCCAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTg                                                                                               <  1:121253/147‑1 (MQ=255)
                                                      gTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGTGTCATTTGTATCAACTCATCACGCGTAGCAATTCGCTAACGTCGCCACCGCAAGGAGCTTAATCATGACTGTACTGATTTTTCTGGGTGGGCTGTTGGGTGGTAGTGCTa                                                                                            >  2:563176/1‑149 (MQ=255)
                                                       tCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTAt                                                                                           <  1:670890/149‑1 (MQ=255)
                                                        cTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGGTTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATc                                                                                          >  2:466374/1‑149 (MQ=255)
                                                           ccccccAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTTTGTTGGGTGGTATTGCTATCGGt                                                                                       <  1:216831/149‑1 (MQ=255)
                                                            ccaccaGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGtt                                                                                      >  2:384718/1‑149 (MQ=255)
                                                             ccccaGCCCCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCCTCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGttt                                                                                     <  1:194931/147‑1 (MQ=255)
                                                              accaGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCATCTCTTCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTg                                                                                    <  1:581955/149‑1 (MQ=255)
                                                                   ccTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTAt                                                                               >  1:327663/1‑149 (MQ=255)
                                                                   ccTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTAt                                                                               >  1:133923/1‑149 (MQ=255)
                                                                           tCGCATTCTTTGAACTGCGTCATTTATATCAGCTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAAGCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTggg                                                                       <  2:287098/149‑1 (MQ=255)
                                                                           tCGCATTCTTTGAACTGCGTCATTTATATCAACCCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTggg                                                                       <  1:194914/149‑1 (MQ=255)
                                                                             gCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGGGCTTAATCATGGCTGTGCTGATTTTACTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCa                                                                     >  1:541725/1‑149 (MQ=255)
                                                                               aTTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCGCt                                                                   >  1:469628/1‑149 (MQ=255)
                                                                                   tttGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCCCtgttg                                                               >  2:87643/1‑149 (MQ=255)
                                                                                   tttGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCGGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGGAcgcttt                                                               >  1:626516/1‑144 (MQ=255)
                                                                                   tttGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTATAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCGGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCGATCGGTTTTGCTATTGCCTGGGCAcgggtg                                                               >  2:174656/1‑144 (MQ=255)
                                                                                    ttGAACTTGGTCATTTATATCAACTCATCACGCGGAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTAGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACtgttgt                                                              <  1:308137/149‑1 (MQ=255)
                                                                                     tGACCTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTATCGTCGCCACCGCTAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACtgttgtt                                                             <  2:92216/149‑1 (MQ=255)
                                                                                      gAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACtgttgttg                                                            <  2:110011/149‑1 (MQ=255)
                                                                                        aCTGCGTCATTTATATCACCTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTtgtg                                                          <  1:394723/149‑1 (MQ=255)
                                                                                         cTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCGCTACCGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGc                                                         <  1:93358/149‑1 (MQ=255)
                                                                                            cGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCgga                                                      >  1:418293/1‑148 (MQ=255)
                                                                                             gTCATTTATATCACCTCATCACGCGTAGCAATTCGCGACCGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCgggg                                                     <  1:214542/149‑1 (MQ=255)
                                                                                                 tttATATCGACTCAGCACGCGTAGCAATTCTCTAACGTCGCACCCGCAAGGGGCTTACTCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGggcggc                                                 <  2:267107/149‑1 (MQ=255)
                                                                                                 tttATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGggcggc                                                 >  1:9827/1‑149 (MQ=255)
                                                                                                  ttATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGgcggcg                                                >  2:218176/1‑149 (MQ=255)
                                                                                                    atatCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGGTGTTGTGCGGGGCGGCGtt                                              >  2:700488/1‑149 (MQ=255)
                                                                                                       tCAACTCCTCCCGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAAt                                           <  1:157874/149‑1 (MQ=255)
                                                                                                       tCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAAt                                           <  1:617299/149‑1 (MQ=255)
                                                                                                          aCTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGGGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGtt                                        <  2:475185/149‑1 (MQ=255)
                                                                                                          aCTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGtt                                        >  2:137476/1‑149 (MQ=255)
                                                                                                            tcatcaCGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTct                                      <  1:148775/149‑1 (MQ=255)
                                                                                                               tcaCGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTctggc                                   <  2:327663/149‑1 (MQ=255)
                                                                                                                caCGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTCATTTTTCTCGTCTGTCTGTTGCGTGGTATTGCTACCGGTTCGCCTATTTCCTGGGCACTGTTGCTGTTCGGGGCGGCGTTAATGTTctggct                                  <  1:563176/149‑1 (MQ=255)
                                                                                                                 gcgcgTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTctggctg                                 <  2:56300/148‑1 (MQ=255)
                                                                                                                 aCGCGTAGCATTTCTCGAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTctggctg                                 <  2:170033/149‑1 (MQ=255)
                                                                                                                 aCGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTctggctg                                 <  1:670185/149‑1 (MQ=255)
                                                                                                                     gTAGCAATTCTCTACCGTCGCCACCGCAAGGGGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACa                             <  2:403165/149‑1 (MQ=255)
                                                                                                                     gTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACa                             >  1:330826/1‑149 (MQ=255)
                                                                                                                            ttCTCTAGCGTCGCCACCGCAAGAAGATTAATCATGGCTGTGCTGATTTTTCCGGGCTGTCTGTTGGGGGGTATTGCTATCGGTTTGCCTATCGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGa                      <  1:244909/149‑1 (MQ=255)
                                                                                                                            ttCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGa                      >  1:300396/1‑149 (MQ=255)
                                                                                                                                ctAACGTCGCCACCGCCAAGAGCTTAATCATGGCCGTGCTGATTTTTCGGGGCTGGCTGGTGTGTGGTGTTGCCATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTTCGGGGCGGCGTTAATGGTCTGGCTGGACCTGTTTGAGGTc                  >  2:180726/1‑149 (MQ=255)
                                                                                                                                 tAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGGTAATGTTCTGGCTGGACATGTTTGATGTcc                 >  1:228696/1‑149 (MQ=255)
                                                                                                                                 caaCCTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTTCCTATTGCCTGGGCACTGTTGCTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTcc                 <  1:176967/148‑1 (MQ=255)
                                                                                                                                  aaCGTCGCCACCGCAAGGGGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCa                <  2:541725/149‑1 (MQ=255)
                                                                                                                                        gctcccGCAGGGGGGTTAATCTTGGCTGTGCTGATTTTTTTGGGCTGTCCGCTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCAt          <  2:380156/145‑1 (MQ=255)
                                                                                                                                        gCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCAt          >  1:651698/1‑149 (MQ=255)
                                                                                                                                         ccACCGCAAGGGGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATg         <  2:573842/149‑1 (MQ=255)
                                                                                                                                          cagcGCAAGGGGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATgg        <  1:320293/146‑1 (MQ=255)
                                                                                                                                          cACCGCAAGGGTCTTGACCCTGGGTGGGTTGATTTTTCTGGGCGGTCTGTTAGGCGGTATTGGTATCGGTATGCCTACTGCCTGGGCACTGTTGTTGTGGGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATgg        <  2:3475/149‑1 (MQ=255)
                                                                                                                                          cACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATgg        <  1:308500/149‑1 (MQ=255)
                                                                                                                                              gCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGcaca    >  2:472275/1‑149 (MQ=255)
                                                                                                                                                aaGGAGCTTAATCATGGCTGTGCTGTTTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACaaa  <  1:466374/149‑1 (MQ=255)
                                                                                                                                                aaGGAGCTTAATCATGGCTGGGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACaaa  <  1:693670/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGTTGTTGTGCGGGGCGGCGTTAATGTTCTGGCTGGACATGTTTGATGTCCAGATCATGGCACAAA  >  NC_000913/3744646‑3744938

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