Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,240,144 T→A 5.2% H49L (CAT→CTT)  PP_3717 ← LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,240,1440TA 5.2% 240.0 / 2.4 76H49L (CAT→CTT) PP_3717LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (41/31);  new base A (2/2);  total (43/33)
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.94e-01

CCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGG  >  NC_002947/4240002‑4240292
                                                                                                                                              |                                                                                                                                                    
ccGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCa                                                                                                                                                <  8:42025/149‑1 (MQ=255)
     gCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGgcg                                                                                                                                           >  5:443328/1‑149 (MQ=255)
     gCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGgc                                                                                                                                            <  7:186538/148‑1 (MQ=255)
      ccGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGgc                                                                                                                                            <  6:427204/147‑1 (MQ=255)
        gTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGgcgcgc                                                                                                                                        >  3:275651/1‑149 (MQ=255)
        gTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGgcgcgc                                                                                                                                        >  7:299821/1‑149 (MQ=255)
         ttCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCt                                                                                                                                       >  4:390788/1‑149 (MQ=255)
           cGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTg                                                                                                                                      <  2:15583/148‑1 (MQ=255)
            ggggATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGc                                                                                                                                     <  2:199402/148‑1 (MQ=255)
             gggATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCaa                                                                                                                                   >  3:166262/1‑149 (MQ=255)
              ggATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAAc                                                                                                                                  <  8:143292/149‑1 (MQ=255)
                    cTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACg                                                                                                                            >  7:252268/1‑149 (MQ=255)
                        cTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGc                                                                                                                         <  7:178626/148‑1 (MQ=255)
                        cTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGc                                                                                                                         >  1:191864/1‑148 (MQ=255)
                         tGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCCCATACCCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCa                                                                                                                       >  4:135498/1‑149 (MQ=255)
                         tGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCa                                                                                                                       >  4:134145/1‑149 (MQ=255)
                           gTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTCGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGcc                                                                                                                        <  4:40264/146‑1 (MQ=255)
                               aGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCActgctg                                                                                                                 >  5:261184/1‑149 (MQ=255)
                                gCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACtgctgc                                                                                                                <  1:367479/149‑1 (MQ=255)
                                 cAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCActgctg                                                                                                                 >  4:115249/1‑147 (MQ=255)
                                  aaCTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAg                                                                                                              >  8:45658/1‑149 (MQ=255)
                                   aCTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAgc                                                                                                             >  6:108266/1‑149 (MQ=255)
                                           gcccgcccCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCAt                                                                                                     >  6:358655/1‑149 (MQ=255)
                                            cccgcccCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGAGCGCGCTCCAACGCCACGTGCCACTGCTGCAGCGCGCCCaat                                                                                                    >  6:186366/1‑147 (MQ=255)
                                              cgcccCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATtgt                                                                                                  >  3:387800/1‑149 (MQ=255)
                                                ccccTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGa                                                                                                <  1:56710/149‑1 (MQ=255)
                                                 cccTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAg                                                                                               >  8:282170/1‑149 (MQ=255)
                                                  ccTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAgg                                                                                              <  8:307642/149‑1 (MQ=255)
                                                            aaGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCggg                                                                                     <  7:282170/148‑1 (MQ=255)
                                                                   atcatcatcGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCgg                                                                              <  2:219241/148‑1 (MQ=255)
                                                                    tcatcatcGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCggg                                                                             <  4:166262/148‑1 (MQ=255)
                                                                    tcatcatcGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCgg                                                                              >  2:263399/1‑147 (MQ=255)
                                                                     catcatcGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACtgctgc                                                                                                                >  3:212756/1‑112 (MQ=255)
                                                                     catcatcGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACtgctgc                                                                                                                <  4:212756/112‑1 (MQ=255)
                                                                       tcatcGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCCCGTGCCACTGCTGCCGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGagag                                                                         >  5:268859/1‑149 (MQ=255)
                                                                         atcGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGc                                                                       >  8:79695/1‑149 (MQ=255)
                                                                            gACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCt                                                                     <  3:62983/148‑1 (MQ=255)
                                                                             aCTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTgg                                                                   <  8:118446/149‑1 (MQ=255)
                                                                             aCTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTgg                                                                   <  3:115249/149‑1 (MQ=255)
                                                                             aCTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTgg                                                                   <  5:145609/149‑1 (MQ=255)
                                                                                   ccTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATg                                                             >  7:390364/1‑149 (MQ=255)
                                                                                       cACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGccc                                                          <  5:108266/148‑1 (MQ=255)
                                                                                       cACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGcc                                                           <  3:119259/147‑1 (MQ=255)
                                                                                       cACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCg                                                         <  4:423781/149‑1 (MQ=255)
                                                                                        aCCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCgg                                                        >  1:324153/1‑149 (MQ=255)
                                                                                         ccAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCgg                                                        >  4:178881/1‑148 (MQ=255)
                                                                                            tgtGGCCGCTGCCATCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAg                                                     >  2:259408/3‑148 (MQ=255)
                                                                                            aCTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGt                                                    >  6:374706/1‑149 (MQ=255)
                                                                                            aCTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGt                                                    >  5:400864/1‑149 (MQ=255)
                                                                                                  cGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGcc                                              >  6:148150/1‑149 (MQ=255)
                                                                                                        cAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAg                                        >  5:432750/1‑149 (MQ=255)
                                                                                                         aGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAgg                                       <  2:191864/149‑1 (MQ=255)
                                                                                                          gCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCa                                                      <  8:216852/133‑1 (MQ=255)
                                                                                                          gCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCa                                                      >  7:216852/1‑133 (MQ=255)
                                                                                                                  aTACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCcgc                              >  4:118275/1‑149 (MQ=255)
                                                                                                                   tACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCcgcg                             <  5:374706/149‑1 (MQ=255)
                                                                                                                      tCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGca                           >  2:192674/1‑148 (MQ=255)
                                                                                                                      tCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGATGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCAAGCCCGAGca                           >  8:225489/1‑148 (MQ=255)
                                                                                                                        cAGCGCCTCGGGGGTTTCACCAAGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGccc                                >  3:265440/1‑141 (MQ=255)
                                                                                                                        cAGCGCCTCGGGGGTTTCACCAAGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGccc                                <  4:265440/141‑1 (MQ=255)
                                                                                                                        cAGCGCCTCGGGGGTTTCACCAAGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGccc                                >  3:270912/1‑141 (MQ=255)
                                                                                                                        cAGCGCCTCGGGGGTTTCACCAAGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGccc                                <  4:270912/141‑1 (MQ=255)
                                                                                                                         aGCGCCTCGGGGGTTTCACAATGGGACAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGc                       >  7:30987/1‑149 (MQ=255)
                                                                                                                              cTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCttt                  <  7:45658/149‑1 (MQ=255)
                                                                                                                                 gggggTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTgcg               >  8:339028/1‑149 (MQ=255)
                                                                                                                                 gggggTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTgc                <  6:443328/148‑1 (MQ=255)
                                                                                                                                    ggTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCa            <  6:432750/149‑1 (MQ=255)
                                                                                                                                     gTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGGGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAg           >  2:428832/1‑149 (MQ=255)
                                                                                                                                     gTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAg           <  4:275651/149‑1 (MQ=255)
                                                                                                                                      tttCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAgg          >  5:290855/1‑149 (MQ=255)
                                                                                                                                       ttCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCcgc                              <  4:163817/128‑1 (MQ=255)
                                                                                                                                       ttCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCcgc                              >  3:163817/1‑128 (MQ=255)
                                                                                                                                         cACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGtt        >  3:121363/1‑148 (MQ=255)
                                                                                                                                           ccATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCAt     >  3:335025/1‑149 (MQ=255)
                                                                                                                                            cATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATc    >  2:404057/1‑149 (MQ=255)
                                                                                                                                            cATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATc    <  6:261184/149‑1 (MQ=255)
                                                                                                                                              tGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCgg  >  1:117339/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
CCGATGCCGTTCGGGGATGCCTTGCTGGTACAGCAACTGCGGCGCCCGCCCCTGGCACTCAAGGCTGATCATCATCGACTCCACCTGCACCAACTGGCCGATGGCAGCGGCCACATACTCCAGCGCCTCGGGGGTTTCACCATGGGCCAGGGCGCGCTGCAACGCCACGTGCCACTGCTGCAGCGCGCCCATTGTGAGGGTAAGCGGGGTATCGGGAGAGGCCTGGGTCATGCCCGGCAGTCTAGCCTGCCAGGCACGCCCGCGCACAGCCCTTTGCGCCAGGTTCATCGG  >  NC_002947/4240002‑4240292

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: