Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,336,227 G→A 6.3% R870H (CGC→CAC)  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,336,2270GA 6.3% 320.0 / 8.0 110R870H (CGC→CAC) PP_4695two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (53/50);  new base A (2/5);  total (55/55)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

Only 100 of 112 total aligned reads displayed.

GGCAGCCACCAGAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAACATCA  >  NC_002947/5336079‑5336373
                                                                                                                                                    |                                                                                                                                                   
ggcAGCCACCAGAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAgcgc                                                                                                                                                    <  1:465165/150‑1 (MQ=255)
     ccaccaGAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGc                                                                                                                                               <  2:1426161/150‑1 (MQ=255)
       accaGAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGccc                                                                                                                                             >  1:1951669/1‑150 (MQ=255)
       accaGAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGccc                                                                                                                                             >  1:1562716/1‑150 (MQ=255)
        ccaGAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCACCTGGGCCa                                                                                                                                            >  2:1396517/1‑150 (MQ=255)
         caGAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCACTTTCGACGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGGGACCCGCAGCGCCTGGCCCAg                                                                                                                                           >  1:2382850/1‑150 (MQ=255)
         caGAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAg                                                                                                                                           <  1:208040/150‑1 (MQ=255)
          aGAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGCTGCCATCGGTCTGCTGGCGTTGAACCGGCGTAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGCACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAgg                                                                                                                                          <  2:2085113/150‑1 (MQ=255)
             aCAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGc                                                                                                                                       <  2:844519/150‑1 (MQ=255)
               aGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTg                                                                                                                                     <  2:2442364/150‑1 (MQ=255)
                gCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGa                                                                                                                                    <  1:1563839/150‑1 (MQ=255)
                gCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTGCGAAGCACAGTTCTTCAACCTTTGCGCCCCGGACCATTGGGCCGCAGGTGACCCGCAGCACCTGGGCCAGGCGCTgc                                                                                                                                    <  2:1331138/150‑2 (MQ=255)
                 cGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGAt                                                                                                                                   >  2:1171408/1‑150 (MQ=255)
                  gaggagCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATc                                                                                                                                  <  2:1600459/150‑1 (MQ=255)
                  gaggagCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCACCTGGGCCAGGCGCTgcgg                                                                                                                                  <  1:808823/150‑4 (MQ=255)
                  gaggagCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCACCTGGGCCAGGCGCTgcgg                                                                                                                                  <  1:1396517/150‑4 (MQ=255)
                    ggagCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGCCGGACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCACCTGGGCCAGGCGCTgcgggt                                                                                                                                <  1:156960/150‑6 (MQ=255)
                       gCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCt                                                                                                                             >  2:2086510/1‑150 (MQ=255)
                       gCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCAACTGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCt                                                                                                                             >  1:2017565/1‑150 (MQ=255)
                                 cTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCTCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAAc                                                                                                                   >  2:2171476/1‑150 (MQ=255)
                                   ggCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGc                                                                                                                 <  2:1031829/150‑1 (MQ=255)
                                     cTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGccc                                                                                                               >  2:2013326/1‑150 (MQ=255)
                                      tGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCcg                                                                                                              >  2:2322769/1‑150 (MQ=255)
                                         aGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACCGTTCCTCACCCTCTGCTATTCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGa                                                                                                           <  1:2120811/150‑1 (MQ=255)
                                         aGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGa                                                                                                           <  2:1131100/150‑1 (MQ=255)
                                          gTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGAt                                                                                                          <  2:2032871/150‑1 (MQ=255)
                                           tGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATg                                                                                                         >  2:250033/1‑150 (MQ=255)
                                               acaAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTc                                                                                                     >  1:727142/1‑150 (MQ=255)
                                               acaAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTc                                                                                                     >  2:349960/1‑150 (MQ=255)
                                                 aaGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGc                                                                                                        <  1:888676/145‑1 (MQ=255)
                                                 aaGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGc                                                                                                        >  2:888676/1‑145 (MQ=255)
                                                  aGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTcgcc                                                                                                  <  1:452227/150‑1 (MQ=255)
                                                   gATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAAGTGCGTGTGAGCCTGCTCTCCAACGCCCGCGATGCCTcgccg                                                                                                 >  2:1773168/1‑150 (MQ=255)
                                                      gCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGccc                                                                                              <  2:1672201/150‑1 (MQ=255)
                                                      gCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGccc                                                                                              <  2:1951669/150‑1 (MQ=255)
                                                        cATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCgg                                                                                            <  1:2322769/150‑1 (MQ=255)
                                                        cATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCgg                                                                                            >  2:104557/1‑150 (MQ=255)
                                                         aTCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGc                                                                                           >  1:223318/1‑150 (MQ=255)
                                                         aTCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGc                                                                                           >  1:454260/1‑150 (MQ=255)
                                                         aTCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGc                                                                                           >  1:1248189/1‑150 (MQ=255)
                                                         aTCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCACCTGGGCCAGGCGCTgcgg                                                                                                                                  <  2:1630616/111‑4 (MQ=255)
                                                         aTCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCACCTGGGCCAGGCGCTgcgg                                                                                                                                  >  1:1630616/1‑108 (MQ=255)
                                                          tCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCa                                                                                          >  2:1394151/1‑150 (MQ=255)
                                                             gTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGTGATGCCTCGCCGCCCGGCAgcg                                                                                       >  2:977747/1‑150 (MQ=255)
                                                             gTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAgcg                                                                                       <  1:1089702/150‑1 (MQ=255)
                                                              tCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGAAgcgc                                                                                      >  1:2051214/1‑150 (MQ=255)
                                                                    ggCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGACAGCGCcgtgcg                                                                                >  2:579818/1‑150 (MQ=255)
                                                                         tGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACGTTTGCGACCCGGACCATTGGGCCGAAGCTCAGCTTCTTCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCcgtgcgcgtgc                                                                           <  1:1223111/150‑1 (MQ=255)
                                                                            aCCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCcgtgc                                                                                 <  1:1385059/141‑1 (MQ=255)
                                                                            aCCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCcgtgc                                                                                 >  2:1385059/1‑141 (MQ=255)
                                                                            aCCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCa                                                                       <  1:570265/150‑1 (MQ=255)
                                                                            aCCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCa                                                                       <  1:2086510/150‑1 (MQ=255)
                                                                             ccGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAg                                                                      >  1:2239467/1‑150 (MQ=255)
                                                                               ggCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCACGCGCacgg                                                                    <  2:223318/150‑4 (MQ=255)
                                                                                gcgcAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGa                                                                   <  1:218701/150‑1 (MQ=255)
                                                                                 cgcAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAg                                                                  <  1:515937/150‑1 (MQ=255)
                                                                                  gcAATTTCGCAGTACAGTTCTTCCACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCACCGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAgg                                                                 <  2:1099214/150‑1 (MQ=255)
                                                                                  gcAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAgg                                                                 <  1:732409/150‑1 (MQ=255)
                                                                                   cAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGt                                                                <  1:1171408/150‑1 (MQ=255)
                                                                                     aTTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTa                                                              <  1:903352/150‑1 (MQ=255)
                                                                                      tttCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTaa                                                             >  1:1555349/1‑150 (MQ=255)
                                                                                       ttCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAAc                                                            <  1:104557/150‑1 (MQ=255)
                                                                                       ttCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAAc                                                            <  1:1394151/150‑1 (MQ=255)
                                                                                          gAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAg                                                         >  1:1386667/1‑150 (MQ=255)
                                                                                           aaGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGc                                                        <  2:2382850/150‑1 (MQ=255)
                                                                                           aaGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGc                                                        >  2:1597239/1‑150 (MQ=255)
                                                                                           aaGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGc                                                        >  2:2175985/1‑150 (MQ=255)
                                                                                            aGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGca                                                       >  1:2354517/1‑150 (MQ=255)
                                                                                                cAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACAcc                                                   >  2:172702/1‑150 (MQ=255)
                                                                                                  gTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGt                                                 >  2:411556/1‑150 (MQ=255)
                                                                                                   ttcttcAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCTCGCGATGCCTCGCCGCGCGGCAGCGCCGTTCGCGTGCGCAGCGAGGTTAACGAGCACACCGTc                                                <  1:47340/150‑1 (MQ=255)
                                                                                                   ttcttcAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCTTCGCCCGCCAGCTCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTc                                                <  1:47143/150‑1 (MQ=255)
                                                                                                      ttcAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGAc                                             <  1:2175985/150‑1 (MQ=255)
                                                                                                      ttcAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCGCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGGTGGCGCAGTGAGGTTAACGAGCACACCGTCGAc                                             <  1:1773168/150‑1 (MQ=255)
                                                                                                       tcAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCcgtgcgcgtgc                                                                           >  1:1038287/1‑120 (MQ=255)
                                                                                                       tcAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCcgtgcgcgtgc                                                                           <  2:1038287/120‑1 (MQ=255)
                                                                                                       tcAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGAcc                                            >  2:80688/1‑150 (MQ=255)
                                                                                                        cAANCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCt                                           >  2:1500643/1‑150 (MQ=255)
                                                                                                         aaCCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTg                                          >  1:498243/1‑150 (MQ=255)
                                                                                                          aCCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGa                                         <  2:2094358/150‑1 (MQ=255)
                                                                                                           ccTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGAt                                        >  2:68176/1‑150 (MQ=255)
                                                                                                            cTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGAtc                                       >  1:1721860/1‑150 (MQ=255)
                                                                                                                gCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGAtcgtcg                                   <  1:2013326/150‑1 (MQ=255)
                                                                                                                  gACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGaa                                 <  2:1846325/150‑1 (MQ=255)
                                                                                                                    cccGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGACAGATTCCGAACGAGCACACCGTCGACCTGATCGTCGAAGa                               >  1:723352/1‑150 (MQ=255)
                                                                                                                    cccGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGa                               >  2:716312/1‑150 (MQ=255)
                                                                                                                      cGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACg                             >  2:2140917/1‑150 (MQ=255)
                                                                                                                              tGGGCCGAAGGTGACCCGCAGCGCCTGGCACAAATACTGTTCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGc                     <  2:2017565/150‑1 (MQ=255)
                                                                                                                              cgggCCGCAGGTGCGCTCCACCGACTGGCCGAGTCGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGc                     <  2:1721860/149‑1 (MQ=255)
                                                                                                                               tggCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCATCACCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCg                    >  2:344986/2‑150 (MQ=255)
                                                                                                                                ggCCTAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCgg                   <  1:1290318/150‑1 (MQ=255)
                                                                                                                                  ccGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGa                 >  1:359893/1‑150 (MQ=255)
                                                                                                                                   cGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAGCGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGCCAAAGGCAGCGGGaa                >  2:1321535/1‑149 (MQ=255)
                                                                                                                                    gAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGAtt               <  1:820029/150‑1 (MQ=255)
                                                                                                                                    gAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGAtt               >  1:696646/1‑150 (MQ=255)
                                                                                                                                         tGACCCGCAGCGCCTGGCCCAGGTGCTGCCCATCCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCgaa          <  2:2136014/150‑1 (MQ=255)
                                                                                                                                             ccGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTTCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAAc      >  2:763453/1‑150 (MQ=255)
                                                                                                                                             ccGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGAAGCGTAGCGGGCTTGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAAc      >  2:920408/1‑150 (MQ=255)
                                                                                                                                                 ggcgcCGTGCCCAGGTGCGGATCAACCTGCTCTCCAACGCCCGCTATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAAcatca  <  2:543608/149‑1 (MQ=255)
                                                                                                                                                 aGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAAcatca  <  1:2325618/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGCAGCCACCAGAACAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAACATCA  >  NC_002947/5336079‑5336373

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: