Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,336,191 T→C 5.3% V874A (GTG→GCG)  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,1910TC 5.3% 229.4 / 2.2 75V874A (GTG→GCG) PP_4695two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (33/38);  new base C (2/2);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.08e-01

CAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAACATCATGGACCGCCTGTT  >  NC_002947/5336045‑5336338
                                                                                                                                                  |                                                                                                                                                    
cAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGCTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGc                                                                                                                                                    >  5:29368/1‑149 (MQ=255)
 aGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCt                                                                                                                                                   <  5:361502/149‑1 (MQ=255)
   cGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGa                                                                                                                                                 <  1:185197/149‑1 (MQ=255)
   cGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGa                                                                                                                                                 <  2:303468/149‑1 (MQ=255)
   cGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGa                                                                                                                                                 <  1:375613/149‑1 (MQ=255)
   cGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGa                                                                                                                                                 <  7:328561/149‑1 (MQ=255)
   cGAGGAGCCCGGGTGCCTGGCTGAAGTGGCACAAGATGCCAGCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCACCTGGGCCAGGCGCTgc                                                                                                                                                 <  2:225908/149‑2 (MQ=255)
      ggagCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTATAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATc                                                                                                                                               >  4:23199/1‑148 (MQ=255)
        agCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAAc                                                                                                                                            <  4:110021/149‑1 (MQ=255)
         gcacGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAATGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAAcc                                                                                                                                           <  3:315283/146‑1 (MQ=255)
            cGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGc                                                                                                                                        >  4:262278/1‑149 (MQ=255)
                  ccTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCACCTGGGCCAGGCGCTgc                                                                                                                                                 <  7:239324/134‑2 (MQ=255)
                  ccTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCACCTGGGCCAGGCGCTgc                                                                                                                                                 >  8:239324/1‑133 (MQ=255)
                      gCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTTGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGc                                                                                                                              <  3:178070/149‑1 (MQ=255)
                        tGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGccc                                                                                                                            >  5:305423/1‑149 (MQ=255)
                            gTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTATAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGa                                                                                                                        <  3:23199/149‑1 (MQ=255)
                            gTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGa                                                                                                                        <  1:144977/149‑1 (MQ=255)
                             tGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGAt                                                                                                                       <  1:336383/149‑1 (MQ=255)
                               gCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGc                                                                                                                     >  3:346308/1‑149 (MQ=255)
                                 acaAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCACCGCCCGCGATGCCt                                                                                                                   >  8:72140/1‑149 (MQ=255)
                                 acaAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCt                                                                                                                   >  5:358210/1‑149 (MQ=255)
                                    aGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTcgc                                                                                                                <  7:23156/149‑1 (MQ=255)
                                        gCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCgccgcc                                                                                                            >  5:90677/1‑149 (MQ=255)
                                         ccATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGccc                                                                                                           <  7:122436/149‑1 (MQ=255)
                                                              aCCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCt                                                                                                                                          <  3:102840/97‑1 (MQ=255)
                                                              aCCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCt                                                                                                                                          >  4:102840/1‑97 (MQ=255)
                                                                cGGCGCAATTTCGAAGTACAGTTCATCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAg                                                                                   >  7:51738/1‑149 (MQ=255)
                                                                  gcgcAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGc                                                                                  >  1:74172/1‑148 (MQ=255)
                                                                  gcgcAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCg                                                                                 <  7:292598/149‑1 (MQ=255)
                                                                   cgcAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGc                                                                                  >  3:233661/1‑147 (MQ=255)
                                                                    gcAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCCTTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAg                                                                               <  7:168847/149‑1 (MQ=255)
                                                                    gcAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAg                                                                               <  7:72140/149‑1 (MQ=255)
                                                                    gcAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAg                                                                               >  4:110932/1‑149 (MQ=255)
                                                                       aTTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGt                                                                             >  3:352476/1‑148 (MQ=255)
                                                                         ttCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGGCCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTaa                                                                          >  5:85368/1‑149 (MQ=255)
                                                                          tCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAAc                                                                         >  7:314511/1‑149 (MQ=255)
                                                                            gAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGa                                                                       <  4:346308/149‑1 (MQ=255)
                                                                            gAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGa                                                                       <  1:338627/149‑1 (MQ=255)
                                                                                    gTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGcaca                                                                  >  4:53467/1‑146 (MQ=255)
                                                                                    gTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACAcc                                                                <  1:287779/148‑1 (MQ=255)
                                                                                    gTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACAcc                                                                >  2:287779/1‑148 (MQ=255)
                                                                                        ttcAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGa                                                           >  5:286014/1‑149 (MQ=255)
                                                                                        ttcAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGa                                                           >  1:78823/1‑149 (MQ=255)
                                                                                             ccTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGGTCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGa                                                      <  1:263312/149‑1 (MQ=255)
                                                                                                 tGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGAAGCGAGGTTAACGAGCACACCGTCGACCTGAtc                                                    <  6:209786/147‑1 (MQ=255)
                                                                                                   cGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGAtcgtcg                                                <  3:53467/149‑1 (MQ=255)
                                                                                                     aCCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGAtcgtcg                                                >  2:107104/1‑147 (MQ=255)
                                                                                                     aCCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGaa                                              <  6:29368/149‑1 (MQ=255)
                                                                                                     aCCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGaa                                              <  6:305423/149‑1 (MQ=255)
                                                                                                     aCCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGaa                                              <  4:84801/149‑1 (MQ=255)
                                                                                                       ccGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGa                                            <  6:90677/149‑1 (MQ=255)
                                                                                                           aCCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGaa                                        <  4:352476/149‑1 (MQ=255)
                                                                                                                 aggCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGc                                  >  2:14479/2‑149 (MQ=255)
                                                                                                                  ggCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCg                                 <  6:358210/149‑1 (MQ=255)
                                                                                                                     cGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGa                              >  6:103944/1‑149 (MQ=255)
                                                                                                                      gAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCCGCGGGAt                             >  7:212/1‑149 (MQ=255)
                                                                                                                      gAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGAt                             <  6:286014/149‑1 (MQ=255)
                                                                                                                       aaGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGAtt                            >  6:46715/1‑149 (MQ=255)
                                                                                                                            gttgatcAGCACCTGGGCCAGGCGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCg                         >  1:246289/7‑147 (MQ=255)
                                                                                                                            gACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCg                         >  7:312723/1‑147 (MQ=255)
                                                                                                                               ccGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCgaagaa                    >  7:147629/1‑149 (MQ=255)
                                                                                                                                 gCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAAca                  <  6:77717/149‑1 (MQ=255)
                                                                                                                                 gCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAAca                  >  5:161237/1‑149 (MQ=255)
                                                                                                                                 gCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAAc                   >  3:300374/1‑148 (MQ=255)
                                                                                                                                  cAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAAcat                 <  2:78823/149‑1 (MQ=255)
                                                                                                                                   aGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAAcatc                <  5:103944/149‑1 (MQ=255)
                                                                                                                                   aGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAAcat                 >  5:239172/1‑148 (MQ=255)
                                                                                                                                      gcatGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCGTTGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAAcatcat              >  1:298893/4‑149 (MQ=255)
                                                                                                                                      gcCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAACATCATg             <  6:22468/149‑1 (MQ=255)
                                                                                                                                       ccTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGCAGAAGGCAGCGGGATTCCCAAGAAaatcagg             >  4:10215/1‑146 (MQ=255)
                                                                                                                                        cTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAACATCATGGa           <  3:262278/149‑1 (MQ=255)
                                                                                                                                             ccAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGcaca                                                                  <  8:107580/89‑1 (MQ=255)
                                                                                                                                             ccAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGcaca                                                                  >  7:107580/1‑89 (MQ=255)
                                                                                                                                                 gTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAACATCATGGACCGCCTGtt  <  4:233661/149‑1 (MQ=255)
                                                                                                                                                 gTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAACATCATGGACCGCCTGtt  <  3:110932/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CAGCGAGGAGCCCGTGTGCCTGGCTGAAGTGGCACAAGATGCCATCGGTCTGCTGGCGTTGAACCGGCGCAATTTCGAAGTACAGTTCTTCAACCTTTGCGACCCGGACCATTGGGCCGAAGGTGACCCGCAGCGCCTGGCCCAGGTGCTGATCAACCTGCTCTCCAACGCCCGCGATGCCTCGCCGCCCGGCAGCGCCG‑TGCGCGTGCGCAGCGAGGTTAACGAGCACACCGTCGACCTGATCGTCGAAGACGAAGGCAGCGGGATTCCGAAGAACATCATGGACCGCCTGTT  >  NC_002947/5336045‑5336338

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