Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 812,056 G→A 5.0% L67L (CTG→TTG)  PP_0698 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947812,0560GA 5.0% 566.4 / 11.8 178L67L (CTG→TTG) PP_0698LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (83/86);  new base A (3/6);  total (86/92)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 182 total aligned reads displayed.

CATGGAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGCTCTCGGTC  >  NC_002947/811909‑812204
                                                                                                                                                   |                                                                                                                                                    
 aTGGAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGc                                                                                                                                                   >  2:5328525/1‑150 (MQ=255)
 aTGGAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGc                                                                                                                                                   >  2:5060051/1‑150 (MQ=255)
   ggAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCAc                                                                                                                                                 <  1:1116557/150‑1 (MQ=255)
   ggAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCAc                                                                                                                                                 <  1:1116459/150‑1 (MQ=255)
   ggAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCAc                                                                                                                                                 <  1:2254908/150‑1 (MQ=255)
    gAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACg                                                                                                                                                >  2:4185419/1‑150 (MQ=255)
    gAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTCGTGCCCTGCAAGGCAcc                                                                                                                                                <  1:1373853/150‑2 (MQ=255)
       gcgcGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGAcc                                                                                                                                             >  1:435258/1‑150 (MQ=255)
        cgcgGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTCGTGCCCTGCAAGGCaccaccc                                                                                                                                            <  1:5034441/150‑2 (MQ=255)
           ggCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCGGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACcggcg                                                                                                                                         >  2:1386093/1‑150 (MQ=255)
           ggCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACcggcg                                                                                                                                         >  2:3643890/1‑150 (MQ=255)
           ggCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTCGTGCCCTGCAAGGCaccacccagg                                                                                                                                         <  2:598174/150‑5 (MQ=255)
            gCAGGCGCGGGCTCAGCCAGCGCATGGCGCTGGTCGGTTCGCAAGACCGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGGTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTTGGTGGTGCCTTGCAGGGCACGACcggcgg                                                                                                                                        <  1:4288964/150‑1 (MQ=255)
             cAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCggg                                                                                                                                       <  2:2118487/150‑1 (MQ=255)
                  gcgGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGa                                                                                                                                  <  2:248120/150‑1 (MQ=255)
                  gcgGGATCAGCCAGCGCAGGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGa                                                                                                                                  <  2:2922444/150‑1 (MQ=255)
                   cgGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAg                                                                                                                                 <  1:1257872/150‑1 (MQ=255)
                           gccagcGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATgg                                                                                                                         >  2:1339346/1‑150 (MQ=255)
                           gccagcGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATgg                                                                                                                         >  2:3265376/1‑150 (MQ=255)
                           gccagcGCATGGCGATGGGCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATgg                                                                                                                         >  2:5372202/1‑150 (MQ=255)
                            ccagcGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGc                                                                                                                        >  1:4246920/1‑150 (MQ=255)
                               gcgcATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGTCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCg                                                                                                                     <  2:49960/150‑1 (MQ=255)
                                 gcATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGt                                                                                                                   >  1:2531499/1‑150 (MQ=255)
                                     ggCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTgcg                                                                                                               <  2:2895798/150‑1 (MQ=255)
                                     ggCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTgcg                                                                                                               >  1:1828142/1‑150 (MQ=255)
                                     ggCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTgcg                                                                                                               <  1:864070/150‑1 (MQ=255)
                                     ggCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTgcg                                                                                                               >  1:1895352/1‑150 (MQ=255)
                                      gCGATGGTCGGGTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTgcgg                                                                                                              <  2:4356433/150‑1 (MQ=255)
                                           ggTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTc                                                                                                         <  2:1811469/150‑1 (MQ=255)
                                             tCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCCGGCAAGGcaccac                                                                                                                                              >  2:3442151/1‑111 (MQ=255)
                                               ggTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGGTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAAc                                                                                                     <  1:3849875/150‑1 (MQ=255)
                                                   cGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGtt                                                                                                 >  1:595770/1‑150 (MQ=255)
                                                   cGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGtt                                                                                                 >  1:5345465/1‑150 (MQ=255)
                                                    gCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGACCAGTTc                                                                                                <  2:712693/150‑1 (MQ=255)
                                                      aaGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCggg                                                                                                                                       >  1:3278268/1‑109 (MQ=255)
                                                      aaGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCggg                                                                                                                                       <  2:3278245/109‑1 (MQ=255)
                                                      aaGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTCGTGCCCTGCAAGGCaccacccag                                                                                                                                          <  2:832899/106‑4 (MQ=39)
                                                        gACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACg                                                                                            >  1:1218769/1‑150 (MQ=255)
                                                               ccagcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGc                                                                                     >  1:3372398/1‑150 (MQ=255)
                                                                cagcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCt                                                                                    >  1:353984/1‑150 (MQ=255)
                                                                cagcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCt                                                                                    >  1:8211/1‑150 (MQ=255)
                                                                  gcacATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCttt                                                                                  >  1:2125152/1‑150 (MQ=255)
                                                                  gcacATTGTCGCGAGCTTGTTGCCGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCttt                                                                                  >  2:1260966/1‑150 (MQ=255)
                                                                   cacaTTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTttc                                                                                 <  2:4948803/150‑1 (MQ=255)
                                                                   cacaTTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTttc                                                                                 >  1:3729202/1‑150 (MQ=255)
                                                                      aTTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGGGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTttcttc                                                                              <  2:3770424/150‑1 (MQ=255)
                                                                        tGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCca                                                                            <  2:3498971/150‑1 (MQ=255)
                                                                            gcgAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCcagcag                                                                        >  1:1165810/1‑150 (MQ=255)
                                                                              gAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGc                                                                      <  1:28950/150‑1 (MQ=255)
                                                                              gAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGc                                                                      >  1:1784900/1‑150 (MQ=255)
                                                                                  ttgttgCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTAt                                                                  >  1:4032367/1‑150 (MQ=255)
                                                                                  ttgttgCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTAt                                                                  >  1:4032203/1‑150 (MQ=255)
                                                                                  ttgttgCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTAt                                                                  >  2:2054000/1‑150 (MQ=255)
                                                                                   tgttgCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTAtt                                                                 <  1:3920686/150‑1 (MQ=255)
                                                                                    gttgCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTAttt                                                                >  1:3146393/1‑150 (MQ=255)
                                                                                     ttgCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTg                                                               >  2:3991923/1‑150 (MQ=255)
                                                                                      tgctgcAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTgg                                                              <  1:2778943/150‑1 (MQ=255)
                                                                                      tgctgcAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTgg                                                              <  1:2778947/150‑1 (MQ=255)
                                                                                      tgctgcAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTgg                                                              <  1:753286/150‑1 (MQ=255)
                                                                                       gctgcAGGCGCTGCACGGCGCCTTCGAGCTGTTCGATGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTggc                                                             >  2:463344/1‑150 (MQ=255)
                                                                                       gctgcAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTggc                                                             <  2:1639119/150‑1 (MQ=255)
                                                                                        ctgcAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTggcg                                                            >  2:3375848/1‑150 (MQ=255)
                                                                                         tgcAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTggcgg                                                           >  1:452323/1‑150 (MQ=255)
                                                                                         tgcAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTggcgg                                                           >  1:4940899/1‑150 (MQ=255)
                                                                                         tgcAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTggcgg                                                           >  1:4987588/1‑150 (MQ=255)
                                                                                             ggCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGa                                                       >  1:2576158/1‑150 (MQ=255)
                                                                                                  tGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCATTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACcgcg                                                  >  1:3583856/1‑150 (MQ=255)
                                                                                                     aCGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGtt                                                                                                          >  1:2986641/1‑91 (MQ=255)
                                                                                                     aCGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTttctt                                                                               <  2:5188/118‑1 (MQ=255)
                                                                                                          gcCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCgtgagt                                          >  1:1583164/1‑150 (MQ=255)
                                                                                                          gcCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCgtgagt                                          <  1:4734980/150‑1 (MQ=255)
                                                                                                            cTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTCTTTGGCGGCTGACCGCGCCgtgagtga                                        >  1:202562/1‑150 (MQ=255)
                                                                                                            cTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCgtgagtga                                        <  2:2369756/150‑1 (MQ=255)
                                                                                                            cTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCgtgagtga                                        <  2:3702461/150‑1 (MQ=255)
                                                                                                             ttCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGAc                                       >  1:3095674/1‑150 (MQ=255)
                                                                                                                gAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTg                                    <  1:5428439/150‑1 (MQ=255)
                                                                                                                 aGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGc                                   >  2:3249237/1‑150 (MQ=255)
                                                                                                                  acagctCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCa                                  >  2:723991/6‑150 (MQ=255)
                                                                                                                    tGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAg                                <  2:2499078/150‑1 (MQ=255)
                                                                                                                      ttCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTg                              >  1:3835335/1‑150 (MQ=255)
                                                                                                                      ttCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTg                              >  1:4617193/1‑150 (MQ=255)
                                                                                                                        cGAAGATCGCCTGGGTGGTGCCTTGCAAGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCt                            >  2:141710/1‑150 (MQ=255)
                                                                                                                         gAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCcagca                                                                         >  2:4811555/1‑104 (MQ=255)
                                                                                                                         gAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCgtga                                            <  2:3131288/133‑1 (MQ=255)
                                                                                                                          aaGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTcg                          >  1:2054344/1‑150 (MQ=255)
                                                                                                                          aaGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTcg                          >  2:4168480/1‑150 (MQ=255)
                                                                                                                           aGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTcgg                         <  2:4246920/150‑1 (MQ=255)
                                                                                                                             aTCGCCTGGGTGGTGACTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTcggcg                       >  2:2059068/1‑150 (MQ=255)
                                                                                                                              tCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTc                                                                                                         <  1:3252536/67‑1 (MQ=255)
                                                                                                                                gCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGgcggcg                    >  2:4029078/1‑150 (MQ=255)
                                                                                                                                 ccTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGgcgc                   <  2:1895352/150‑1 (MQ=255)
                                                                                                                                 ccTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGgcgc                   <  2:2497039/150‑1 (MQ=255)
                                                                                                                                    gggTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCa                <  2:1229635/150‑1 (MQ=255)
                                                                                                                                         gtgCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAg           >  2:5041764/1‑150 (MQ=255)
                                                                                                                                          tgCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGc          >  1:34286/1‑150 (MQ=255)
                                                                                                                                             cTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGctct       <  2:1072414/150‑1 (MQ=255)
                                                                                                                                             cTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGctct       <  1:2942662/150‑1 (MQ=255)
                                                                                                                                             cTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGctct       <  2:1072329/150‑1 (MQ=255)
                                                                                                                                              ttGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGCtctc      <  2:4940899/150‑1 (MQ=255)
                                                                                                                                                  aGGGCACGACCGGCGGGGGTGAGGAACATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGCTCTCGGTc  >  2:1931565/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CATGGAAGCGCGGCAGGCGCGGGATCAGCCAGCGCATGGCGATGGTCGGTTCGCAAGACAGCACCAGCACATTGTCGCGAGCTTGTTGCTGCAGGCGCTGCACGGCGCCTTCGAGCTGTTCGAAGATCGCCTGGGTGGTGCCTTGCAGGGCACGACCGGCGGGGGTGAGGAAAATGGCCCGGTTGCGGCGTTCGAACAGTTCCACGCCAAGGCTTTCTTCCAGCAGGCGTATTTGGCGGCTGACCGCGCCGTGAGTGACGTGCAAGTGCTCGGCGGCGCGCACGAAGCTCTCGGTC  >  NC_002947/811909‑812204

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