Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,719,164 A→G 6.2% E621E (GAA→GAG)  PP_5020 → methyl‑accepting chemotaxis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,719,1640AG 6.2% 238.9 / 4.3 81E621E (GAA→GAG) PP_5020methyl‑accepting chemotaxis protein
Reads supporting (aligned to +/- strand):  ref base A (38/38);  new base G (3/2);  total (41/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 = 9.43e-01

GATCCAGGCGATGAACCAGCAGATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAATCTCCACCCGCCGTTCACCATT  >  NC_002947/5719017‑5719311
                                                                                                                                                   |                                                                                                                                                   
gATCCAGGCGATGAACCAGCAGATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGaa                                                                                                                                                     <  6:96493/148‑1 (MQ=255)
     aGGCGATGAACCAGCAGATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGgcg                                                                                                                                               >  7:30198/1‑149 (MQ=255)
         gATGAACCAGCAGATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGt                                                                                                                                           >  2:272340/1‑149 (MQ=255)
         gATGAACCAGCAGATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGt                                                                                                                                           >  3:86093/1‑149 (MQ=255)
          aTGAACCAGCAGATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGtt                                                                                                                                          <  5:222548/149‑1 (MQ=255)
           tGAACCAGCAGATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTg                                                                                                                                         <  4:369357/149‑1 (MQ=255)
             aaCCAGCAGATTGCCACGGCGGGAGAGGAGAAGACGGCGGTGGCCAAAGAAATCAACCGCAGCATGATGAATGTGCGTAATATTTCCGACAAGACCACGGCGGCCAGTGAGGAGACGGCCAGCTCAAGCGTGGAACTGGCGCGGTTggt                                                                                                                                       >  7:56077/1‑148 (MQ=255)
             aaCCAGCAGATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTgg                                                                                                                                        <  1:75698/148‑1 (MQ=255)
             aaCCAGCAGATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTgg                                                                                                                                        <  4:136682/148‑1 (MQ=255)
             aaCCAGCAGATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTgg                                                                                                                                        <  4:267220/148‑1 (MQ=255)
                cagcagATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTAc                                                                                                                                    <  3:273783/149‑1 (MQ=255)
                cagcagATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGAGGCCAGTGAGGAGACGGCCAGCTCCAGCTTGGAACTGGCGCGGTTGGGTAc                                                                                                                                    >  5:222361/1‑149 (MQ=255)
                 agcagATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTAcc                                                                                                                                   >  7:355001/1‑149 (MQ=255)
                      aTTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCAttt                                                                                                                              >  2:292984/1‑149 (MQ=255)
                         gCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGtt                                                                                                                                          <  4:356458/134‑1 (MQ=255)
                         gCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGtt                                                                                                                                          >  3:356458/1‑134 (MQ=255)
                          cgaCGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCaa                                                                                                                          <  6:117817/147‑1 (MQ=255)
                          ccACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTgg                                                                                                                            >  3:20562/1‑146 (MQ=255)
                          ccACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGc                                                                                                                            >  7:239188/1‑147 (MQ=255)
                           cACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCGAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAg                                                                                                                         <  3:46559/149‑1 (MQ=255)
                            aCGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAgg                                                                                                                        <  1:301245/149‑1 (MQ=255)
                              ggcggcCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGggt                                                                                                                      <  1:157675/149‑1 (MQ=255)
                                   ccGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGggttggtt                                                                                                                 <  4:88282/149‑1 (MQ=255)
                                   ccGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGATCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGggttggtt                                                                                                                 >  8:97641/1‑149 (MQ=255)
                                    cGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGgttggttg                                                                                                                >  4:188205/1‑149 (MQ=255)
                                          gCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGt                                                                                                          <  3:188205/149‑1 (MQ=255)
                                             gACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGtt                                                                                                         >  2:6058/1‑147 (MQ=255)
                                               cggcggTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAgg                                                                                                     <  8:289127/149‑1 (MQ=255)
                                                   ggTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCtgt                                                                                                 <  6:224902/149‑1 (MQ=255)
                                                        ccGAAGAAATGAACCGCAGCGTGATGAATGTGCGTGATATTTCCGAGCAGACCTGGGCGGCCAGTGTGGAGACGGCCAGGTCCAGCGTGGAACTGGCGCGGTTGGGTAACCATTTGCAAGTGTTTGTTGGGCGGTTCAGGCTGTGATCt                                                                                            <  3:213587/149‑1 (MQ=255)
                                                        ccGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCt                                                                                            >  1:365525/1‑149 (MQ=255)
                                                         cGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCt                                                                                            >  7:315023/1‑148 (MQ=255)
                                                         cGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGAACAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCt                                                                                            >  8:211484/1‑148 (MQ=255)
                                                           aagaaATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGtt                                                                                                         <  7:67497/133‑1 (MQ=255)
                                                           aagaaATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGtt                                                                                                         >  8:67497/1‑133 (MQ=255)
                                                            agaaATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAAc                                                                                        <  7:97641/149‑1 (MQ=255)
                                                               aaTCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCACACCGCGGCGGCCAGGGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAg                                                                                     <  1:181536/149‑1 (MQ=255)
                                                                aTCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGc                                                                                    <  1:272340/149‑1 (MQ=255)
                                                                 tCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCg                                                                                   >  5:273003/1‑149 (MQ=255)
                                                                 tCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCg                                                                                   >  8:113241/1‑149 (MQ=255)
                                                                 tCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCg                                                                                   >  6:191592/1‑149 (MQ=255)
                                                                     ccGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCg                                                                               >  3:358534/1‑149 (MQ=255)
                                                                      cGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGc                                                                              <  7:122096/149‑1 (MQ=255)
                                                                           cGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGt                                                                         <  5:330356/149‑1 (MQ=255)
                                                                             tgatgaATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGCACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCgg                                                                          >  7:346607/1‑146 (MQ=255)
                                                                             tgatgaATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGCACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCgg                                                                          <  8:346607/146‑1 (MQ=255)
                                                                               atgaATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCgg                                                                      <  7:121611/148‑1 (MQ=255)
                                                                               atgaATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCgg                                                                      <  7:113241/148‑1 (MQ=255)
                                                                               atgaATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCgg                                                                      >  8:121611/1‑148 (MQ=255)
                                                                                       gCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTcg                                                             >  7:197578/1‑149 (MQ=255)
                                                                                        cGTGATATTTCCGACCAGACCTCGTCGGCCAGTGATGAGACGGCCATCTCCAGCGTGGAACTGTCGCGGTTTGGTACCCATTTTCAAGGGTTGGTTTGGTGGTTCAGGCTGTGATCTGAACTAGCGGTCTTCCGGTTGGGCCTCTTcgc                                                            <  8:56077/149‑1 (MQ=255)
                                                                                         gTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTcgcg                                                           >  1:65706/1‑149 (MQ=255)
                                                                                          tGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCgg                                                          <  8:239188/149‑1 (MQ=255)
                                                                                                ttCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAAc                                                    >  8:139241/1‑149 (MQ=255)
                                                                                                 tCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTTGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAAc                                                    <  1:292984/148‑1 (MQ=255)
                                                                                                 tCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAAc                                                    <  6:273003/148‑1 (MQ=255)
                                                                                                   cGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCg                                                 >  2:136009/1‑149 (MQ=255)
                                                                                                    gACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGc                                                >  5:11405/1‑149 (MQ=255)
                                                                                                     aCCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGc                                                >  3:17872/1‑148 (MQ=255)
                                                                                                     aCCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGc                                                >  1:362944/1‑148 (MQ=255)
                                                                                                         gACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATg                                           >  8:42408/1‑149 (MQ=255)
                                                                                                          aCCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTATTgg                                          <  1:6058/149‑1 (MQ=255)
                                                                                                            cTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGc                                        <  6:11405/149‑1 (MQ=255)
                                                                                                            cTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGc                                        >  8:27957/1‑149 (MQ=255)
                                                                                                            cTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGc                                        >  6:298998/1‑149 (MQ=255)
                                                                                                             tCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGc                                        >  7:264646/1‑148 (MQ=255)
                                                                                                                gcggcCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAg                                    >  3:209089/1‑149 (MQ=255)
                                                                                                                     cAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGaa                                                                                         >  8:176376/1‑91 (MQ=255)
                                                                                                                     cAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGaa                                                                                         <  7:176376/91‑1 (MQ=255)
                                                                                                                     cAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGCCGCGGTTGGGTACCCATTTGCAAGGTTTGGTTGGGCGGTTCAGGCTGAGAACTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACACGAGAATGGGCCAAGCCgg                                >  8:320561/1‑148 (MQ=255)
                                                                                                                              gaCGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAAtc                      <  8:30198/149‑1 (MQ=255)
                                                                                                                                gcgCCAGTTCCACGCTGGAGCTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGCGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAAtctc                    >  6:221856/3‑149 (MQ=255)
                                                                                                                                 ggCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAAtctc                    >  3:287482/1‑148 (MQ=255)
                                                                                                                                  gCCAGTTCCACGCTGGAGCTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGcc                                       >  6:57304/1‑128 (MQ=255)
                                                                                                                                  gCCAGTTCCACGCTGGAGCTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGcc                                       <  5:57304/128‑1 (MQ=255)
                                                                                                                                  gCCAGTTCCACGCTGGAGCTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAATCTCCa                  >  5:364838/1‑149 (MQ=255)
                                                                                                                                          cacgctGGAGCTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAATCTCCACCCGCCGt          <  6:364838/144‑1 (MQ=255)
                                                                                                                                          cAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGGTGGGCGGTGCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAATCTCCACccgcc            <  5:74668/147‑1 (MQ=255)
                                                                                                                                            gCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAATCTCCACCCGCCGTTc        <  2:362944/149‑1 (MQ=255)
                                                                                                                                            gCGTGGAACTGGCGCGGGTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAATCTCCACCCGCCGTTc        <  4:17872/149‑1 (MQ=255)
                                                                                                                                             cGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAATCTCCACCCGCCGTTCa       <  4:209089/149‑1 (MQ=255)
                                                                                                                                             cGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAATCTCCACCCGCCGTTCa       <  2:365525/149‑1 (MQ=255)
                                                                                                                                                  aaCTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGTCGGCTGGCCAATCTCCACCCGCCGTTCACCAtt  <  7:42408/149‑1 (MQ=255)
                                                                                                                                                  aaCTGGCGCGGGTGGGTAACCATTTGGAATGGTTGGTTGGGGGGGTCATGGTGGGTTCTTAAATTTCGGTCGGCCGGTCGGGGCTGTTTGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAATCTCCACCCGCCGTTCACCAtt  <  8:315023/149‑1 (MQ=25)
                                                                                                                                                   |                                                                                                                                                   
GATCCAGGCGATGAACCAGCAGATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAATCTCCACCCGCCGTTCACCATT  >  NC_002947/5719017‑5719311

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: ‑