Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,335,729 T→A 6.8% L704Q (CTG→CAG)  PP_4695 → two‑component system sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,335,7290TA 6.8% 419.8 / 14.5 130L704Q (CTG→CAG) PP_4695two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (60/61);  new base A (2/7);  total (62/68)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 131 total aligned reads displayed.

CCTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGACCTTACCGAAACCCAGGCACTGGAAGACAAGCTGGTGCACTCCGAGCGC  >  NC_002947/5335580‑5335877
                                                                                                                                                     |                                                                                                                                                    
ccTGATGTGGAACAAGGCCATGCAAGACCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCt                                                                                                                                                      <  2:1507996/150‑1 (MQ=255)
 cTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTGCAAAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTg                                                                                                                                                     >  1:737387/1‑150 (MQ=255)
   gATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGgc                                                                                                                                                   <  2:800482/150‑1 (MQ=255)
   gATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCATGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCAGgc                                                                                                                                                   <  2:1182655/150‑1 (MQ=255)
     tgtgGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGgcgc                                                                                                                                                 <  2:2402104/150‑1 (MQ=255)
      gtgGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCGTTGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCt                                                                                                                                                >  1:781024/1‑150 (MQ=255)
         gAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGa                                                                                                                                             >  2:856708/1‑150 (MQ=255)
          aaCAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCAGGCGCTgctt                                                                                                                                            <  2:1007417/150‑4 (MQ=255)
           aCAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTg                                                                                                                                                     >  2:2212083/1‑140 (MQ=255)
           aCAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTg                                                                                                                                                     <  1:2212083/140‑1 (MQ=255)
           aCAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACg                                                                                                                                           >  1:249429/1‑150 (MQ=255)
           aCAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCAGGCGCTgcttg                                                                                                                                           <  1:1205227/150‑5 (MQ=255)
               ggCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCAGGCGCTgctt                                                                                                                                            >  2:2407315/1‑142 (MQ=255)
               ggCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCAGGCGCTgctt                                                                                                                                            <  1:2407315/145‑4 (MQ=255)
                 ccATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGc                                                                                                                                     <  2:1434073/150‑1 (MQ=255)
                  cATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCCTCGACGAGCCTTGGCGCGGCCTGCTGCACGGCTTCATCAACGTGCCCGACGAGCACCTGCACCAGCAGCGCCTGGCGCTGGACGGCCAGcc                                                                                                                                    >  1:469357/1‑150 (MQ=255)
                  cATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGcc                                                                                                                                    >  1:514820/1‑150 (MQ=255)
                     ggAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCcgcg                                                                                                                                 >  1:441442/1‑150 (MQ=255)
                      gAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGGGGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCgcgc                                                                                                                                >  1:1999626/1‑150 (MQ=255)
                       aaGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCgct                                                                                                                               <  1:2214507/150‑1 (MQ=255)
                        agagCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCgctg                                                                                                                              <  1:856708/150‑1 (MQ=255)
                           gCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGGCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCgctggct                                                                                                                           <  2:1999626/150‑1 (MQ=255)
                               aCCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAAc                                                                                                                       >  1:317859/1‑150 (MQ=255)
                               aCCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAAc                                                                                                                       <  2:1107280/150‑1 (MQ=255)
                               aCCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAAc                                                                                                                       >  2:2148012/1‑150 (MQ=255)
                                 cGGCATCGCCGCCACGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTAATCAAAGGACCCGACGAGCACCTGCACAAGCAGCGCCTAGCGCTGGACGGCCAGCCGCGCTGTCTGAACCt                                                                                                                     >  2:686601/1‑150 (MQ=255)
                                 cGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCt                                                                                                                     >  2:209513/1‑150 (MQ=255)
                                       cgccgcCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACaa                                                                                                               <  2:441442/150‑1 (MQ=255)
                                       cgccgcCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACaa                                                                                                               <  2:249429/150‑1 (MQ=255)
                                        gccgccAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAg                                                                                                              >  1:627080/1‑150 (MQ=255)
                                             cAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGc                                                                                                         >  1:1033095/1‑150 (MQ=255)
                                             cAAGCATGTGGTCGGCGCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAgg                                                                                                         <  2:1658607/150‑2 (MQ=255)
                                                gCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCAt                                                                                                      >  2:625839/1‑150 (MQ=255)
                                                 cATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATc                                                                                                     <  2:1195412/150‑1 (MQ=255)
                                                      ggTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATcgacga                                                                                                <  2:457772/150‑1 (MQ=255)
                                                        tCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGc                                                                                              <  1:1491502/150‑1 (MQ=255)
                                                         cGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACg                                                                                                                                           <  1:1181441/104‑1 (MQ=255)
                                                         cGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACg                                                                                                                                           >  2:1181441/1‑104 (MQ=255)
                                                            cTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAATGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCAtt                                                                                          <  1:2368165/150‑1 (MQ=255)
                                                            cTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCAtt                                                                                          <  1:625839/150‑1 (MQ=255)
                                                             tCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAAGGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTg                                                                                         >  2:1566022/1‑150 (MQ=255)
                                                             tCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTg                                                                                         >  2:712280/1‑150 (MQ=255)
                                                                 gcCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCGTCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCa                                                                                                           <  1:1527878/128‑1 (MQ=255)
                                                                  ccTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGcccc                                                                                    <  1:1610027/150‑1 (MQ=255)
                                                                   cTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGccccc                                                                                   <  1:1760931/150‑1 (MQ=255)
                                                                    tGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCg                                                                                  >  1:569307/1‑150 (MQ=255)
                                                                     ggTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCgg                                                                                 >  1:2225292/1‑150 (MQ=255)
                                                                       tCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCATCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTa                                                                               >  1:82956/1‑150 (MQ=255)
                                                                           cATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAgg                                                                                                             <  2:805758/116‑1 (MQ=255)
                                                                           cATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAgg                                                                                                             >  1:805758/1‑116 (MQ=255)
                                                                             tCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCg                                                                         <  2:1033095/150‑1 (MQ=255)
                                                                               gacgaGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATcgac                                                                                                  <  1:641169/123‑1 (MQ=255)
                                                                                acgaGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGc                                                                      >  1:1679584/1‑150 (MQ=255)
                                                                                acgaGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCAGGc                                                                      <  2:569307/150‑1 (MQ=255)
                                                                                 gaaGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCt                                                                     >  1:1070161/3‑150 (MQ=255)
                                                                                    gCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGc                                                                                                                   <  2:1243487/101‑1 (MQ=255)
                                                                                    gCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGc                                                                                                                   >  1:1243487/1‑101 (MQ=255)
                                                                                        tGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTActg                                                              <  2:317859/150‑1 (MQ=255)
                                                                                         ggCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTActgc                                                             <  2:475282/150‑1 (MQ=255)
                                                                                           gacgGCCTGCGGCAAGGTTGCATGAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTActgctg                                                           <  2:1342937/148‑1 (MQ=255)
                                                                                             cgGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCa                                                                                                           <  1:297092/100‑1 (MQ=255)
                                                                                             cgGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCtggt                                                         <  2:2161462/150‑1 (MQ=255)
                                                                                                   gctgcAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGa                                                   <  1:209513/150‑1 (MQ=255)
                                                                                                     tgcAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGAcc                                                 >  2:2247842/1‑150 (MQ=255)
                                                                                                      gcAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGACCt                                                >  2:1853207/1‑150 (MQ=255)
                                                                                                      gcAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGACCt                                                >  2:2018354/1‑150 (MQ=255)
                                                                                                       cAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCgg                                                                        >  2:1778551/1‑125 (MQ=255)
                                                                                                       cAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCAGTAGGAAGACCtt                                               >  1:2175915/1‑150 (MQ=255)
                                                                                                        aaGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCAtt                                                                                          <  2:599304/106‑1 (MQ=255)
                                                                                                        aaGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCAtt                                                                                          >  1:599304/1‑106 (MQ=255)
                                                                                                        aaGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTgg                                                                                        >  1:1973707/1‑108 (MQ=255)
                                                                                                        aaGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTgg                                                                                        <  2:1973707/108‑1 (MQ=255)
                                                                                                           gCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGACCt                                                >  2:78371/1‑145 (MQ=255)
                                                                                                           gCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGACCt                                                <  1:78371/145‑1 (MQ=255)
                                                                                                               catcaACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGACCTTACCGAAAc                                       <  2:1283017/150‑1 (MQ=255)
                                                                                                                atcaACGTGCCCGACGAGCACCTGCACCAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGACCTTACCGAAAcc                                      >  1:676708/1‑150 (MQ=255)
                                                                                                                 tcaACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGACCTTACCGAAAccc                                     >  1:1566141/1‑150 (MQ=255)
                                                                                                                  caACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGACCTTACCGAAACCCa                                    <  1:312381/150‑1 (MQ=255)
                                                                                                                         cccGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGACCTTACCGAAACCCAGGCACTg                             >  2:869278/1‑150 (MQ=255)
                                                                                                                           cgatgaGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGACAGCCGCGCTGGCTGAACCTGCACAAAGCAGCCATCGACGAGCCATTGGCCCCCGTTAACAGCGGGCCGGTACTGGTGGTGGAAGACCTGACCGAAACCCAGGCACTggc                           >  2:881741/1‑149 (MQ=255)
                                                                                                                           cgacgaGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGACCTTACCGAAACCCAGGCACTGGa                           <  2:1522825/150‑1 (MQ=255)
                                                                                                                            gacgaGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGACCTTACCGAAACCCAGGCACTGGaa                          <  1:1145652/150‑1 (MQ=255)
                                                                                                                            gacgaGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGACCTTACCGAAACCCAGGCACTGGaa                          <  1:1566022/150‑1 (MQ=255)
                                                                                                                            gacgaGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGACCTTACCGAAACCCAGGCACTGGaa                          >  1:1509233/1‑150 (MQ=255)
                                                                                                                                  cACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGACCTTACCGAAACCCa                                    >  1:1275118/1‑134 (MQ=255)
                                                                                                                                    ccTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGACCTTACCGAAACCCAGGCACTGGAAGACAAGCt                  >  1:1349284/1‑150 (MQ=255)
                                                                                                                                    ccTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGACCTTACCGAAACCCAGGCACTGGAAGACAAGCt                  <  1:1359827/150‑1 (MQ=255)
                                                                                                                                     cTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAg                                                                                                              >  1:2042744/1‑57 (MQ=255)
                                                                                                                                     cTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACa                                                                            <  2:2009806/91‑1 (MQ=255)
                                                                                                                                     cTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACa                                                                            >  1:2009806/1‑91 (MQ=255)
                                                                                                                                     cTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGACCTTACCGAAACCCAGGCACTGGAAGACAAGCTg                 <  2:1679584/150‑1 (MQ=255)
                                                                                                                                      tGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGACCTTACCGAAACCCAGGCACTGGAAGACAAGCTgg                <  2:1349284/150‑1 (MQ=255)
                                                                                                                                       gCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGACCTTACCGAAACCCAGGCACTGGAAGACAAGCTGGt               <  2:648124/150‑1 (MQ=255)
                                                                                                                                       gCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGACCTTACCGAAACCCAGGCACTGGAAGACAAGCTGGt               <  1:387330/150‑1 (MQ=255)
                                                                                                                                           aaGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGaa                                                     >  2:1640740/1‑108 (MQ=255)
                                                                                                                                           aaGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGACCTTACCGAAACCCAGGCACTGGAAGACAAGCTGGTGCAc           >  1:224954/1‑150 (MQ=255)
                                                                                                                                                gcgcCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGACCTTACCGAAACCCAGGCACTGGAAGACAAGCTGGTGCACTCCGa      >  1:1284476/1‑150 (MQ=255)
                                                                                                                                                    cTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGACCTTACCGAAACCCAGGCACTGGAAGACAAGCTGGTGCACTCCGAgcgc  >  1:1699505/1‑150 (MQ=255)
                                                                                                                                                    cTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGACCTTACCGAAACCCAGGCACTGGAAGACAAGCTGGTGCACTCCGAgcgc  >  2:2242455/1‑150 (MQ=255)
                                                                                                                                                    cTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGCACTGCTGGTGGAAGACCTTACCGAAACCCAGGCACTGGAAGACAAGCTGGTGCACTCCGAgcgc  >  1:1491717/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAAGCAGCGCCTGGCGCTGGACGGCCAGCCGCGCTGGCTGAACCTGCACAAGGCAGCCATCGACGAGCCATTGGCCCCCGGTAACAGCGGGCTGGTACTGCTGGTGGAAGACCTTACCGAAACCCAGGCACTGGAAGACAAGCTGGTGCACTCCGAGCGC  >  NC_002947/5335580‑5335877

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: ‑