Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,291,297 T→A 7.6% Y283N (TAT→AAT)  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,2970TA 7.6% 241.2 / 8.7 79Y283N (TAT→AAT) PP_4664LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (36/37);  new base A (4/2);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

GGGGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGCCTATGGTGACCCGTC  >  NC_002947/5291156‑5291446
                                                                                                                                             |                                                                                                                                                     
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                                                                                                                          <  2:946793/147‑1 (MQ=255)
                        cTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTAcc                                                                                                                          >  1:946793/1‑147 (MQ=255)
                             agcagTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGAcg                                                                                                                  >  2:681863/1‑150 (MQ=255)
                              gcagTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGAcgc                                                                                                                 <  2:696773/150‑1 (MQ=255)
                                agTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGc                                                                                                                       <  1:700083/142‑1 (MQ=255)
                                agTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGc                                                                                                                       >  2:700083/1‑142 (MQ=255)
                                  tACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGCGCTCAACCATTTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCaa                                                                                                             >  2:37211/1‑150 (MQ=255)
                                         cTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGca                                                                                                      <  1:1099970/150‑1 (MQ=255)
                                            gCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGTCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCacac                                                                                                   >  2:116632/1‑150 (MQ=255)
                                            gCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAgg                                                                                                           <  1:1295282/142‑1 (MQ=255)
                                            gCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAgg                                                                                                           >  2:1295282/1‑142 (MQ=255)
                                                ggCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACt                                                                                               <  1:548902/150‑1 (MQ=255)
                                                   aTGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGt                                                                                                                                       >  2:1424772/1‑107 (MQ=255)
                                                   aTGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGt                                                                                                                                       <  1:1424772/107‑1 (MQ=255)
                                                     gctgcCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGa                                                                                          >  2:685601/1‑150 (MQ=255)
                                                     gctgcCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGa                                                                                          <  2:33675/150‑1 (MQ=255)
                                                      ctgcCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAg                                                                                         <  1:1265717/150‑1 (MQ=255)
                                                       tgcCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACt                                                                                               <  2:649771/143‑1 (MQ=255)
                                                       tgcCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACt                                                                                               >  1:649771/1‑143 (MQ=255)
                                                       tgcCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAgg                                                                                        <  1:47994/150‑1 (MQ=255)
                                                       tgcCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAgg                                                                                        >  2:1113942/1‑150 (MQ=255)
                                                       tgcCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAgg                                                                                        <  1:1214220/150‑1 (MQ=255)
                                                       tgcCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGCCTGACGCGCAAGGGCGCACACCAATGTCGAgg                                                                                        >  1:296544/1‑150 (MQ=255)
                                                          cGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCaa                                                                                     >  2:1402140/1‑150 (MQ=255)
                                                           gTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGTCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCaaa                                                                                    <  1:116632/150‑1 (MQ=255)
                                                            ttGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCa                                                                                                                                               <  1:205736/90‑1 (MQ=255)
                                                            ttGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCa                                                                                                                                               >  2:205736/1‑90 (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                                                                                        <  1:805700/113‑1 (MQ=255)
                                                                                            cAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAgg                                                                                        >  2:805700/1‑113 (MQ=255)
                                                                                            cAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACg                                                   >  1:512059/1‑150 (MQ=255)
                                                                                              aCCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACAcc                                                         >  2:455863/1‑142 (MQ=255)
                                                                                              aCCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACAcc                                                         <  1:455863/142‑1 (MQ=255)
                                                                                                cgccgGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGc                                                                                >  2:413093/1‑117 (MQ=255)
                                                                                                cgccgGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGc                                                                                <  1:413093/117‑1 (MQ=255)
                                                                                                cgccgGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGa                                               <  1:222664/150‑1 (MQ=255)
                                                                                                cgccgGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGa                                               <  1:789447/150‑1 (MQ=255)
                                                                                                cgccgGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGa                                               <  1:1113942/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                                 <  2:927605/134‑1 (MQ=255)
                                                                                                                              gCGGCATAGCCCGGCAATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTgc                                 >  1:927605/1‑134 (MQ=255)
                                                                                                                                ggCATAGCCCGGCAATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGCCt               >  1:677076/1‑150 (MQ=255)
                                                                                                                                   aTTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGCCTATg            <  1:1102987/150‑1 (MQ=255)
                                                                                                                                      gCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGCCTATGGTg         <  1:681863/150‑1 (MQ=255)
                                                                                                                                             tATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTCTTTCATGGCCGTTGGGCCTATGGTGACCCGTc  <  1:37211/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                     
GGGGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGCCTATGGTGACCCGTC  >  NC_002947/5291156‑5291446

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: