Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,861,202 C→A 5.2% R27R (CGC→CGA PP_3413 → two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,861,2020CA 5.2% 246.2 / 2.1 77R27R (CGC→CGAPP_3413two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base C (41/32);  new base A (2/2);  total (43/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

AGACCGAGTCACTGGTGGAAATGGTCGACCTGGCCAAGCGCAACAGCGTCTATTGAGCGCTCGGCTGCCATGAAACAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGTGCGCGCCCCGATTACCCCCCGTTCGAAATCAACGTCAGCCAGCATGA  >  NC_002947/3861053‑3861348
                                                                                                                                                      |                                                                                                                                                  
aGACCGAGTCACTGGTGGAAATGGTCGACCTGGCCAAGCGCAACAGCGTCTATTGAGCGCTCGGCTGCCATGAAACAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGc                                                                                                                                                    <  1:2024694/150‑1 (MQ=255)
      aGTCACTGGTGGAAATGGTCGACCTGGCCAAGCGCAACAGCGTCTATTGAGCGCTCGGCTGCCATGAAACAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATctg                                                                                                                                              >  1:851657/1‑150 (MQ=255)
      aGTCACTGGAGGAAATGGTCGACCTGGCCAAGCGCAACAGCGTCTATTGAGCGCTCGGCTGCCATGAAACAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATctg                                                                                                                                              >  1:2069988/1‑150 (MQ=255)
         cACTGGTGGAAATGGTCGACCTGGCCAAGCGCAACAGCGTCTATTGAGCGCTCGGCTGCCATGAAACAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATctgctg                                                                                                                                           <  1:1409818/150‑1 (MQ=255)
          aCTGGTGGAAATGGTCGACCTGGCCAAGCGCAACAGCGTCTATTGAGCGCTCGGCTGCCATGAAACAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTgg                                                                                                                                          >  1:1957168/1‑150 (MQ=255)
           cTGGTGGAAATGGTCGACCTGGCCAAGCGCAACAGCGTCTATTGAGCGCTCGGCTGCCATGAAACAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGgc                                                                                                                                         <  1:702289/150‑1 (MQ=255)
               tggAAATGGTCGACCTGGCCAAGCGCAACAGCGTCTATTGAGCGCTCGGCTGCCATGAAACAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAAGTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGgcgcgc                                                                                                                                     >  2:1564241/1‑150 (MQ=255)
                    aTGGTCGACCTGGCCAAGCGCAACAGCGTCTATTGAGCGCTCGGCTGCCATGAAACAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGAt                                                                                                                                >  1:1360889/1‑150 (MQ=255)
                     tGGTCGACCTGGCCAAGCGCAACAGCGTCTATTGAGCGCTCGGCTGCCATGAAACAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATc                                                                                                                               <  2:2069988/150‑1 (MQ=255)
                       gTCGACCTGGCCAAGCGCAACAGCGTCTATTGAGCGCTCGGCTGCCATGAAACAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAg                                                                                                                             >  2:928986/1‑150 (MQ=255)
                       gTCGACCTGGCCAAGCGCAACAGCGTCTATTGAGCGCTCGGCTGCCATGAAACAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAg                                                                                                                             >  1:487236/1‑150 (MQ=255)
                                gCCAAGCGCAACAGCGTCTATTGAGCGCTCGGCTGCCATGAAACAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGAGCTCTATCTGTGCCGCTGc                                                                                                                    <  2:688906/150‑1 (MQ=255)
                                gCCAAGCGCAACAGCGTCTATTGAGCGCTCGGCTGCCATGAAACAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGAAGGCAAAACGTGAACCACGCCAACGTCTGGCGCGCGCGTGCAGTGCCGCTGc                                                                                                                    >  2:1777101/1‑150 (MQ=255)
                                               gTCTATTGAGCGCTCGGCTGCCATGAAACAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATc                                                                                                     >  1:1268188/1‑150 (MQ=255)
                                                  tATTGACCGCTCTGCTACCATGAAAGAGTATCAGGCGGGCGTCCTTCTGCCTATCGGTCTTGCCAAAAACACGGAACTGCAGGCCAGCAGTAAACCACACCATCTGTTGGCGCGGTCGAACAGTGGCGCTGCCCGACTGAAGTTATccag                                                                                                  >  2:517168/1‑148 (MQ=255)
                                                     tGAGCGCTCGGCTGCCATGAAACAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCga                                                                                               <  2:1268188/150‑1 (MQ=255)
                                                     tGAGCGCTCGGCTGCCATGAAACAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCga                                                                                               <  2:2290880/150‑1 (MQ=255)
                                                          gcTCGGCTGCCATGAAACAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGc                                                                                                                          >  2:301216/1‑118 (MQ=255)
                                                          gcTCGGCTGCCATGAAACAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGc                                                                                                                          <  1:301216/118‑1 (MQ=255)
                                                          gcTCGGCTGCCATGAAACAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATc                                                                                          >  2:1515902/1‑150 (MQ=255)
                                                            tCGGCTGCCATGAAACAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGa                                                                                                                                 >  2:173690/1‑109 (MQ=255)
                                                            tCGGCTGCCATGAAACAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGa                                                                                                                                 <  1:173690/109‑1 (MQ=255)
                                                            tCGGCTGCCATGAAACAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCTCTGACGTTATCCGGCGATGATCGc                                                                                        <  1:1594073/150‑1 (MQ=255)
                                                                 tGCCATGAAACAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTTACGTTATCCGGCGATGATCGCCACTg                                                                                   >  2:821229/1‑150 (MQ=255)
                                                                    cATGAAACAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCt                                                                                <  1:1564241/150‑1 (MQ=255)
                                                                        aaaCAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCCCGATCAGTGCCGCTGCCCCACTGCCGTTATCCGGCGATGATCGCCACTGGCTGaaa                                                                            >  2:659513/1‑150 (MQ=255)
                                                                            aGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAaacaac                                                                        >  1:1426927/1‑150 (MQ=255)
                                                                              ttGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCTACTGGCTGAAACAAcgc                                                                      <  1:1940400/150‑1 (MQ=255)
                                                                                  tGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCCCGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAAc                                                                  >  1:1980206/1‑150 (MQ=255)
                                                                                         cGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCAATGGCTGAAGCAACGCGAACGCCTGGt                                                           >  1:273746/1‑150 (MQ=255)
                                                                                           tCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTaa                                                         >  2:1420685/1‑150 (MQ=255)
                                                                                             cTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATg                                                       >  1:455652/1‑150 (MQ=255)
                                                                                              ttCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATgg                                                      >  1:1662189/1‑150 (MQ=255)
                                                                                              ttCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATgg                                                      >  1:188972/1‑150 (MQ=255)
                                                                                               tCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTTACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATggg                                                     <  1:821229/150‑1 (MQ=255)
                                                                                               tCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAAAAACGCGAAAGCCTGATAATggg                                                     <  2:1744337/150‑1 (MQ=255)
                                                                                                 tGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGCCAGCGAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCa                                                   <  1:1777101/150‑1 (MQ=255)
                                                                                                  ggCTATCGGTCTGGCCAACTCTACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAg                                                  <  2:538360/150‑1 (MQ=255)
                                                                                                   gCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGc                                                                                        >  2:2060748/1‑111 (MQ=255)
                                                                                                   gCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGc                                                                                        <  1:2060748/111‑1 (MQ=255)
                                                                                                   gCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGACATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGt                                                 >  2:587836/1‑150 (MQ=255)
                                                                                                   gCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGACATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGt                                                 >  2:586502/1‑150 (MQ=255)
                                                                                                     tATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGTgc                                               >  1:1883525/1‑150 (MQ=255)
                                                                                                            cTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGTGCGCGcccc                                        <  1:1515902/150‑1 (MQ=255)
                                                                                                             tGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATgg                                                      >  1:335395/1‑135 (MQ=255)
                                                                                                             tGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATgg                                                      <  2:335395/135‑1 (MQ=255)
                                                                                                             tGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGt                                                 <  1:922155/140‑1 (MQ=255)
                                                                                                             tGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGt                                                 >  2:922155/1‑140 (MQ=255)
                                                                                                                ccAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATggg                                                     >  2:666741/1‑133 (MQ=255)
                                                                                                                ccAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATggg                                                     <  1:666741/133‑1 (MQ=255)
                                                                                                                     tCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGTGCGCGCCCCGATTAcccc                               <  2:1360889/150‑1 (MQ=255)
                                                                                                                     tCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGTGCGCGCCCCGATTAcccc                               >  1:1428682/1‑150 (MQ=255)
                                                                                                                       cACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGTGCGCGCCCCGATTAcccccc                             >  2:1473589/1‑150 (MQ=255)
                                                                                                                       cACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGTGCGCGCCCCGATTAcccccc                             >  1:1948020/1‑150 (MQ=255)
                                                                                                                          ggAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCCATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGTGCGCGCCCCGATTACCCCCCGtt                          <  1:646184/150‑1 (MQ=255)
                                                                                                                           gAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGTGCGCGCCCCGATTACCCCCCGTTc                         <  1:68746/150‑1 (MQ=255)
                                                                                                                           gAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGTGCGCGCCCCGATTACCCCCCGTTc                         <  2:455652/150‑1 (MQ=255)
                                                                                                                               tGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACg                                                                 <  2:1554261/106‑1 (MQ=255)
                                                                                                                               tGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACg                                                                 >  1:1554261/1‑106 (MQ=255)
                                                                                                                               tGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGTGCGCGCCCCGATTACCCCCCGTTCGaaa                     >  2:1033149/1‑150 (MQ=255)
                                                                                                                               tGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGCCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGTGCGCGCCCCGATTACCCCCCGTTCGaaa                     <  2:851657/150‑1 (MQ=255)
                                                                                                                                gCAGGCCAGCAGTGAACCACGACATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGTGCGCGCCCCGATTACCCCCCGTTCGAAAt                    <  1:587836/150‑1 (MQ=255)
                                                                                                                                gCAGGCCAGCAGTGAACCACGACATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGTGCGCGCCCCGATTACCCCCCGTTCGAAAt                    <  1:586502/150‑1 (MQ=255)
                                                                                                                                    gCCAGCAGAGAACCACGCCATCTGCTGGCGAGCTCGCACCGTCACCCAGTCGCTGTGACGTTATCCGGCGATGAGCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGTGCGCGCCCCGATTACCCCCCGTTCGAAATCAAc                >  1:45191/1‑150 (MQ=255)
                                                                                                                                      cagcagTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTgg                                                            <  2:2046502/104‑1 (MQ=255)
                                                                                                                                      cagcagTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTgg                                                            >  1:2046502/1‑104 (MQ=255)
                                                                                                                                      cagcagTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTg                                                             <  1:1326200/103‑1 (MQ=255)
                                                                                                                                      cagcagTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTg                                                             >  2:1326200/1‑103 (MQ=255)
                                                                                                                                      cagcagTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGTGCGCGCCCCGATTACCCCCCGTTCGAAATCAACGt              <  1:1473589/150‑1 (MQ=255)
                                                                                                                                           agTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGTGCGCGCCCCGATTACCCCCCGTTCGAAATCAACGTcagc          <  1:928986/150‑1 (MQ=255)
                                                                                                                                             tGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGTGCGCGCCCCGATCACCCCCCGTTCGAAATCAACGTcagcca        >  2:192988/1‑150 (MQ=255)
                                                                                                                                               aaCCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGTGCGCGCCCCGATTACCCCCCGTTCGAAATCAACGTcagccagc      >  1:2294410/1‑150 (MQ=255)
                                                                                                                                                aCCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGTGCGCGCCCCGAt                                     >  1:930496/1‑118 (MQ=255)
                                                                                                                                                aCCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGTGCGCGCCCCGAt                                     <  2:930496/118‑1 (MQ=255)
                                                                                                                                                 ccACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGTGCGCGCCCCGATTACCCCCCGTTCGAAATCAACGTCAGCCAGCAt    >  1:2127602/1‑150 (MQ=255)
                                                                                                                                                   aCGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCTACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGTGCGCGCCCCGATTACCCCCCGTTCGAAATCAACGTCAGCCAGCATGa  >  2:1348701/1‑150 (MQ=255)
                                                                                                                                                   aCGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGTGCGCGCCCCGATTACCCCCCGTTCGAAATCAACGTCAGCCAGCATGa  <  2:1426927/150‑1 (MQ=255)
                                                                                                                                                   aCGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGTGCGCGCCCCGATTACCCCCCGTTCGAAATCAACGTCAGCCAGCATGa  >  1:242334/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
AGACCGAGTCACTGGTGGAAATGGTCGACCTGGCCAAGCGCAACAGCGTCTATTGAGCGCTCGGCTGCCATGAAACAGTTGCTGGCGAGCGTCCTTCTGGCTATCGGTCTGGCCAACTCCACGGAACTGCAGGCCAGCAGTGAACCACGCCATCTGCTGGCGCGCTCGATCAGTGCCGCTGCCCCACTGACGTTATCCGGCGATGATCGCCACTGGCTGAAACAACGCGAACGCCTGGTAATGGGCAGTGCGCGCCCCGATTACCCCCCGTTCGAAATCAACGTCAGCCAGCATGA  >  NC_002947/3861053‑3861348

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