Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,005,621 T→G 8.2% L184R (CTG→CGG) ‡ PP_0867 → FecA‑like outer membrane receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,005,6210TG 8.2% 300.9 / 14.2 98L184R (CTG→CGG) ‡PP_0867FecA‑like outer membrane receptor
Reads supporting (aligned to +/- strand):  ref base T (44/46);  new base G (4/4);  total (48/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

GAGCGTGGTACAGAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATTTCTTCTGGCAACCTC  >  NC_002947/1005475‑1005760
                                                                                                                                                  |                                                                                                                                           
gagCGTGGTACAGAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGt                                                                                                                                           >  4:283049/1‑149 (MQ=255)
 agCGTGGTACAGAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTc                                                                                                                                          <  4:203737/149‑1 (MQ=255)
 agCGTGGTACAGAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTc                                                                                                                                          >  2:167063/1‑149 (MQ=255)
  gCGTGGTACAGAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTc                                                                                                                                          <  4:138592/148‑1 (MQ=255)
  gCGTGGTACAGAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCg                                                                                                                                         >  2:228284/1‑149 (MQ=255)
  gCGTGGTACAGAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCg                                                                                                                                         <  5:263135/149‑1 (MQ=255)
       gTACAGAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGGCTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGaa                                                                                                                                    >  8:41429/1‑149 (MQ=255)
       gTACAGAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGaa                                                                                                                                    <  2:159208/149‑1 (MQ=255)
        tACAGAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACAGGTCGCTgcca                                                                                                                                   <  1:410125/149‑4 (MQ=255)
         aCAGAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACg                                                                                                                                  >  4:120747/1‑149 (MQ=255)
            gAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGTGTTCAAGTGCAGGACTCCAACGGGTCCGGTGGCAGCGACCTGTCGCTGAACGTgg                                                                                                                               <  2:26830/149‑1 (MQ=255)
            gAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTgg                                                                                                                               <  6:444103/149‑1 (MQ=255)
             aaCCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTggg                                                                                                                              >  4:208684/1‑149 (MQ=255)
               ccaccCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCg                                                                                                                            >  3:371082/1‑149 (MQ=255)
                  cccTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTAc                                                                                                                         <  4:371082/149‑1 (MQ=255)
                   ccTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTAcg                                                                                                                        >  8:143032/1‑149 (MQ=255)
                       gggCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCgg                                                                                                                     <  8:385722/148‑1 (MQ=255)
                       gggCAAGCACTGTGGTTCGCCGCGAGGCAATGGTTGAGAAGGGCGCAATGAACGTGCTCGACGTGCTGCGCGGCATGGCCGGCGTGCAAGTGCAGTACTCCAACGGGACCGTTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGc                                                                                                                    <  2:241430/149‑1 (MQ=255)
                        ggCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGcc                                                                                                                   <  2:116915/149‑1 (MQ=255)
                            cgcacTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTAc                                                                                                               <  6:411771/149‑1 (MQ=255)
                            cgcacTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTAc                                                                                                               <  3:170575/149‑1 (MQ=255)
                             gcacTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACt                                                                                                              >  5:107626/1‑149 (MQ=255)
                              cacTGT‑GTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTc                                                                                                            <  1:116079/149‑1 (MQ=255)
                                   tgGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTcgcg                                                                                                         <  8:16298/148‑1 (MQ=255)
                                   gggTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCgcgc                                                                                                        <  1:46367/148‑1 (MQ=255)
                                      tACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTg                                                                                                     >  5:141295/1‑149 (MQ=255)
                                        cgccgcGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTc                                                                                                   >  4:251415/1‑149 (MQ=255)
                                            gcgAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTcgccg                                                                                               >  1:67398/1‑149 (MQ=255)
                                            gcgAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCACTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTcgccg                                                                                               <  8:95216/149‑1 (MQ=255)
                                             cgAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTcgccgc                                                                                              >  1:234531/1‑149 (MQ=255)
                                                   aaTGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCa                                                                                        <  7:256870/149‑1 (MQ=255)
                                                     tggtggAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACg                                                                                      >  5:76405/1‑149 (MQ=255)
                                                     tggtggAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACg                                                                                      <  1:228284/149‑1 (MQ=255)
                                                      ggtggAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACgg                                                                                     >  6:203780/1‑149 (MQ=255)
                                                      ggtggAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACgg                                                                                     >  6:210054/1‑149 (MQ=255)
                                                       gtggAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACgg                                                                                     >  7:265782/1‑148 (MQ=255)
                                                         ggAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTAc                                                                                  >  2:369611/1‑149 (MQ=255)
                                                          gagaAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTc                                                                                                   >  4:288443/1‑131 (MQ=255)
                                                          gagaAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTc                                                                                                   <  3:288443/131‑1 (MQ=255)
                                                                  cgcAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGGCTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATg                                                                         <  7:41429/149‑1 (MQ=255)
                                                                   gcAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATgg                                                                        >  1:276196/1‑149 (MQ=255)
                                                                     aaTGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCa                                                                      >  2:415827/1‑149 (MQ=255)
                                                                     aaTGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCa                                                                      >  3:44066/1‑149 (MQ=255)
                                                                       tGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATc                                                                    <  7:143032/149‑1 (MQ=255)
                                                                       tGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATc                                                                    <  2:234531/149‑1 (MQ=255)
                                                                         aaCGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATccc                                                                  >  3:249229/1‑149 (MQ=255)
                                                                         aaCGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATcc                                                                   <  1:167063/148‑1 (MQ=255)
                                                                          aCGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCg                                                                 >  8:129359/1‑149 (MQ=255)
                                                                           cGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCgg                                                                <  6:76405/149‑1 (MQ=255)
                                                                              gcgcGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGGGCTCCACGGTACTGATCGATGGCATccc                                                                  <  1:242257/144‑1 (MQ=255)
                                                                              gcgcGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATccc                                                                  >  2:242257/1‑144 (MQ=255)
                                                                                gcgcCGTGCTGCGCTGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGTACGTGGGCGTTCGCGGCCTTACTTCGCGCCTGTCGCCGCGCGCCACGGTACTGATCGATGGCATCCCGgccgcc                                                           <  6:328349/149‑1 (MQ=255)
                                                                                     gTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTGCTGATCGATGGCATCCCGGCCGCCTTcgc                                                      <  6:141295/149‑1 (MQ=255)
                                                                                       gctgcGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCgcgc                                                    >  7:346502/1‑149 (MQ=255)
                                                                                        ctgcGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGcc                                                   <  2:67398/149‑1 (MQ=255)
                                                                                         tgcGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCg                                                  >  3:103619/1‑149 (MQ=255)
                                                                                           cgcgGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGt                                                 <  5:203780/148‑1 (MQ=255)
                                                                                           cgcgGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGt                                                 <  5:210054/148‑1 (MQ=255)
                                                                                           cgcgGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGt                                                 <  3:208684/148‑1 (MQ=255)
                                                                                           cgcgGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTa                                                >  5:345121/1‑149 (MQ=255)
                                                                                           cgcgGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTa                                                <  6:107626/149‑1 (MQ=255)
                                                                                             cgGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATg                                              <  1:145574/149‑1 (MQ=255)
                                                                                             cgGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATg                                              >  8:388859/1‑149 (MQ=255)
                                                                                             cgGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATg                                              <  3:120747/149‑1 (MQ=255)
                                                                                               gCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGgc                                            <  3:251415/149‑1 (MQ=255)
                                                                                                  tCCCCGGCGTTCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGgccag                                         <  1:384020/149‑1 (MQ=255)
                                                                                                  tCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGgc                                            >  7:248114/1‑146 (MQ=255)
                                                                                                  tCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGgc                                            <  8:248114/146‑1 (MQ=255)
                                                                                                    cccGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGgccagc                                        >  7:179545/1‑148 (MQ=255)
                                                                                                             cAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCt                                 >  1:208717/1‑146 (MQ=255)
                                                                                                             cAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCt                                 <  2:208717/146‑1 (MQ=255)
                                                                                                             cAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTc                              >  3:28513/1‑149 (MQ=255)
                                                                                                                gTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGAt                           <  2:75436/149‑1 (MQ=255)
                                                                                                                 tGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATg                          <  1:415827/149‑1 (MQ=255)
                                                                                                                 tGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATg                          <  4:249229/149‑1 (MQ=255)
                                                                                                                 tGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATg                          <  8:346502/149‑1 (MQ=255)
                                                                                                                   cAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCCCAGCTGTCGATgg                         >  6:409129/1‑148 (MQ=255)
                                                                                                                      gACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCTATgg                         <  7:337772/145‑1 (MQ=255)
                                                                                                                      gACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATgg                         >  8:337772/1‑145 (MQ=255)
                                                                                                                         tCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGAt                  >  6:33618/1‑149 (MQ=255)
                                                                                                                              cGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATttct              >  4:78329/1‑148 (MQ=255)
                                                                                                                               gggACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATttcttc            >  3:367232/1‑149 (MQ=255)
                                                                                                                                 gACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCGACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCCGCCACAGCTGTCGATGGCACCGATTTCTTCTg          >  4:178599/1‑149 (MQ=255)
                                                                                                                                 gACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATTTCTTCTg          <  4:28513/149‑1 (MQ=255)
                                                                                                                                  aCCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATTTCTTCTgg         >  1:294697/1‑149 (MQ=255)
                                                                                                                                   ccacgttcAGCGACAGGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATTTCTTCTg          <  6:158944/140‑1 (MQ=255)
                                                                                                                                   ccacgttcAGCGACAGGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATTTCTTCTg          >  5:158944/8‑147 (MQ=255)
                                                                                                                                    cGGTGGCAGCGACCTGTCTATGAACGTGGGCGTACGAGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGgcca                                          <  6:339101/114‑1 (MQ=255)
                                                                                                                                    cGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTg                                <  3:60146/124‑1 (MQ=255)
                                                                                                                                    cGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTg                                >  4:60146/1‑124 (MQ=255)
                                                                                                                                    cGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATAGATGGCATCCCGGCCGCCTTCGCGCCGTATGgcca                                          >  5:339101/1‑114 (MQ=255)
                                                                                                                                     acgttcAGCGACAGGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATTTCTTCTg          <  4:122980/140‑1 (MQ=255)
                                                                                                                                     acgttcAGCGACAGGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATTTCTTCTg          >  3:122980/6‑145 (MQ=255)
                                                                                                                                     acgttcAGCGACAGGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATTTCTTCTGGCaa      >  1:169732/6‑149 (MQ=255)
                                                                                                                                       tGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATTTCTTCTGGCAAcc    >  5:198767/1‑149 (MQ=255)
                                                                                                                                       gttcAGCGACAGGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGAt                  <  2:102854/132‑1 (MQ=255)
                                                                                                                                       gttcAGCGACAGGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGAt                  >  1:102854/4‑135 (MQ=255)
                                                                                                                                         gCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATTTCTTCTGGCAACCTc  >  4:296975/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
GAGCGTGGTACAGAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATTTCTTCTGGCAACCTC  >  NC_002947/1005475‑1005760

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