Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,910,411 G→A 6.4% P134L (CCG→CTG)  PP_1710 ← phthalate permease family MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,910,4110GA 6.4% 184.7 / 2.5 62P134L (CCG→CTG) PP_1710phthalate permease family MFS transporter
Reads supporting (aligned to +/- strand):  ref base G (29/29);  new base A (2/2);  total (31/31)
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.46e-01

AGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGCCCTGGGCCACGGTGACCGCCGACCAGCTGATCAGCGACACGCCATAGGTCAGGCGTGAGCCGAAGCGGTCGGCAATGAAGC  >  NC_002947/1910265‑1910555
                                                                                                                                                  |                                                                                                                                                
aGTCATAGGTGGTCGCCAGCCACAGCAGCGCACAGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCggg                                                                                                                                                >  1:265763/1‑149 (MQ=255)
      aGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGACGGCCGGGGCCTc                                                                                                                                           >  3:114762/1‑148 (MQ=255)
      aGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAAACGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTc                                                                                                                                           >  7:272572/1‑148 (MQ=255)
       ggtggtCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCa                                                                                                                                         <  3:285207/149‑1 (MQ=255)
        gtggtCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCAt                                                                                                                                        >  2:253178/1‑149 (MQ=255)
           gtCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGgc                                                                                                                                     <  3:174624/149‑1 (MQ=255)
            tCGCCAGCCACAGCAGCTCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGgcg                                                                                                                                    <  6:280146/149‑1 (MQ=255)
              gccagccaCAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGcc                                                                                                                                  <  7:105744/149‑1 (MQ=255)
               ccagccaCAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCg                                                                                                                                 >  2:187226/1‑149 (MQ=255)
                         cagcGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGgc                                                                                                                       <  3:299554/149‑1 (MQ=255)
                         cagcGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGgc                                                                                                                       >  2:309244/1‑149 (MQ=255)
                         cagcGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGgc                                                                                                                       <  1:95775/149‑1 (MQ=255)
                         cagcGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCGGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGgc                                                                                                                       <  3:267523/149‑1 (MQ=255)
                          agcGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGgcg                                                                                                                      <  2:181193/149‑1 (MQ=255)
                             gcACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAg                                                                                                                   <  8:272572/149‑1 (MQ=255)
                                  ggTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCg                                                                                                               >  3:304001/1‑148 (MQ=255)
                                   gTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGaa                                                                                                             >  5:130006/1‑149 (MQ=255)
                                   gTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGaa                                                                                                             >  6:188316/1‑149 (MQ=255)
                                   gTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGaa                                                                                                             >  4:312658/1‑149 (MQ=255)
                                   gTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGaa                                                                                                             >  3:169689/1‑149 (MQ=255)
                                       aCAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGt                                                                                                         <  5:188316/149‑1 (MQ=255)
                                       aCAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGt                                                                                                         <  1:253178/149‑1 (MQ=255)
                                        cAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGt                                                                                                         >  7:234743/1‑148 (MQ=255)
                                        cAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTg                                                                                                        <  4:12354/149‑1 (MQ=255)
                                           cgcgGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACa                                                                                                           >  3:202224/1‑143 (MQ=255)
                                           cgcgGTGCCCAGGTACTGGCAACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACa                                                                                                           <  4:202224/143‑1 (MQ=255)
                                            gcgGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCa                                                                                                                          >  1:274000/1‑127 (MQ=255)
                                            gcgGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCa                                                                                                                          <  2:274000/127‑1 (MQ=255)
                                              ggTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGa                                                                                                  <  2:265763/149‑1 (MQ=255)
                                                  cccAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGcc                                                                                              >  6:319102/1‑149 (MQ=255)
                                                  cccAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGcc                                                                                              >  6:260053/1‑149 (MQ=255)
                                                    cAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAg                                                                                                      >  2:312406/1‑139 (MQ=255)
                                                    cAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAg                                                                                                      <  1:312406/139‑1 (MQ=255)
                                                    cAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGc                                                                                            >  7:26196/1‑149 (MQ=255)
                                                     aGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACa                                                                                                           >  6:236330/1‑133 (MQ=255)
                                                     aGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACa                                                                                                           <  5:236330/133‑1 (MQ=255)
                                                        tACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGc                                                                                        <  8:234743/149‑1 (MQ=255)
                                                                                    ccATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGCCCTGGGCCACGGTGACCGCCGAc                                                            <  1:187226/149‑1 (MQ=255)
                                                                                         ccGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGCCCTGGGCCACGGTGACCGCCGACCAGc                                                        <  3:312658/148‑1 (MQ=255)
                                                                                         ccGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGCCCTGGGCCACGGTGACCGCCGACCAGc                                                        <  2:114343/148‑1 (MQ=255)
                                                                                           gcgcTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGCCCTGGGCCACGGTGACCGCCGACCAGCTGa                                                     <  8:26196/149‑1 (MQ=255)
                                                                                             gcTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGCCCTGGGCCACGGTGACCGCCGACCAGCTGAt                                                    >  4:196270/1‑148 (MQ=255)
                                                                                              cTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGCCCTGGGCCACGGTGACCGCCGACCAGCTGATCa                                                  <  1:181314/149‑1 (MQ=255)
                                                                                              cTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGCCCTGGGCCACGGTGACCGCCGACCAGCTGATCa                                                  <  1:309244/149‑1 (MQ=255)
                                                                                                       ggggAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGCCCTGGGCCACGGTGACCGCCGACCAGCTGATCAGCGACACGc                                         <  6:130006/149‑1 (MQ=255)
                                                                                                       ggggAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGCCCTGGGCCACGGTGACCGCCGACCAGCTGATCAGCGACACGc                                         <  4:169689/149‑1 (MQ=255)
                                                                                                       ggggAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGCCCTGGGCCACGGTGACCGCCGACCAGCTGATCAGCGACACGc                                         >  7:17576/1‑149 (MQ=255)
                                                                                                             ccATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCAGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGTCCTGGGCCACGGTGACCGCCGACCAGCTGATCAGCGACACGCCATAgg                                   >  7:312182/1‑149 (MQ=255)
                                                                                                             ccATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCAGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGTCCTGGGCCACGGTGACCGCCGACCAGCTGATCAGCGACACGCCATAgg                                   >  7:312178/1‑149 (MQ=255)
                                                                                                                   cGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGCCCTGGGCCACGGTGACCGCCGACCAGCTGATCAGCGACACGCCATAGGTCAGGc                             >  8:54098/1‑149 (MQ=255)
                                                                                                                                cTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGCCCTGGGCCACGGTGACCGCCGACCAGCTGATCAGCGACACGCCATAGGTCAGGCGTGAGCCGAAGCg                <  4:304001/149‑1 (MQ=255)
                                                                                                                                  gTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGCCCTGGGCCACGGTGACCGCCGACCAGCTGATCAGCGACACGCCATAGGTCAGGCGTGAGCCGAAGCGGt              >  8:234141/1‑149 (MQ=255)
                                                                                                                                      gCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGCCCTGGGCCACGGTGACCGCCGACCAGCTGATCAGCGACACGCCATAGGTCAGGCGTGAGCCGAAGCGGTCgg           >  8:29727/1‑148 (MQ=255)
                                                                                                                                       ccGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGCCCTGGGCCACGGTGACCGCCGACCAg                                                         <  3:16075/101‑1 (MQ=255)
                                                                                                                                       ccGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGCCCTGGGCCACGGTGACCGCCGACCAg                                                         >  4:16075/1‑101 (MQ=255)
                                                                                                                                       ccGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGCCCTGGGCCACGGTGACCGCCGACCAGCTGATCAGCGACACGCCATAGGTCAGGCGTGAGCCGAAGCGGTCGGc          >  8:306304/1‑148 (MQ=255)
                                                                                                                                         gggAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGCCCTGGGCCACGGTGACCGCCGACCAGCTGATCAGCGACACGCCATAGGTCAGGCGTGAGCCGAAGCGGTCGGCAAt       >  2:90422/1‑149 (MQ=255)
                                                                                                                                          ggAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGCCCTGGGCCACGGTGACCGCCGACCAGCTGATCAGCGACACGCCATAGGTCAGGCGTGAGCCGAAGCGGTCGGCAAt       >  2:164549/1‑148 (MQ=255)
                                                                                                                                            aaGGCCAGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGTCCTGGGCCACGGTGACCGCCGACCAGCTGATCAGCGACACGCCATAGGTCAGGCGTGAGCCGAAGCGGTCGGCAATGaa    <  8:312178/149‑1 (MQ=255)
                                                                                                                                            aaGGCCAGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGTCCTGGGCCACGGTGACCGCCGACCAGCTGATCAGCGACACGCCATAGGTCAGGCGTGAGCCGAAGCGGTCGGCAATGaa    <  8:312182/149‑1 (MQ=255)
                                                                                                                                              ggCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGCCCTGGGCCACGGTGACCGCCGACCAGCTGATCAGCGACACGCCATAGGTCAGGCGTGAGCCGAAGCGGTCgg           <  2:219570/140‑1 (MQ=255)
                                                                                                                                              ggCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGCCCTGGGCCACGGTGACCGCCGACCAGCTGATCAGCGACACGCCATAGGTCAGGCGTGAGCCGAAGCGGTCgg           >  1:219570/1‑140 (MQ=255)
                                                                                                                                                ccGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGCCCTGGGCCACGGTGACCGCCGACCAGCTGATCAGCGACACGCCATAGGTCAGGCGTGAGCCGAAGCGGTCGGCAATGAagc  >  1:179566/1‑147 (MQ=255)
                                                                                                                                                  ggggCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGCCCTGGGCCACGGTGACCGCCGACCAGCTGATCAGCGACACGCCATAGGTCAGGCGTGAGCCGAAGCGGTCGGCAATGAagc  >  2:264194/1‑145 (MQ=255)
                                                                                                                                                  ggggCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGCCCTGGGCCACGGTGACCGCCGACCAGCTGATCAGCGACACGCCATAGGTCAGGCGTGAGCCGAAGCGGTCGGCAATGAagc  <  1:264194/145‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
AGTCATAGGTGGTCGCCAGCCACAGCAGCGCACCGGTGAACAGCGCGGTGCCCAGGTACTGGCCACATACGTAGATGCTGCTGGCCATGCCGCGCTCGCGAGCGGGGAACCATACCGTAACCGCACGGCTGTTGGCCGGGAAGGCCGGGGCCTCCATGGCGCCGACCGCCAGGCGCAGGCCGAACAGTGAGGCGAAGCCGCTGGCCAGGCCCTGGGCCACGGTGACCGCCGACCAGCTGATCAGCGACACGCCATAGGTCAGGCGTGAGCCGAAGCGGTCGGCAATGAAGC  >  NC_002947/1910265‑1910555

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