Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,783,375 C→G 6.8% V598V (GTC→GTG kefA → intermediate mechanosensitive channel protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,783,3750CG 6.8% 391.9 / 15.1 127V598V (GTC→GTGkefAintermediate mechanosensitive channel protein
Reads supporting (aligned to +/- strand):  ref base C (57/60);  new base G (5/4);  total (62/65)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 131 total aligned reads displayed.

TTGGCCTGGCGCTGTGCAGCTATGCCTTGCAGATCGATGCTCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAGCCTTCGGCACTGGCAGACGACGTGTTGGGCATCGGCGTGGT  >  NC_002947/5783227‑5783524
                                                                                                                                                       |                                                                                                                                                     
ttGGCCTGGCGCTGTACAGCTATGCCTTGCAGATCGATGCTCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTcc                                                                                                                                                      <  1:234361/150‑1 (MQ=255)
     cTGGCGCTGTGCAGCTATGCCTTGCAGATCGATGCTCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTc                                                                                                                                                 >  1:2188974/1‑150 (MQ=255)
      tGGCGCTGTGCAGCTATGCCTTGCAGATCGATGCTCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTcc                                                                                                                                                >  2:1375044/1‑150 (MQ=255)
       ggCGCTGTGCAGCTATGCCTTGCAGATCGATGCTCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCg                                                                                                                                               <  1:2241094/150‑1 (MQ=255)
         cgcTGTGCAGCTATGCCTTGCAGATCGATGCTCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCt                                                                                                                                             <  1:98932/150‑1 (MQ=255)
           ttgtgCAGCTATGCCTTGCCGATCGATGCTCGCGGGCAAACCGCCAACTTGTGCGCCGCCCTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCCCTTCCGCTgg                                                                                                                                           <  2:1013174/149‑1 (MQ=255)
           cTGTGCAGCTATGCCTTGCAGATCGATGCTCGCGGGCAAAACGCCCACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTgg                                                                                                                                           <  1:1358975/150‑1 (MQ=255)
           cTGTGCAGCTATGCCTTGCAGATCGATGCTCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTgg                                                                                                                                           <  1:713556/150‑1 (MQ=255)
           cTGTGCAGCTATGCCTTGCAGAGCGCTGATCCGATCATATTGTCCTACCTGGGCGCTGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTCGCGTGGCGGAAGTGGACTTCCGccac                                                                                                                                           <  2:1203792/150‑4 (MQ=255)
                 aGCTATGCCTTGCAGATCGATGCTCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACaaa                                                                                                                                     >  1:547342/1‑150 (MQ=255)
                 aGCTATGCCTTGCAGATCGATGCTCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACaaa                                                                                                                                     >  1:80758/1‑150 (MQ=255)
                 aGCTATGCCTTGCAGATCGATGCTCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACaaa                                                                                                                                     >  2:2349435/1‑150 (MQ=255)
                   cTATGCCTTGCAGATCGATGCTCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGAATAACGTTTACACCGCCTACCGTATCCTCGCAGAACATGAACGCGGGGAATGCCACTTACGCTGGCACAAAcc                                                                                                                                   >  1:2253152/1‑150 (MQ=255)
                        ccTTGCAGATCGATGCTCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAgg                                                                                                                              <  2:547342/150‑1 (MQ=255)
                        ccTTGCAGATCGATGCTCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAgg                                                                                                                              >  1:2241986/1‑150 (MQ=255)
                          ttGCAGATCGATGCTCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTc                                                                                                                            <  1:579781/150‑1 (MQ=255)
                          ttGCAGATCGATGCTCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTc                                                                                                                            >  1:1767269/1‑150 (MQ=255)
                          ttGCAGATCGATGCTCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTc                                                                                                                            >  2:201126/1‑150 (MQ=255)
                          ttGCAGATCGATGCTCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGATGGCACAAACCGCAGGTc                                                                                                                            >  2:1636833/1‑150 (MQ=255)
                            gCAGATCGATGCTCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGa                                                                                                                          <  2:1363550/150‑1 (MQ=255)
                                   gATGCTCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTg                                                                                                                   >  2:2424804/1‑150 (MQ=255)
                                   gATGCTCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTg                                                                                                                   <  2:219298/150‑1 (MQ=255)
                                    aTGCTCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTGGACTTCCGc                                                                                                                                              <  2:1079640/122‑1 (MQ=255)
                                    aTGCTCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTGGACTTCCGc                                                                                                                                              >  1:1078759/1‑122 (MQ=255)
                                    aTGCTCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTGGACTTCCGc                                                                                                                                              <  2:1078759/122‑1 (MQ=255)
                                     tGCTCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTgcg                                                                                                                 >  1:2170251/1‑150 (MQ=255)
                                        tCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCa                                                                                                                                                     >  2:1449068/1‑111 (MQ=255)
                                        tCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCa                                                                                                                                                     <  1:1449068/111‑1 (MQ=255)
                                        tCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCAc                                                                                                                                                    <  2:1055715/112‑1 (MQ=255)
                                        tCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCAc                                                                                                                                                    >  1:1055715/1‑112 (MQ=255)
                                        tCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCgg                                                                                                              >  2:304471/1‑150 (MQ=255)
                                          gcgGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCt                                                                                                            <  2:2188974/150‑1 (MQ=255)
                                           cgGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCAGGTGTTCTACACCTCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTg                                                                                                           <  2:1298101/150‑1 (MQ=255)
                                            gggCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCACTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTgg                                                                                                          <  2:1087257/150‑1 (MQ=255)
                                              gCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGt                                                                                                        <  1:1097902/150‑1 (MQ=255)
                                              gCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGt                                                                                                        <  1:2349435/150‑1 (MQ=255)
                                               cAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGtt                                                                                                       <  2:80758/150‑1 (MQ=255)
                                                aaaaCGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTc                                                                                                      >  2:1795319/1‑150 (MQ=255)
                                                  aaCGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTcgc                                                                                                    >  1:840190/1‑150 (MQ=255)
                                                        cAACCTGTGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACCCACTACCGTACGATCTTCTGCGCGGCGGGGTTCGCCGCCt                                                                                               >  2:21992/1‑150 (MQ=255)
                                                         aaCCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTCACGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTg                                                                                              <  1:1636833/150‑1 (MQ=255)
                                                             tGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCa                                                                                          <  1:1874784/150‑1 (MQ=255)
                                                                  gccgcACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACCTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACtgtg                                                                                     >  1:1124101/1‑150 (MQ=255)
                                                                     gcACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACtgtg                                                                                     <  2:2344947/147‑1 (MQ=255)
                                                                     gcACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACtgtg                                                                                     >  1:2344947/1‑147 (MQ=255)
                                                                     gcACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTgtggtg                                                                                  <  2:532986/150‑1 (MQ=255)
                                                                       aCTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCt                                                                                >  2:2392069/1‑150 (MQ=255)
                                                                        ctctGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTg                                                                               >  1:2000928/1‑150 (MQ=255)
                                                                         tctGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTgg                                                                              <  1:1763739/150‑1 (MQ=255)
                                                                               aGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCACGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTgg                                                                        >  2:346116/1‑150 (MQ=255)
                                                                                 cTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTc                                                                      <  1:2430310/150‑1 (MQ=255)
                                                                                  tGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCa                                                                            <  2:2114442/143‑1 (MQ=255)
                                                                                  tGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCa                                                                            >  1:2114442/1‑143 (MQ=255)
                                                                                    gCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGc                                                                   >  1:1314128/1‑150 (MQ=255)
                                                                                    gCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGc                                                                   >  1:765619/1‑150 (MQ=255)
                                                                                    gCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGc                                                                   <  2:2120358/150‑1 (MQ=255)
                                                                                      ccAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTGGACTTCCGcca                                                                                                                                            <  1:709853/75‑3 (MQ=255)
                                                                                      ccAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTGGACTTCCGcca                                                                                                                                            >  2:709853/1‑73 (MQ=255)
                                                                                      ccAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTcgc                                                                                                    <  2:2085822/115‑1 (MQ=255)
                                                                                      ccAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTcgc                                                                                                    >  1:2085822/1‑115 (MQ=255)
                                                                                         ggCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCAc                                                                           >  2:2127466/1‑137 (MQ=255)
                                                                                         ggCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCAc                                                                           <  1:2127466/137‑1 (MQ=255)
                                                                                                  tggTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTggtggcggtgg                                                      >  2:1173090/1‑150 (MQ=255)
                                                                                                   ggTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTggtggcggtggc                                                     >  2:1745197/1‑150 (MQ=255)
                                                                                                      gtTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGa                                                  >  1:1098074/1‑150 (MQ=255)
                                                                                                       ttCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTg                                                                               >  2:728654/1‑120 (MQ=255)
                                                                                                       ttCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAg                                                 >  1:568136/1‑150 (MQ=255)
                                                                                                       ttCNACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAg                                                 >  2:1503163/1‑150 (MQ=255)
                                                                                                          tACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCAc                                              >  2:1311289/1‑150 (MQ=255)
                                                                                                            cacCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCa                                            >  1:263604/1‑150 (MQ=255)
                                                                                                            cacCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCa                                            >  1:1586628/1‑150 (MQ=255)
                                                                                                             acCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAg                                           >  2:2271429/1‑150 (MQ=255)
                                                                                                               cGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAGcc                                         >  1:1439744/1‑150 (MQ=255)
                                                                                                                gCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAGCCt                                        >  1:1228182/1‑150 (MQ=255)
                                                                                                                    aCCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAGCCTTCgg                                    <  2:589945/150‑1 (MQ=255)
                                                                                                                    aCCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAGCCTTCgg                                    <  2:2000928/150‑1 (MQ=255)
                                                                                                                    aCCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAGCCTTCgg                                    <  2:1439744/150‑1 (MQ=255)
                                                                                                                     ccGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCAc                                                                           >  1:2405364/1‑110 (MQ=255)
                                                                                                                     ccGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTggtggcg                                                          >  1:1035306/1‑127 (MQ=255)
                                                                                                                     ccGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGa                                                  <  1:1818200/135‑1 (MQ=255)
                                                                                                                     ccGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACCTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGa                                                  >  2:1818200/1‑135 (MQ=255)
                                                                                                                        tATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAGCCTTCGGCa                                  >  1:1500670/1‑148 (MQ=255)
                                                                                                                         aTCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAGCCTTCGGCACTg                               <  1:2006786/150‑1 (MQ=255)
                                                                                                                            cTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAGCCTTc                                      >  2:2397887/1‑140 (MQ=255)
                                                                                                                            cTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAGCCTTc                                      <  1:2397887/140‑1 (MQ=255)
                                                                                                                             tCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAGCCTTCGGCACTGGCAg                           >  1:1051876/1‑150 (MQ=255)
                                                                                                                             tCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAGCCTTCGGCACTGGCAg                           >  2:1665354/1‑150 (MQ=255)
                                                                                                                             tCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAGCCTTCGGCACTGGCAg                           >  2:1572914/1‑150 (MQ=255)
                                                                                                                               gCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAGCCTTCGGCACTGGCAgac                         <  1:681742/150‑1 (MQ=255)
                                                                                                                               gCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAGCCTTCGGCACTGGCAgac                         <  1:1903841/150‑1 (MQ=255)
                                                                                                                                   ccGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAGCCTTCGGCACTGGCAGacgacg                     >  2:1266371/1‑150 (MQ=255)
                                                                                                                                        ggcgtggcgGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGCTTCCCCTCCGGGGCCCTATGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAGCCTTCGGCACTGGCAGACGACGTGTTg                <  1:766217/150‑1 (MQ=255)
                                                                                                                                          gcgtggcgGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAGCCTTCGGCACTGGCAGACGACGTGTTgg               <  1:528543/150‑1 (MQ=255)
                                                                                                                                            gtggcgGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAGCCTTCGGCACTGGCAGACGACGTGTTGGGc             >  1:1963964/1‑150 (MQ=255)
                                                                                                                                                  gAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAGCCTTCGGCACTGGCAGACGACGTGTTGGGCATCGGc       <  1:1511566/150‑1 (MQ=255)
                                                                                                                                                     gTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAGCCTTCGGCACTgg                              >  2:2350825/1‑124 (MQ=255)
                                                                                                                                                     gTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAGCCTTCGGCACTGGCAGACGACGTGTTGGGCATCGGCgtg    >  1:1444559/1‑150 (MQ=255)
                                                                                                                                                       ccACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAGCCTTCg                                     <  2:2390003/115‑1 (MQ=255)
                                                                                                                                                       ccACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAGCCTTCg                                     >  1:2390003/1‑115 (MQ=255)
                                                                                                                                                       ccACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAGCCTTCGGCACTGGCAGACGACGTGTTGGGCATCGGCgtggt  >  1:345429/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
TTGGCCTGGCGCTGTGCAGCTATGCCTTGCAGATCGATGCTCGCGGGCAAAACGCCAACCTGGGCGCCGCACTCTGGCAGCTGGCCCAGGCCTGGCTGGTGTTCTACACCGCCTACCGTATCCTCGCACCCGGTGGCGTGGCGGAAGTCCACTTCCGCTGGCACAAACCGCAGGTCGAGTTCCTGCGCGGCTGGGTTCGCCGCCTGGGCACTGTGGTGCTGGCACTGGTCGGCGTGGTGGCGGTGGCCGAGCACCAGCCTTCGGCACTGGCAGACGACGTGTTGGGCATCGGCGTGGT  >  NC_002947/5783227‑5783524

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: