Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,074,213 A→C 6.8% intergenic (+248/‑182) PP_2683 → / → PP_2684 two‑component system sensor histidine kinase/response regulator/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,074,2130AC 6.8% 223.7 / 6.2 74intergenic (+248/‑182)PP_2683/PP_2684two‑component system sensor histidine kinase/response regulator/membrane protein
Reads supporting (aligned to +/- strand):  ref base A (35/34);  new base C (3/2);  total (38/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 = 8.82e-01

AAGAAGTTAGTAAACTAACTAATACGGCTTTAAATACCGCTGGTCGGTGCTTTTGCTCAAGAAATGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACTGTGTGCAATTGAGTAACCCTTTGCTGTTTCAATTAACGATATG  >  NC_002947/3074066‑3074360
                                                                                                                                                   |                                                                                                                                                   
aagaagTTAGTAAACTAACTAATACGGCTTTAAATACCGCTGGTCGGTGCTTTTGCTCAAGAAATGAACGATTTGCCAGTGGGTCGTTCATATATGCCAGCGAATGTGTTCTTCAAGTACAACACTTCGCACATCTCCCCAGAACGTAc                                                                                                                                                    >  3:263900/1‑149 (MQ=255)
   aagTTAGTAAACTAACTAATACGGCTTTAAATACCGCTGGTCGGTGCTTTTGCTCAAGAAATGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTgg                                                                                                                                                 <  6:189460/149‑1 (MQ=255)
         gTAAACTAACTAATACGGCTTTAAATACCGCTGGTCGGTGCTTTTGCTCAAGAAATGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGtt                                                                                                                                            >  5:170156/1‑148 (MQ=255)
         gTAAACTAACTAATACGGCTTTAAATACCGCTGGTCGGTGCTTTTGCTCAAGAAATGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGtt                                                                                                                                            <  2:30384/148‑1 (MQ=255)
         gTAAACTAACTAATACGGCTTTAAATACCGCTGGTCGGTGCTTTTGCTCAAGAAATGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTc                                                                                                                                           <  6:170156/149‑1 (MQ=255)
          tAAACTAACTAATACGGCTTTAAATACCGCTGGTCGGTGCTTTTGCTCAAGAAATGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTc                                                                                                                                           >  7:276590/1‑148 (MQ=255)
           aaaCTAACTAATACGGCTTTAAATACCGCTGGTCGGTGCTTTTGCTCAAGAAATGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCa                                                                                                                                          >  1:227680/1‑148 (MQ=255)
            aactaactATTACGGCGTTAGATACCGCTTGTCGGTGCTTTTGCTGAATAAATGAACGATTTGCCAGTGGGTGGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTAGCACACTTCGCACGTCTCGCCAGAAGGTACAGTACGTTctg                                                                                                                                         <  6:287102/148‑3 (MQ=255)
               taactaaTACGGCTTTAAATACCGCTGGTCGGTGCTTTTGCTCAAGAAATGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTAAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACtt                                                                                                                                      >  5:209842/1‑148 (MQ=255)
                aactaaTACGGCTTTAAATACCGCTGGTCGGTGCTTTTGCTCAAGAAATGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGa                                                                                                                                    >  3:208003/1‑149 (MQ=255)
                   taaTACGGCTTTAAATACCGCTGGTCGGTGCTTTTGCTCAAGAAATGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGAttt                                                                                                                                 >  7:107386/1‑149 (MQ=255)
                      tACGGCTTTAAATACCGCTGGTCGGTGCTTTTGCTCAAGAAATGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGaa                                                                                                                              <  4:185924/149‑1 (MQ=255)
                                  tACCGCTGGTCGGTGCTTTTGCTCAAGAAATGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATg                                                                                                                   >  3:221103/1‑148 (MQ=255)
                                     cGCTGGTCGGTGCTTTTGCTCAAGAAATGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCa                                                                                                                         <  1:1093/139‑1 (MQ=255)
                                     cGCTGGTCGGTGCTTTTGCTCAAGAAATGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCa                                                                                                                         >  2:1093/1‑139 (MQ=255)
                                              gTGCTTTTGCTCAAGAAATGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCtttt                                                                                                      >  3:171055/1‑149 (MQ=255)
                                               tGCTTTTGCTCAAGAAATGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCttttt                                                                                                     <  1:248835/149‑1 (MQ=255)
                                                  ttttGCTCAAGAAATGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTAc                                                                                                  >  5:162563/1‑149 (MQ=255)
                                                    ttGCTCAAGAAATGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTAc                                                                                                  <  5:169463/147‑1 (MQ=255)
                                                    ttGCTCAAGAAATGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTAc                                                                                                  <  4:192785/147‑1 (MQ=255)
                                                      gCTCAAGAAATGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACctc                                                                                               <  8:276590/148‑1 (MQ=255)
                                                       cTCAAGAAATGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACctct                                                                                              <  4:208003/148‑1 (MQ=255)
                                                        tCAAGAAATGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGa                                                                                            >  7:197738/1‑149 (MQ=255)
                                                          aaGAAATGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGAc                                                                                           >  5:34500/1‑148 (MQ=255)
                                                          aaGAAATGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGAc                                                                                           >  3:56152/1‑148 (MQ=255)
                                                          aaGAAATGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGAc                                                                                           >  8:291886/1‑148 (MQ=255)
                                                            gAAATGCACGAGTTGCCAGTGGGTCGTTCCTATTTGCCAGCTAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGATTCGCAGGGATGCACACGGGCTTTTTTACCTCTGTCAAt                                                                                        <  4:20062/149‑1 (MQ=255)
                                                             aaaTGATCGTTTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATg                                                                                       <  4:263900/149‑1 (MQ=255)
                                                             aaaTGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATg                                                                                       <  7:291886/149‑1 (MQ=255)
                                                              aaTTAACGTTTTGCCTGTGGGGCGTTCCTATTTGCCATGGTATTTGTTCTTCAAGTACCACACTTCGCAACTCTCCCCAGTACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATg                                                                                       <  8:229650/148‑1 (MQ=255)
                                                               aTGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTAAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGa                                                                                     <  6:209842/149‑1 (MQ=255)
                                                                 gAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAg                                                                                    >  3:20586/1‑148 (MQ=255)
                                                                    cGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAAc                                                                                <  4:56152/149‑1 (MQ=255)
                                                                     gATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATg                                                                                                                   <  1:110955/113‑1 (MQ=255)
                                                                     gATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATg                                                                                                                   >  2:110955/1‑113 (MQ=255)
                                                                     gATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACg                                                                               >  7:16930/1‑149 (MQ=255)
                                                                             aGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTc                                                                        >  7:251392/1‑148 (MQ=255)
                                                                                                ccAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCaa                                                          >  1:336462/1‑143 (MQ=255)
                                                                                                ccAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCaa                                                          <  2:336462/143‑1 (MQ=255)
                                                                                                ccAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGtt                                                     >  8:249005/1‑148 (MQ=255)
                                                                                                 cAGCGAATGTGTTCTTCGAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATTGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGtt                                                     >  5:64543/1‑147 (MQ=255)
                                                                                                     gAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGGCa                                               >  4:254090/1‑149 (MQ=255)
                                                                                                      aaTGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCAc                                              >  8:186230/1‑149 (MQ=255)
                                                                                                        tgtgTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCAcag                                            >  5:160971/1‑147 (MQ=255)
                                                                                                          tgtTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACtgt                                           <  8:197738/148‑1 (MQ=255)
                                                                                                          tgtTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACtgt                                           <  6:34500/148‑1 (MQ=255)
                                                                                                           gtTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATg                                                                                       >  1:1386/1‑103 (MQ=255)
                                                                                                           gtTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATg                                                                                       <  2:1386/103‑1 (MQ=255)
                                                                                                           gtTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACtgtgt                                         >  8:177225/1‑149 (MQ=255)
                                                                                                                tcAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACTGTGTGCAAt                                    >  4:80355/1‑149 (MQ=255)
                                                                                                                 cAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACTGTGTGCAAtt                                   >  8:141428/1‑149 (MQ=255)
                                                                                                                   aGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATTGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACTGTGTGCAATTGa                                 <  6:64543/149‑1 (MQ=255)
                                                                                                                   aGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACTGTGTGCAATTGa                                 <  7:249005/149‑1 (MQ=255)
                                                                                                                     tACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACTGTGTGCAATTGAg                                >  5:29649/1‑148 (MQ=255)
                                                                                                                     tACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACTGTGTGCAATTGAGt                               <  2:65822/149‑1 (MQ=255)
                                                                                                                       ccACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAg                                                                >  4:94386/1‑114 (MQ=255)
                                                                                                                       ccACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAg                                                                <  3:94386/114‑1 (MQ=255)
                                                                                                                        cacaCTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACTGTGTGCAATTGAGTAAc                            <  8:16930/149‑1 (MQ=255)
                                                                                                                                  aCCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACTGTGTGCAATTGAGTAACCCTTTGCTGt                  <  7:186230/149‑1 (MQ=255)
                                                                                                                                    ctctCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACTGTGTGCAATTGAGAAACCCTTTGCTGttt                <  4:221103/149‑1 (MQ=255)
                                                                                                                                      ctcCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACTGTGTGCAATTGAGTAACCCTTTGCTGTTTc               >  2:98896/1‑148 (MQ=255)
                                                                                                                                        ccccAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACTGTGTGCAATTGAGTAACCCTTTGCTGTTTCAAt            <  4:20586/149‑1 (MQ=255)
                                                                                                                                        aagttgAACGTCCAGTACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACTGTGTGCAATTGAGTAACCCTTTGCTGTTTCAAt            >  8:180447/6‑149 (MQ=255)
                                                                                                                                        aagttgAACGTCCAGTACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCAATGTCCAACCGAGGCTCAACGGTTGAGCCCCTGTGTGAAATTGAGTAACCCTTTTGTGTTTCAAt            >  8:185520/6‑149 (MQ=255)
                                                                                                                                         cccAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACTGTGTGCAATTGAGTAACCCTTTGCTGTTTCAAtt           <  6:29649/149‑1 (MQ=255)
                                                                                                                                         cccAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACTGTGTGCAATTGAGTAACCCTTTGCTGTTTCAAt            >  5:336696/1‑148 (MQ=255)
                                                                                                                                         cccAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTAAACGGTTGAGCCACTGTGTGCAATTGAGTAACCCTTTGCTGTTTCAAtt           >  8:321148/1‑149 (MQ=255)
                                                                                                                                         agttgAACGTCCAGTACGTTCAACTTGATTTGAATCGCAGGTATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACTGTGTGCAATTGAGTAACCCTTTGCTGTTTCAAtt           <  7:185520/145‑1 (MQ=255)
                                                                                                                                         agttgAACGTCCAGTACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACTGTGTGCAATTGAGTAACCCTTTGCTGTTTCAAtt           <  7:180447/145‑1 (MQ=255)
                                                                                                                                          gttgAACGTCCAGTACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACTGTGTGCAATTGAGTAACCCTTTGCTGTTTCAAtt           >  2:115463/4‑148 (MQ=255)
                                                                                                                                               aCGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACTGTGTGCAATTGa                                 <  3:119171/121‑1 (MQ=255)
                                                                                                                                               aCGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACTGTGTGCAATTGa                                 >  4:119171/1‑121 (MQ=255)
                                                                                                                                               aCGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACTGTGTGCAATTGAGTAACCCTTTGCTGTTTCAATTAACGat     <  3:300094/149‑1 (MQ=255)
                                                                                                                                               aCGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCCCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACCGTGTGCAATTGAGTAACCCTTTGCTGTTTCAATTAACGat     <  2:227680/149‑1 (MQ=255)
                                                                                                                                                   aCTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACTGTGTGCAATTGAGTAACCCTTTGCTGTTTCAATTAACGATATg  >  2:119990/1‑148 (MQ=255)
                                                                                                                                                   aCTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACTGTGTGCAATTGAGTAACCCTTTGCTGTTTCAATTAACGATATg  >  7:140322/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAGAAGTTAGTAAACTAACTAATACGGCTTTAAATACCGCTGGTCGGTGCTTTTGCTCAAGAAATGAACGATTTGCCAGTGGGTCGTTCCTATTTGCCAGCGAATGTGTTCTTCAAGTACCACACTTCGCACCTCTCCCCAGAACGTACTGGACGTTCAACTTGATTTGAATCGCAGGGATGCACACGGGCTTTTTTACCTCTGACAATGGAGTAACGCCATGTCCAACCGAGGCTCAACGGTTGAGCCACTGTGTGCAATTGAGTAACCCTTTGCTGTTTCAATTAACGATATG  >  NC_002947/3074066‑3074360

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: