Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,291,289 T→A 7.7% L280* (TTG→TAG)  PP_4664 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,291,2890TA 7.7% 247.1 / 8.8 80L280* (TTG→TAG) PP_4664LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (35/38);  new base A (4/2);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

ATGGTCGCGGGCGGGGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGCCTATG  >  NC_002947/5291144‑5291436
                                                                                                                                                 |                                                                                                                                                   
aTGGTCGCGGGCGGGGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCg                                                                                                                                                 <  2:203937/150‑1 (MQ=255)
 tGGTCGCGGGCGGGGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCgg                                                                                                                                                <  2:1242569/150‑1 (MQ=255)
            ggggCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCa                                                                                                                                     <  2:1166110/150‑1 (MQ=255)
               gCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGaa                                                                                                                                  >  1:33675/1‑150 (MQ=255)
                 cGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACt                                                                                                                                >  2:548902/1‑150 (MQ=255)
                  ggAATCGGGGTGGTGGCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTg                                                                                                                               >  1:696773/1‑150 (MQ=255)
                  ggAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTg                                                                                                                               >  1:721471/1‑150 (MQ=255)
                   gAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTgg                                                                                                                              <  2:721471/150‑1 (MQ=255)
                   gAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATCGGCCACCCGCCGGTTGCTGGTGTGCGCACGGCATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTgg                                                                                                                              <  2:797107/150‑1 (MQ=255)
                         gggTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCa                                                                                                                        >  1:756000/1‑150 (MQ=255)
                            tggtgCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCt                                                                                                                     <  2:756000/150‑1 (MQ=255)
                              gtgCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCttt                                                                                                                   >  2:459858/1‑150 (MQ=255)
                               tgCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTa                                                                                                                  >  2:1411037/1‑150 (MQ=255)
                                    cTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTAcc                                                                                                                >  1:946793/1‑147 (MQ=255)
                                    cTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTAcc                                                                                                                <  2:946793/147‑1 (MQ=255)
                                         agcagTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGAcg                                                                                                        >  2:681863/1‑150 (MQ=255)
                                          gcagTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGAcgc                                                                                                       <  2:696773/150‑1 (MQ=255)
                                            agTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGc                                                                                                             >  2:700083/1‑142 (MQ=255)
                                            agTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGc                                                                                                             <  1:700083/142‑1 (MQ=255)
                                              tACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGCGCTCAACCATTTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCaa                                                                                                   >  2:37211/1‑150 (MQ=255)
                                                     cTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGca                                                                                            <  1:1099970/150‑1 (MQ=255)
                                                        gCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGTCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCacac                                                                                         >  2:116632/1‑150 (MQ=255)
                                                        gCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAgg                                                                                                 >  2:1295282/1‑142 (MQ=255)
                                                        gCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAgg                                                                                                 <  1:1295282/142‑1 (MQ=255)
                                                            ggCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACt                                                                                     <  1:548902/150‑1 (MQ=255)
                                                               aTGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGt                                                                                                                             <  1:1424772/107‑1 (MQ=255)
                                                               aTGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGt                                                                                                                             >  2:1424772/1‑107 (MQ=255)
                                                                 gctgcCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGa                                                                                >  2:685601/1‑150 (MQ=255)
                                                                 gctgcCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGa                                                                                <  2:33675/150‑1 (MQ=255)
                                                                  ctgcCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAg                                                                               <  1:1265717/150‑1 (MQ=255)
                                                                   tgcCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACt                                                                                     >  1:649771/1‑143 (MQ=255)
                                                                   tgcCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACt                                                                                     <  2:649771/143‑1 (MQ=255)
                                                                   tgcCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAgg                                                                              <  1:1214220/150‑1 (MQ=255)
                                                                   tgcCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAgg                                                                              <  1:47994/150‑1 (MQ=255)
                                                                   tgcCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAgg                                                                              >  2:1113942/1‑150 (MQ=255)
                                                                   tgcCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGCCTGACGCGCAAGGGCGCACACCAATGTCGAgg                                                                              >  1:296544/1‑150 (MQ=255)
                                                                      cGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCaa                                                                           >  2:1402140/1‑150 (MQ=255)
                                                                       gTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGTCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCaaa                                                                          <  1:116632/150‑1 (MQ=255)
                                                                        ttGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCa                                                                                                                                     >  2:205736/1‑90 (MQ=255)
                                                                        ttGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCa                                                                                                                                     <  1:205736/90‑1 (MQ=255)
                                                                        ttGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAg                                                                               <  2:876267/144‑1 (MQ=255)
                                                                        ttGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAg                                                                               >  1:876267/1‑144 (MQ=255)
                                                                            aCTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGcc                                                                     >  1:196559/1‑150 (MQ=255)
                                                                            aCTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATAGCCCGGCAATGcc                                                                                                                                        >  2:866943/1‑83 (MQ=39)
                                                                            aCTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATAGCCCGGCAATGcc                                                                                                                                        <  1:866943/83‑1 (MQ=39)
                                                                                  tGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCAt                                                               >  2:1102987/1‑150 (MQ=255)
                                                                                      aCTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAgcg                                                           >  1:493034/1‑150 (MQ=255)
                                                                                      aCTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAgcg                                                           >  2:586511/1‑150 (MQ=255)
                                                                                       cTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAgcgc                                                          <  1:586511/150‑1 (MQ=255)
                                                                                        agCAGGAGCCATGGGCCAGCCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGcgcg                                                         <  1:21463/149‑1 (MQ=255)
                                                                                            ggAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCa                                                     <  2:196559/150‑1 (MQ=255)
                                                                                             gAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCCATTTGGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAg                                                    >  1:1224939/1‑150 (MQ=255)
                                                                                                 cATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGacac                                                >  2:1129915/1‑150 (MQ=255)
                                                                                                    gggCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGAGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTg                                             <  1:685601/150‑1 (MQ=255)
                                                                                                     cagcAACCGGCGGTTGCTGGTGTGCGCACGGGATTTTGCGGACTTGCTGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTgg                                            >  1:1097577/4‑150 (MQ=255)
                                                                                                        cAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAgg                                                                              >  2:805700/1‑113 (MQ=255)
                                                                                                        cAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAgg                                                                              <  1:805700/113‑1 (MQ=255)
                                                                                                        cAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACg                                         >  1:512059/1‑150 (MQ=255)
                                                                                                          aCCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACAcc                                               <  1:455863/142‑1 (MQ=255)
                                                                                                          aCCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACAcc                                               >  2:455863/1‑142 (MQ=255)
                                                                                                            cgccgGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGc                                                                      >  2:413093/1‑117 (MQ=255)
                                                                                                            cgccgGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGc                                                                      <  1:413093/117‑1 (MQ=255)
                                                                                                            cgccgGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGa                                     <  1:789447/150‑1 (MQ=255)
                                                                                                            cgccgGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGa                                     <  1:1113942/150‑1 (MQ=255)
                                                                                                            cgccgGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGa                                     <  1:222664/150‑1 (MQ=255)
                                                                                                             gccgGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAg                                    <  1:1411037/150‑1 (MQ=255)
                                                                                                                        gtgtGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACt                         >  1:171301/1‑150 (MQ=255)
                                                                                                                         tgtgCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGc                                                                      >  1:604792/1‑104 (MQ=255)
                                                                                                                         tgtgCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGc                                                                      <  2:604792/104‑1 (MQ=255)
                                                                                                                         tgtgCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTg                        >  1:215393/1‑150 (MQ=255)
                                                                                                                           tgCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTgcg                      <  2:493034/150‑1 (MQ=255)
                                                                                                                              gcACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGcaca                                                                                          <  1:560242/79‑1 (MQ=255)
                                                                                                                              gcACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGcaca                                                                                          >  2:560242/1‑79 (MQ=255)
                                                                                                                                        ttGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTgg         <  2:1097577/150‑1 (MQ=255)
                                                                                                                                         tGCGGCATAGCCCGGCAATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTggg        >  1:105406/1‑150 (MQ=255)
                                                                                                                                          gCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGc       <  1:687346/150‑1 (MQ=255)
                                                                                                                                          gCGGCATAGCCCGGCAATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTgc                       >  1:927605/1‑134 (MQ=255)
                                                                                                                                          gCGGCATAGCCCGGCAATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTgc                       <  2:927605/134‑1 (MQ=255)
                                                                                                                                            ggCATAGCCCGGCAATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGCCt     >  1:677076/1‑150 (MQ=255)
                                                                                                                                               aTTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGCCTATg  <  1:1102987/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ATGGTCGCGGGCGGGGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGCCTATG  >  NC_002947/5291144‑5291436

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