Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,622,190 G→C 6.9% E368D (GAG→GAC) ‡ PP_1421 → sensor histidine kinase TctE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,622,1900GC 6.9% 223.5 / 5.5 72E368D (GAG→GAC) ‡PP_1421sensor histidine kinase TctE
Reads supporting (aligned to +/- strand):  ref base G (34/33);  new base C (2/3);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

AGGCGGTGCGCAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCG  >  NC_002947/1622043‑1622336
                                                                                                                                                    |                                                                                                                                                  
aGGCGGTGCGCAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAgt                                                                                                                                                   <  6:218078/149‑1 (MQ=255)
  gCGGTGCGCAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAgttg                                                                                                                                                 >  6:79647/1‑149 (MQ=255)
    ggTGCGCAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAgttgtt                                                                                                                                               <  5:79647/149‑1 (MQ=255)
    ggTGCGCAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAgttgtt                                                                                                                                               <  6:223188/149‑1 (MQ=255)
     gTGCGCAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGttgttg                                                                                                                                              <  7:96746/149‑1 (MQ=255)
      tGCGCAGGGTCTGGACTTGAGTCTGTTTGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGa                                                                                                                                             <  2:269107/149‑1 (MQ=255)
      tGCGCAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGa                                                                                                                                             <  6:19552/149‑1 (MQ=255)
        cgcAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGGGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACAACTCGTTGAg                                                                                                                                            <  4:155338/148‑1 (MQ=255)
           aGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAAc                                                                                                                                        >  4:75451/1‑149 (MQ=255)
                cTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCtggt                                                                                                                                   >  2:358017/1‑149 (MQ=255)
                 tGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCtggtg                                                                                                                                  >  6:27253/1‑149 (MQ=255)
                   gACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCtggt                                                                                                                                   <  3:214448/146‑1 (MQ=255)
                   gACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGa                                                                                                                                <  2:368533/149‑1 (MQ=255)
                   gACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGa                                                                                                                                >  3:146029/1‑149 (MQ=255)
                     cTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGc                                                                                                                                           <  6:273891/136‑1 (MQ=255)
                     cTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGc                                                                                                                                           >  5:273891/1‑136 (MQ=255)
                      ttGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACaa                                                                                                                             >  4:83600/1‑149 (MQ=255)
                        gAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAAt                                                                                                                            <  1:67239/148‑1 (MQ=255)
                        gAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACAACTCGTTGAGCa                                                                                                                                          <  1:221494/134‑1 (MQ=255)
                        gAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACAACTCGTTGAGCa                                                                                                                                          >  2:221494/1‑134 (MQ=255)
                          gTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACaa                                                                                                                             >  2:56815/1‑145 (MQ=255)
                          gTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACaa                                                                                                                             <  1:56815/145‑1 (MQ=255)
                             aGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTg                                                                                                                      <  3:331186/149‑1 (MQ=255)
                              gTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACAACTCGTTGAGCaatg                                                                                                                                       <  3:179326/131‑3 (MQ=255)
                                    ccGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATa                                                                                                               <  4:325776/149‑1 (MQ=255)
                                     cgcgAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATAc                                                                                                              >  3:23735/1‑149 (MQ=255)
                                      gcgAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACg                                                                                                             >  1:104561/1‑149 (MQ=255)
                                       cgAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGc                                                                                                            <  5:27253/149‑1 (MQ=255)
                                         aGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCa                                                                                                          <  2:330300/149‑1 (MQ=255)
                                              ggCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCgg                                                                                                      <  3:351449/148‑1 (MQ=255)
                                               gCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGCGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGt                                                                                                    >  6:254158/1‑149 (MQ=255)
                                               gCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGt                                                                                                    >  6:314653/1‑149 (MQ=255)
                                                catggccatggcGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTg                                                                                                   >  6:255507/1‑149 (MQ=255)
                                                 atggccatggcGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTgg                                                                                                  <  3:83600/149‑1 (MQ=255)
                                                 atggccatggcGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTgg                                                                                                  <  2:104561/149‑1 (MQ=255)
                                                       atggcGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAAtgt                                                                                            <  8:33511/149‑1 (MQ=255)
                                                            gcCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGAttt                                                                                       <  8:184164/149‑1 (MQ=255)
                                                              cTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGCGGCCGAGTCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCACCGAGTTGTGGCGCAACCTGGTGGCCAATGACGTGGCGCATACGCAAGCGGGGGGCAATGTGAttttt                                                                                     >  7:351272/1‑148 (MQ=255)
                                                              cTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTg                                                                                     <  4:313675/149‑1 (MQ=255)
                                                              cTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTg                                                                                     >  2:129043/1‑149 (MQ=255)
                                                               tttGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAAGGTGATTTTGc                                                                                    >  5:220713/1‑149 (MQ=255)
                                                                   gcgcACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAAt                                                                                              <  4:265418/135‑1 (MQ=255)
                                                                   gcgcACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAAt                                                                                              >  3:265418/1‑135 (MQ=255)
                                                                     gcACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTaa                                                                              >  6:177692/1‑149 (MQ=255)
                                                                      cacaAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAAt                                                                             <  4:146029/149‑1 (MQ=255)
                                                                      cacaAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAAt                                                                             >  4:91531/1‑149 (MQ=255)
                                                                        caAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATgg                                                                           <  8:240791/149‑1 (MQ=255)
                                                                           gcgcGGGGTGGCGCTGGCGCTGGAGGCCGGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCTTACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGgcg                                                                         >  1:22435/1‑149 (MQ=255)
                                                                                   tggcgctggcgcTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTg                                                                <  3:75451/149‑1 (MQ=255)
                                                                                       gctggcgctggAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtgg                                                            <  3:220755/149‑1 (MQ=255)
                                                                                        ctggcgctggAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtgga                                                           <  8:8937/149‑1 (MQ=255)
                                                                                           gcgctggAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtggaagt                                                        <  2:182070/149‑1 (MQ=255)
                                                                                            cgctggAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtggaagtg                                                       <  4:23735/149‑1 (MQ=255)
                                                                                              ctggAGGCCGAGGCGCCGGTGTGGTTGAAGGGGTATAATACATTGCTCAACGAGTTGTTGAGCAACCTTGTGGACAATGCCTTTGAGCATACGCAAGAGGGTGGAAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtggaagtgtt                                                     <  4:304494/149‑3 (MQ=255)
                                                                                                 gAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCtggaagtggaa                                                    <  3:91531/147‑1 (MQ=255)
                                                                                                 gAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAga                                                  <  8:206869/149‑1 (MQ=255)
                                                                                                 gAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAga                                                  >  2:78024/1‑149 (MQ=255)
                                                                                                    gCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACTAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAgacga                                               >  6:273647/1‑149 (MQ=255)
                                                                                                        gAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACgg                                            >  3:276417/1‑149 (MQ=255)
                                                                                                        gAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACgg                                            <  5:254158/149‑1 (MQ=255)
                                                                                                               cGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAgg                                      >  7:134959/1‑148 (MQ=255)
                                                                                                                ggTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTa                                    >  6:364160/1‑149 (MQ=255)
                                                                                                                                gAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAgcg                    >  6:82628/1‑149 (MQ=255)
                                                                                                                                 aGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAgcgc                   >  6:284140/1‑149 (MQ=255)
                                                                                                                                  gCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAgcgc                   >  2:130479/1‑148 (MQ=255)
                                                                                                                                     gACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCGGTGGCCGAGCGCGAAc               >  7:128358/1‑149 (MQ=255)
                                                                                                                                       cATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACgg             <  5:364160/149‑1 (MQ=255)
                                                                                                                                       cATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTCGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACgg             >  3:260558/1‑149 (MQ=255)
                                                                                                                                         ttGCTCAACAACTCGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGgt           >  6:335602/1‑149 (MQ=255)
                                                                                                                                           gCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGgtgt         >  7:108137/1‑149 (MQ=255)
                                                                                                                                               aaCGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGt                                                                                                    >  4:355538/1‑54 (MQ=255)
                                                                                                                                               aaCGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGt                                                                                                    <  3:355538/54‑1 (MQ=255)
                                                                                                                                               aaCGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGACCGGGTGTTc       >  6:328703/1‑147 (MQ=255)
                                                                                                                                               aaCGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTc       >  3:306842/1‑147 (MQ=255)
                                                                                                                                               aaCGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTc       <  1:358017/147‑1 (MQ=255)
                                                                                                                                                  gAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCg  >  8:192177/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGGCGGTGCGCAGCGTCTGGACTTGAGTCAGTTGGCCCGCGAGCTGGGCATGGCCATGGCGCCTTTGGCGCACAAGCGCGGGGTGGCGCTGGCGCTGGAGGCCGAGGCGCCGGTGTGGTTGAAGGGGGAGCCGACATTGCTCAACGAGTTGTTGAGCAACCTGGTGGACAATGCCTTGGCGCATACGCCAGCGGGTGGCAATGTGATTTTGCGGGTAATGGCGCCGGCAGTGCTGGAAGTGGAAGACGACGGGCCAGGTATTCCCGTGGCCGAGCGCGAACGGGTGTTCGAGCG  >  NC_002947/1622043‑1622336

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