Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,305,150 A→T 5.8% L355Q (CTG→CAG)  PP_2906 ← sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,305,1500AT 5.8% 206.4 / 2.4 66L355Q (CTG→CAG) PP_2906sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (27/34);  new base T (2/2);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTTCGCTCTGCTTCGGGAATACCAGGCCCGCTATCCTCGATCGCCATGTGTACGCCGCTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCATGCAGTTCTACCAGCAGCTCATGCGGGTCAAGCTGGATGAAGCCCGCGCG  >  NC_002947/3305001‑3305295
                                                                                                                                                     |                                                                                                                                                 
ccGTTCGCTCTGCTTCGGGAATACCAGGCCCGCTATCCTCGATCGCCATGTGTACGCCGCTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAAca                                                                                                                                                   <  1:26433/150‑1 (MQ=255)
        tctGCTTCGGGAATACCAGGCCCGCTATCCTCGATCGCCATGTGTACGCCGCTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTcg                                                                                                                                           <  2:1554354/150‑1 (MQ=255)
          tgctgcGGGAATACCAGGGCCGCTGTCCGCGATCGCCATGTGTACGCCGCTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCGGCTGCGGCAAGGCCTCGAACAGCAGTTcgcg                                                                                                                                         <  1:1276844/150‑1 (MQ=255)
          tGCTTCGGGAATACCAGGCCCGCTATCCTCGATCGCCATGTGTACGCCGCTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTcgcg                                                                                                                                         <  2:846164/150‑1 (MQ=255)
                 ggAATACCAGGCCCGCTATCCTCGATCGCCATGTGTACGCCGCTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTCCCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCg                                                                                                                                  >  1:760412/1‑150 (MQ=255)
                  gcatagaaGGACCGCTCTCCTAGCTCGCTATGTGTACGCCGCTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGt                                                                                                                                 <  2:760412/143‑1 (MQ=255)
                        cAGGCCCGCTATCCTCGATCGCCATGTGTACGCCGCTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACg                                                                                                                           >  1:50544/1‑150 (MQ=255)
                        cAGGCCCGCTATCCTCGATCGCCATGTGTACGCCGCTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACg                                                                                                                           >  1:214874/1‑150 (MQ=255)
                         aGGCCCGCTATCCTCGATCGCCATGTGTACGCCGCTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACTGCTGTTcgaggc                                                                                                                                       >  1:1449031/1‑133 (MQ=255)
                         aGGCCCGCTATCCTCGATCGCCATGTGTACGCCGCTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACTGCTGTTCGAggc                                                                                                                                       <  2:1449031/137‑2 (MQ=255)
                         aGGCCCGCTATCCTCGATCGCCATGTATACGCCGCTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGc                                                                                                                          >  1:808714/1‑150 (MQ=255)
                                 tATCCTCGCTCGCCATGTGTACGCCGCTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCt                                                                                                                  <  1:63865/150‑1 (MQ=255)
                                      tCGATCGCCATGTGTACGCCGCTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTTGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCa                                                                                                             >  1:134951/1‑150 (MQ=255)
                                      tCGATCGCCATGTGTACGCCGCTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACTGCTGTTCGAggccttgg                                                                                                                                  <  1:45217/129‑7 (MQ=255)
                                      tCGATCGCCATGTGTACGCCGCTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACACTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACTGCTGTTCGAggccttgg                                                                                                                                  >  2:45217/1‑123 (MQ=255)
                                        gATCGCCATGTGTACGCCGCTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAgg                                                                                                           >  2:302838/1‑150 (MQ=255)
                                           cGCCATGTGTACGCCGCTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTg                                                                                                        <  2:243035/150‑1 (MQ=255)
                                             ccATGTGTACGCCGCTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGt                                                                                                      >  1:842250/1‑150 (MQ=255)
                                             ccATGTGTACGCCGCTGTCATCCTGCGTGGCCGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCCGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGt                                                                                                      <  2:336396/150‑1 (MQ=255)
                                              cATGTGTACGCCGCTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTg                                                                                                     >  1:1154392/1‑150 (MQ=255)
                                                 gtatACGCCGCTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCaa                                                                                                  <  2:808714/147‑1 (MQ=255)
                                                  tgtACGCCGCTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCCCCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGCAACGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAAt                                                                                                 <  1:1275342/150‑1 (MQ=255)
                                                    tACGCCGCTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTg                                                                                               >  1:826927/1‑150 (MQ=255)
                                                       gccgcTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAg                                                                                            >  2:345341/1‑150 (MQ=255)
                                                       gccgcTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGAGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAg                                                                                            >  2:1625098/1‑150 (MQ=255)
                                                            tGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGaa                                                                                       <  1:791222/150‑1 (MQ=255)
                                                               cATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTTGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATgc                                                                                    <  2:134951/150‑1 (MQ=255)
                                                                   cTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTg                                                                                                     >  1:511927/1‑129 (MQ=255)
                                                                   cTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTg                                                                                                     <  2:511927/129‑1 (MQ=255)
                                                                     gCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGtt                                                                                          >  1:129725/1‑138 (MQ=255)
                                                                     gCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGtt                                                                                          <  2:129725/138‑1 (MQ=255)
                                                                           cGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCgg                                                                        <  2:394356/150‑1 (MQ=255)
                                                                            ggTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCgg                                                                                                                                        <  1:1144521/85‑1 (MQ=255)
                                                                            ggTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCgg                                                                                                                                        >  2:1144521/1‑85 (MQ=255)
                                                                            ggTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGc                                                                                                                          <  1:1387240/99‑1 (MQ=255)
                                                                            ggTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGc                                                                                                                          >  2:1387240/1‑99 (MQ=255)
                                                                               gATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCaa                                                                    <  1:1003148/150‑1 (MQ=255)
                                                                                  gcgcACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATgg                                                                 >  1:495847/1‑150 (MQ=255)
                                                                                    gcACCCGACCACTCTCGTGCGTGAACTTGATCGCTTTACCCACAAGATTTGCCAAGTCCTCGTACAGCAGTCCGCGGTCGCCGTGCCATGGCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTGCTTTGCCGCCAGCCTCGGCCAATGGCa                                                               <  2:605222/150‑1 (MQ=255)
                                                                                    gcACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCa                                                               <  1:597210/150‑1 (MQ=255)
                                                                                    ccACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCa                                                               <  1:189032/149‑1 (MQ=255)
                                                                                     cACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAg                                                              <  1:1244563/150‑1 (MQ=255)
                                                                                      aCCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAgg                                                             >  2:1164646/1‑150 (MQ=255)
                                                                                       cccGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGt                                                            >  1:1230718/1‑150 (MQ=255)
                                                                                               cccTCGGGCGTGAGCTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCa                                                    <  2:598791/150‑1 (MQ=255)
                                                                                                ccTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCAt                                                   >  2:1161951/1‑150 (MQ=255)
                                                                                                 cTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGCATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCATg                                                  <  2:842250/150‑1 (MQ=255)
                                                                                                             cTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCATGCAg                                               <  2:1047976/141‑1 (MQ=255)
                                                                                                             cTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCATGCAg                                               >  1:1047976/1‑141 (MQ=255)
                                                                                                                  tCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCATGCAGTTCTACCAGCAGCt                                 >  1:988688/1‑150 (MQ=255)
                                                                                                                          cccACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCATGCAGtt                                             <  2:963704/130‑1 (MQ=255)
                                                                                                                          cccACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCATGCAGtt                                             >  1:963704/1‑130 (MQ=255)
                                                                                                                            caccaGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCATGCAg                                               <  2:389920/126‑1 (MQ=255)
                                                                                                                            caccaGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCATGCAg                                               >  1:389920/1‑126 (MQ=255)
                                                                                                                            caccaGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCATGCAGTTCTACCAGCAGCTCATGCGGGTc                       >  2:261864/1‑150 (MQ=255)
                                                                                                                             accaGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGtt                                                                                                         <  2:470455/67‑1 (MQ=255)
                                                                                                                             accaGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGtt                                                                                                         >  1:470455/1‑67 (MQ=255)
                                                                                                                             accaGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCATGCAGTTCTACCAGCAGCTCATGCgg                          >  1:840965/1‑146 (MQ=255)
                                                                                                                             accaGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCATGCAGTTCTACCAGCAGCTCATGCgg                          <  2:840965/146‑1 (MQ=255)
                                                                                                                                aGATTCGCCAAGGCCTCGAACAGCAGTTCGAGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCATGCAGTTCTACCAGCAGCTCATGCGGGTCAAGc                   <  1:1625098/150‑1 (MQ=255)
                                                                                                                                 gATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCATGCAGTTCTACCAGCAGCTCATGCGGGTCAAGCt                  <  1:345341/150‑1 (MQ=255)
                                                                                                                                    tCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCATGCAGTTCTACCAGCAGCTCATGCTGGTCAAGCTGGa               <  1:1161951/150‑1 (MQ=255)
                                                                                                                                    tCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCATGCAGTTCTACCAGCAGCTCATGCGGGTCAAGCTGGa               <  2:214874/150‑1 (MQ=255)
                                                                                                                                        cAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCAt                                                   <  1:884166/110‑1 (MQ=255)
                                                                                                                                        cAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCAt                                                   >  2:884166/1‑110 (MQ=255)
                                                                                                                                         aaGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAATTCATGCAGTTCTACCAGCAGCTCATGCGGGTCAAGCTGGATGAAg          <  1:1164646/150‑1 (MQ=255)
                                                                                                                                                 gAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCATGCAGTTCTACCAGCAGCTCATGCGGGTCAAGCTGGATGAAGCCcgcgcg  <  1:942218/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CCGTTCGCTCTGCTTCGGGAATACCAGGCCCGCTATCCTCGATCGCCATGTGTACGCCGCTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCATGCAGTTCTACCAGCAGCTCATGCGGGTCAAGCTGGATGAAGCCCGCGCG  >  NC_002947/3305001‑3305295

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