Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,005,623 C→A 7.2% L184M (CTG→ATG) ‡ 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,6230CA 7.2% 299.1 / 6.6 96L184M (CTG→ATG) ‡PP_0867FecA‑like outer membrane receptor
Reads supporting (aligned to +/- strand):  ref base C (50/39);  new base A (2/5);  total (52/44)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

AGCGTGGTACAGAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATTTCTTCTGGCAACCTCGACAGCATC  >  NC_002947/1005479‑1005772
                                                                                                                                                 |                                                                                                                                                      
agCGTGGTACAGAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGACAGACCCTGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCg                                                                                                                                                   <  1:1357975/150‑1 (MQ=255)
  cGTGGTACAGAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCCACGGGACCGGTGGCAGCGACCTGTCGCt                                                                                                                                                 <  2:1251028/150‑1 (MQ=255)
    tGGTACAGAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGa                                                                                                                                               >  1:2348807/1‑150 (MQ=255)
    tGGTACAGAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGa                                                                                                                                               >  1:2432649/1‑150 (MQ=255)
    tGGTACAGAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGa                                                                                                                                               >  1:637265/1‑150 (MQ=255)
    tGGTACAGAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACAGGTCGCTgc                                                                                                                                               <  2:566522/150‑2 (MQ=255)
    gggTACAGAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCCGCGACAGGTCGCTgc                                                                                                                                               <  2:1673792/149‑2 (MQ=255)
     ggTACAGAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGaa                                                                                                                                              >  1:1507536/1‑150 (MQ=255)
       tACAGAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCGACGGGACCGGTGGCAGCGACCTGTCGCTGAACg                                                                                                                                            >  2:562088/1‑150 (MQ=255)
       tACAGAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGACAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACg                                                                                                                                            >  2:2021169/1‑150 (MQ=255)
         cAGAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTg                                                                                                                                          >  2:2112723/1‑150 (MQ=255)
          aGAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGGCGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTgg                                                                                                                                        >  1:1747664/1‑150 (MQ=255)
           gAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTggg                                                                                                                                       >  1:940194/1‑150 (MQ=255)
            aaCCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACAGGTCGCTgc                                                                                                                                               <  2:1856962/142‑2 (MQ=255)
            aaCCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACAGGTCGCTgc                                                                                                                                               >  1:1856962/1‑141 (MQ=255)
            aaCCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACAGGTCGCTgc                                                                                                                                               >  1:862559/1‑141 (MQ=255)
            aaCCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACAGGTCGCTgc                                                                                                                                               <  2:862559/142‑2 (MQ=255)
              ccaccCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGt                                                                                                                                    >  1:354379/1‑150 (MQ=255)
               caccCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTa                                                                                                                                   >  1:813510/1‑150 (MQ=255)
                 cccTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTAcg                                                                                                                                 >  2:449348/1‑150 (MQ=255)
                  ccTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTAcgc                                                                                                                                >  1:388450/1‑150 (MQ=255)
                        gcacgcacTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGAACAGTGCTGCGCGGATCGGCCGGGGTGCAGGTGCAGGACTCCAACGGGACCGGTGGAAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCtt                                                                                                                          >  1:1485358/1‑150 (MQ=255)
                         cacgcacTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTa                                                                                                                         >  2:820688/1‑150 (MQ=255)
                         cacgcacTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAATGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTa                                                                                                                         >  2:1211925/1‑150 (MQ=255)
                               cTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTcgc                                                                                                                   <  1:2040479/150‑1 (MQ=255)
                                tgtgGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTcgcg                                                                                                                  <  1:939069/150‑1 (MQ=255)
                                tgtgGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGCGGCCACCAGCGATCGCAGAACGCTTCCGAACCCGGCCTTACTTcgcg                                                                                                                  >  2:1842414/1‑150 (MQ=255)
                                        gccgcGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTcgc                                                                                                          >  2:267187/1‑150 (MQ=255)
                                         ccgcGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTcgcc                                                                                                         >  1:2074724/1‑150 (MQ=255)
                                          cgcgAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCt                                                                                                                           <  2:2307975/131‑1 (MQ=255)
                                          cgcgAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCt                                                                                                                           >  1:2307975/1‑131 (MQ=255)
                                                gCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAAc                                                                                                                                             <  2:1398359/108‑1 (MQ=255)
                                                gCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAAc                                                                                                                                             >  1:1398359/1‑108 (MQ=255)
                                                   aTGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACg                                                                                               <  2:2398039/150‑1 (MQ=255)
                                                   aTGGCGGAGAGGGGCGAAATGACCGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACCTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACg                                                                                               <  1:1027524/150‑1 (MQ=255)
                                                         gagaAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACTGCTGACGACGGGCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTg                                                                                         >  1:146858/1‑150 (MQ=255)
                                                              gggCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGa                                                                                    >  2:1123603/1‑150 (MQ=255)
                                                                gcgcAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATg                                                                                  <  2:637265/150‑1 (MQ=255)
                                                                 cgcAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATgg                                                                                 <  1:282088/150‑1 (MQ=255)
                                                                 cgcAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATgg                                                                                 <  2:940194/150‑1 (MQ=255)
                                                                    aaTGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCAt                                                                              <  2:56717/150‑1 (MQ=255)
                                                                         aCGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCgg                                                                         <  1:2021169/150‑1 (MQ=255)
                                                                              cgcgACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGgccgcc                                                                    >  2:1259617/1‑150 (MQ=255)
                                                                                 gACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTc                                                                 <  2:2143733/150‑1 (MQ=255)
                                                                                 gACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTc                                                                 <  1:2126658/150‑1 (MQ=255)
                                                                                 gACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCATTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCNTTc                                                                 <  2:1507536/150‑1 (MQ=255)
                                                                                  aCGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTcg                                                                >  1:1536558/1‑150 (MQ=255)
                                                                                     tgctgcGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCgcgc                                                             >  2:532775/1‑150 (MQ=255)
                                                                                      gctgcGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGcc                                                            >  1:597399/1‑150 (MQ=255)
                                                                                        tgcGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGt                                                          <  1:1259617/150‑1 (MQ=255)
                                                                                        tgcGCCGCATCCCCGGCGCGCAAGTGCAGGACTCCAACGGGACAGGTAGTAAGACAAGACATTTATTGAGTGGCATCCGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGt                                                          <  1:2112723/150‑1 (MQ=255)
                                                                                         gcgcGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTgg                                                                                                                                        <  2:2116842/71‑1 (MQ=255)
                                                                                         gcgcGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTgg                                                                                                                                        >  1:2116842/1‑71 (MQ=255)
                                                                                         gcgcGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTa                                                         >  1:1201685/1‑150 (MQ=255)
                                                                                           gcgGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGAAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATg                                                       <  2:1485358/150‑1 (MQ=255)
                                                                                             ggCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGgc                                                     >  2:2274478/1‑150 (MQ=255)
                                                                                                 atacacGGCGTGAAAGTTCAAGAAGCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGAGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGgccag                                                  >  2:279412/6‑150 (MQ=255)
                                                                                                    cccGGCGTGCAATTGCAGGACTCCAACGGGAGCGGTGGGAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTATCGCCGCGCTCCACGGCACTGAACCATGGCAGCCCGGCCGCCTTCGCGCCGTATGgccagcca                                               >  2:534414/1‑150 (MQ=255)
                                                                                                    cccGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGgccagcca                                               >  1:267660/1‑150 (MQ=255)
                                                                                                     ccGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCcac                                              >  1:1573783/1‑150 (MQ=255)
                                                                                                     ccGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCcac                                              >  1:2045392/1‑150 (MQ=255)
                                                                                                      cGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATGCCGGCCGCCTTCGCGCCGTATGGCCAGCcaca                                             >  1:1819534/1‑150 (MQ=255)
                                                                                                      cGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATGACGATCGCCTTCGCGCCGTATGGCCAGCcaca                                             >  1:1819871/1‑150 (MQ=255)
                                                                                                      cGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCcaca                                             <  2:2348807/150‑1 (MQ=255)
                                                                                                         cGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCt                                          <  2:1483896/150‑1 (MQ=255)
                                                                                                              aaGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGa                                     <  2:561665/150‑1 (MQ=255)
                                                                                                               aGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGAt                                    <  2:1201685/150‑1 (MQ=255)
                                                                                                               aGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGAt                                    <  1:2031256/150‑1 (MQ=255)
                                                                                                                 tGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATgg                                  <  1:820688/150‑1 (MQ=255)
                                                                                                                   cacgACCCCAACGGGACCGGTGGCAGCGACCTGTCGCTCAACTTGGGCGAACGCGGCCTTACTTCGCCCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCTGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCa                                <  1:778124/147‑1 (MQ=255)
                                                                                                                    aGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCAc                               <  2:1536558/150‑1 (MQ=255)
                                                                                                                      gACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTg                                                                                                                                          <  2:1582443/42‑1 (MQ=255)
                                                                                                                      gACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTg                                                                                                                                          >  1:1582443/1‑42 (MQ=255)
                                                                                                                       aCTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGa                            >  2:1581996/1‑150 (MQ=255)
                                                                                                                       aCTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGAATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGa                            <  1:2274478/150‑1 (MQ=255)
                                                                                                                        cTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGCTGGCACCGAt                           >  1:1072217/1‑150 (MQ=255)
                                                                                                                         tCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTcg                                                                >  1:439616/1‑112 (MQ=255)
                                                                                                                         tCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTcg                                                                <  2:439616/112‑1 (MQ=255)
                                                                                                                         tCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGAtt                          >  1:2350760/1‑150 (MQ=255)
                                                                                                                         tCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGAtt                          >  2:1869943/1‑150 (MQ=255)
                                                                                                                           cAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCAATGGCACCGATttc                        >  1:518538/1‑150 (MQ=255)
                                                                                                                           cAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCACGGCACGGGACTGATCGATGGCATCCCGAAAGCCTTCACGCCGTATGGCCAGCCACAGCTGTCGATGCCACCGATttc                        >  1:533773/1‑150 (MQ=255)
                                                                                                                              cGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATttcttc                     <  2:597399/150‑1 (MQ=255)
                                                                                                                              cGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATttcttc                     <  1:1211925/150‑1 (MQ=255)
                                                                                                                              cGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATttcttc                     <  2:267660/150‑1 (MQ=255)
                                                                                                                              cGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATttcttc                     <  2:2074724/150‑1 (MQ=255)
                                                                                                                               gggACCGGTGGCTTCGACCTGTCGCTGAACGTGGGCGTCCCCTGCCGTACTGCAAGTTTTCCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATTtcttct                    >  2:1263911/1‑150 (MQ=255)
                                                                                                                               gggACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCACCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATTtcttct                    <  2:2432649/150‑1 (MQ=255)
                                                                                                                                ggACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATTTCTTCTg                   <  2:354379/150‑1 (MQ=255)
                                                                                                                                 gACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATTtcttct                    <  1:1684791/148‑1 (MQ=255)
                                                                                                                                 gACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATTtcttct                    >  2:1684791/1‑148 (MQ=255)
                                                                                                                                    cGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATgg                                  >  2:2294285/1‑131 (MQ=255)
                                                                                                                                    cGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATgg                                  <  1:2294285/131‑1 (MQ=255)
                                                                                                                                         gCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATTTCTTCTGGCAACCTCg          >  1:766555/1‑150 (MQ=255)
                                                                                                                                          cAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCCTCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATTTCTTCTGGCAACCTCGa         <  1:1123603/150‑1 (MQ=255)
                                                                                                                                          cAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATTTCTTCTGGCAACCTCGa         >  2:380879/1‑150 (MQ=255)
                                                                                                                                                 cTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGAt                                    >  1:972373/1‑116 (MQ=255)
                                                                                                                                                 cTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGAt                                    <  2:972373/116‑1 (MQ=255)
                                                                                                                                                 cTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATTTCTTCTGGCAACCTCGACAGCATc  <  1:2294215/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
AGCGTGGTACAGAACCACCCTGGGGCACGCACTGTGGTACGCCGCGAGGCAATGGTGGAGAAGGGCGCAATGAACGTGCGCGACGTGCTGCGCGGCATCCCCGGCGTGCAAGTGCAGGACTCCAACGGGACCGGTGGCAGCGACCTGTCGCTGAACGTGGGCGTACGCGGCCTTACTTCGCGCCTGTCGCCGCGCTCCACGGTACTGATCGATGGCATCCCGGCCGCCTTCGCGCCGTATGGCCAGCCACAGCTGTCGATGGCACCGATTTCTTCTGGCAACCTCGACAGCATC  >  NC_002947/1005479‑1005772

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