Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,119,244 G→C 6.3% E324Q (GAG→CAG)  PP_4505 → sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,119,2440GC 6.3% 244.5 / 5.3 79E324Q (GAG→CAG) PP_4505sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (42/32);  new base C (2/3);  total (44/35)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

CCTGGAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTGATCAGC  >  NC_002947/5119099‑5119393
                                                                                                                                                 |                                                                                                                                                     
ccTGGAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGgc                                                                                                                                                   <  2:1100390/150‑1 (MQ=255)
 cTGGAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGgcg                                                                                                                                                  >  2:1400430/1‑150 (MQ=255)
   ggAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCa                                                                                                                                                >  2:324937/1‑150 (MQ=255)
    gAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAg                                                                                                                                               >  2:171519/1‑150 (MQ=255)
       gACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGCGGATGCCGAACAAGCCCCTGCCGAACCGGTCGACCTCAATGCCCTTCTCGGCAGCGTGCGCAACGACGCCCCGCCGAGCGAGCCGGATCCGGACCTGCGCCTCCAGGCGCacgcc                                                                                                                                            <  2:798901/150‑5 (MQ=255)
        aCCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAg                                                                                                                                           <  2:1193071/150‑1 (MQ=255)
         ccTGATCAGCGAAATCCTCCCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTCCAGGCGCacgtcct                                                                                                                                          <  2:1253472/150‑7 (MQ=255)
          cTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGt                                                                                                                                         >  2:841104/1‑150 (MQ=255)
          cTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGt                                                                                                                                         >  1:700020/1‑150 (MQ=255)
            gATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGttt                                                                                                                                       <  1:132676/150‑1 (MQ=255)
                  cGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACCAGCCTAACCAGCACCGGTCGACCTCAATGCCCTATTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCtt                                                                                                                                 >  1:388216/1‑150 (MQ=255)
                   gAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTg                                                                                                                                <  2:773812/150‑1 (MQ=255)
                    aaaTCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGc                                                                                                                               >  1:992408/1‑150 (MQ=255)
                       tCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAgg                                                                                                                            >  1:767808/1‑150 (MQ=255)
                       tCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAgg                                                                                                                            <  1:841104/150‑1 (MQ=255)
                         cTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCTGGATCAGGACGTGCGCCTCCAGGCGCacgtcc                                                                                                                                           >  1:924714/1‑128 (MQ=255)
                         cTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCTGGATCAGGACGTGCGCCTCCAGGCGCacgtcc                                                                                                                                           <  2:924714/133‑6 (MQ=255)
                         cTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGt                                                                                                                          <  2:700020/150‑1 (MQ=255)
                           cGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTg                                                                                                                        <  1:171519/150‑1 (MQ=255)
                            gCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTgg                                                                                                                       >  1:553339/1‑150 (MQ=255)
                              ccTGGCCCGGGTGGATGCCGAACAAGCCCATACCGAACCGGTCTACCTCAATGCCCTACTCGGCAGCGTGCGCAAGTACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTgggc                                                                                                                     >  2:195897/1‑148 (MQ=255)
                                 ggCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAAc                                                                                                                  <  1:1400430/150‑1 (MQ=255)
                                  gCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTTGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTACTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACg                                                                                                                 >  1:660790/1‑150 (MQ=255)
                                    ccGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTAAATGCACTACTAGGCAAGGTGGGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGGCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAAccat                                                                                                               >  2:194504/1‑147 (MQ=255)
                                     cGGGTGGATGCCGCACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCa                                                                                                                                                >  2:1059751/1‑116 (MQ=255)
                                     cGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCa                                                                                                                                                <  1:1059751/116‑1 (MQ=255)
                                     cGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGTGTTGAGCTTGCAGGGTTGGCCAACGCTg                                                                                                              <  1:750047/150‑1 (MQ=255)
                                          ggATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGa                                                                                                         <  2:992408/150‑1 (MQ=255)
                                                cGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGc                                                                                                   >  1:1358721/1‑150 (MQ=255)
                                                   aCAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACGCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGttt                                                                                                                                       <  1:362255/111‑1 (MQ=255)
                                                   aCAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACACGGCAGCGTGCGCAAGGGCGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCGGGAGGCGCAGCCAGGttt                                                                                                                                       >  2:362255/1‑111 (MQ=255)
                                                         ccATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACaa                                                                                          <  2:1358902/150‑1 (MQ=255)
                                                          cATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAAc                                                                                         <  1:913398/150‑1 (MQ=255)
                                                                gAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGa                                                                                                             >  2:995984/1‑124 (MQ=255)
                                                                gAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGa                                                                                                             <  1:995984/124‑1 (MQ=255)
                                                                gAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACctgctg                                                                                   >  2:622740/1‑150 (MQ=255)
                                                                 aaCCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCtgctgc                                                                                  <  2:1358721/150‑1 (MQ=255)
                                                                           ccTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACctgct                                                                                    <  2:557340/138‑1 (MQ=255)
                                                                           ccTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACctgct                                                                                    >  1:557340/1‑138 (MQ=255)
                                                                           ccTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCt                                                                        <  2:767808/150‑1 (MQ=255)
                                                                           ccTCAATGCCCTACTCGCCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCt                                                                        >  2:730537/1‑150 (MQ=255)
                                                                               aaTGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACCCCctgct                                                                                    <  2:1065246/134‑1 (MQ=255)
                                                                               aaTGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACctgct                                                                                    >  1:1065246/1‑134 (MQ=255)
                                                                               aaTGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATg                                                                            <  1:242811/142‑1 (MQ=255)
                                                                               aaTGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATg                                                                            >  2:242811/1‑142 (MQ=255)
                                                                                     cTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAAc                                                              >  1:1350561/1‑150 (MQ=255)
                                                                                      tACTCGGCAGCGTGCGCAAGGACGCCCTGCTTAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAAcc                                                             >  1:617019/1‑150 (MQ=255)
                                                                                      tACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAAcc                                                             >  2:540870/1‑150 (MQ=255)
                                                                                      tACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAAcc                                                             >  1:1201522/1‑150 (MQ=255)
                                                                                      tACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAAcc                                                             >  1:160376/1‑150 (MQ=255)
                                                                                         tCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCcgg                                                          >  1:567471/1‑150 (MQ=255)
                                                                                           ggCAGCGTCCGCAAGGACGCCCTGCTGAGCGCGCCGTATCAGGACGTGCGCCTGGAGGCGCAGCCAGCTTTGAGCTTGCAGGGTTGGCGAACGCTGATCGAGCGTGCCGTTGACAATCTGCTGCGCAATGCCCTGCGCTTCAACCcggcg                                                        <  2:494374/150‑1 (MQ=255)
                                                                                            gCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACg                                                                                                                 >  1:659478/1‑92 (MQ=255)
                                                                                            gCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACg                                                                                                                 <  2:659478/92‑1 (MQ=255)
                                                                                                    cgcAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCg                                               <  2:1201522/150‑1 (MQ=255)
                                                                                                     gcAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCcg                                                           <  1:1316784/137‑1 (MQ=255)
                                                                                                     gcAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCcg                                                           >  2:1316784/1‑137 (MQ=255)
                                                                                                     gcAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGACCAGCCgg                                              >  2:289031/1‑150 (MQ=255)
                                                                                                      cAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGt                                             >  1:561697/1‑150 (MQ=255)
                                                                                                       aaGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTc                                            <  2:553339/150‑1 (MQ=255)
                                                                                                         ggACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGa                                          >  2:1394270/1‑150 (MQ=255)
                                                                                                          gACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAg                                         <  1:540870/150‑1 (MQ=255)
                                                                                                           aCGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAgg                                        <  1:730537/150‑1 (MQ=255)
                                                                                                             gCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACctgc                                                                                     >  1:740167/1‑103 (MQ=255)
                                                                                                             gCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACctgc                                                                                     <  2:740167/103‑1 (MQ=255)
                                                                                                             gCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTc                                      >  1:256131/1‑150 (MQ=255)
                                                                                                              cccTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCa                                     >  2:1286827/1‑150 (MQ=255)
                                                                                                              cccTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGCCGAGGTcc                                     >  2:1088842/1‑149 (MQ=255)
                                                                                                                    tgAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGcgccg                               >  1:1369397/1‑150 (MQ=255)
                                                                                                                           gcCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCTGGCCCGCCgggt                                            <  1:207317/130‑3 (MQ=255)
                                                                                                                           gcCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCTGGCCCGCCgggt                                            >  2:207317/1‑128 (MQ=255)
                                                                                                                            ccGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAg                       >  1:360467/1‑150 (MQ=255)
                                                                                                                              ggATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCa                     <  1:1394270/150‑1 (MQ=255)
                                                                                                                               gATCAGGACGTGCGCCTGGAGGCGGAGCCAGGTTTGAGCTTGCAGGGTTGGGCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGCCGAGTTCAGCGCCGCGCGCgaggag                    >  2:7966/1‑150 (MQ=255)
                                                                                                                                aTCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAgg                   >  2:884657/1‑150 (MQ=255)
                                                                                                                                 tCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGa                  <  1:718374/150‑1 (MQ=255)
                                                                                                                                    ggACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCg               >  1:1208891/1‑150 (MQ=255)
                                                                                                                                         tGCGCCTCCAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCAAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCg          >  2:760602/1‑150 (MQ=255)
                                                                                                                                            gcCTGGAGGAGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTGa       <  1:1286827/150‑1 (MQ=255)
                                                                                                                                                 gAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAAGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTGATCAgc  <  2:561697/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CCTGGAAGACCTGATCAGCGAAATCCTCGCCCTGGCCCGGGTGGATGCCGAACAAGCCCATGCCGAACCGGTCGACCTCAATGCCCTACTCGGCAGCGTGCGCAAGGACGCCCTGCTGAGCGCGCCGGATCAGGACGTGCGCCTGGAGGCGCAGCCAGGTTTGAGCTTGCAGGGTTGGCCAACGCTGATCGAGCGTGCCGTTGACAACCTGCTGCGCAATGCCCTGCGCTTCAACCCGGCGGGCCAGCCGGTCGAGGTCAGCGCCGCGCGCGAGCAGGACCGCATCGTGATCAGC  >  NC_002947/5119099‑5119393

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