Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,305,153 A→T 6.0% L354Q (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,1530AT 6.0% 209.2 / 3.2 66L354Q (CTG→CAG) PP_2906sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (27/35);  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 = 9.96e-01

TCTGCTTCGGGAATACCAGGCCCGCTATCCTCGATCGCCATGTGTACGCCGCTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCATGCAGTTCTACCAGCAGCTCATGCGGGTCAAGCTGGATGAAGCCCGCGCG  >  NC_002947/3305009‑3305295
                                                                                                                                                |                                                                                                                                              
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                                                                                                                                        >  2:1144521/1‑85 (MQ=255)
                                                                    ggTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCgg                                                                                                                                        <  1:1144521/85‑1 (MQ=255)
                                                                    ggTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGc                                                                                                                          >  2:1387240/1‑99 (MQ=255)
                                                                    ggTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGc                                                                                                                          <  1:1387240/99‑1 (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                          <  2:840965/146‑1 (MQ=255)
                                                                                                                     accaGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCATGCAGTTCTACCAGCAGCTCATGCgg                          >  1:840965/1‑146 (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)
                                                                                                                                              gcagTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCAt                                                   >  2:619843/1‑96 (MQ=255)
                                                                                                                                              gcagTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCAt                                                   <  1:619843/96‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
TCTGCTTCGGGAATACCAGGCCCGCTATCCTCGATCGCCATGTGTACGCCGCTGTCATCCTGCGTGGCGGTGATGCGCACCCGACCACCCTCGGGCGTGAACTTGATCGCGTTACCCACCAGATTCGCCAAGGCCTCGAACAGCAGTTCGCGGTCGCCGTGCAACGCCGGTAGCTGCGCAGGTTGGTGCAATTGCAGTTGAATGCCGCCATCCTCGGCCAATGGCAGGTAAAAATCATGCAGTTCTACCAGCAGCTCATGCGGGTCAAGCTGGATGAAGCCCGCGCG  >  NC_002947/3305009‑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: ‑