Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,786,374 A→G 5.6% L293P (CTG→CCG)  mhpT ← 3‑(3‑hydroxy‑phenyl)propionate transporter MhpT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,786,3740AG 5.6% 251.7 / 2.2 82L293P (CTG→CCG) mhpT3‑(3‑hydroxy‑phenyl)propionate transporter MhpT
Reads supporting (aligned to +/- strand):  ref base A (37/40);  new base G (2/3);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGTTGCGCGCACCTGGGTCGAATACAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGCGGCCGGCTGAAGCCTTGGCCGAT  >  NC_002947/3786227‑3786519
                                                                                                                                                    |                                                                                                                                                 
cccGTTGCGCGCACCTGGGTCGAATACAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCGCGACAGCCCCAGCCCCGCCAACGACGCCGGc                                                                                                                                                 <  2:193710/150‑1 (MQ=255)
  cGTTGCGCGCACCTGGGTCGAATACAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCAt                                                                                                                                               >  1:918946/1‑150 (MQ=255)
    ttGCGCGCACCTGGGTCGAATACAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGc                                                                                                                                             <  2:1549866/150‑1 (MQ=255)
            aCCTGGGTCGAATACAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAg                                                                                                                                     <  1:923388/150‑1 (MQ=255)
             ccTGGGTCGAATACAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGGACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAgg                                                                                                                                    <  2:238538/150‑1 (MQ=255)
                      aaTACAGCGTCGCCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATa                                                                                                                           <  2:1478933/150‑1 (MQ=255)
                       aTACAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATAc                                                                                                                          <  1:1548244/150‑1 (MQ=255)
                          cAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATAccgc                                                                                                                       >  2:1612466/1‑150 (MQ=255)
                          cAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATAccgc                                                                                                                       <  1:1724187/150‑1 (MQ=255)
                           aGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATAccgcc                                                                                                                      >  2:1529130/1‑150 (MQ=255)
                                        gcgcGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAAcc                                                                                                         >  1:1409212/1‑150 (MQ=255)
                                           cgTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCcgc                                                                                                      <  2:154623/150‑1 (MQ=255)
                                              aCAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGt                                                                                                   >  1:1506488/1‑150 (MQ=255)
                                                 gCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCca                                                                                                <  2:1696815/150‑1 (MQ=255)
                                                  cACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCcat                                                                                               >  2:108385/1‑150 (MQ=255)
                                                  cACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCcat                                                                                               >  2:1888452/1‑150 (MQ=255)
                                                   aCCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCcatc                                                                                              <  2:1759069/150‑1 (MQ=255)
                                                    ccAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCcatca                                                                                             <  1:199157/150‑1 (MQ=255)
                                                    ccAACGGGCCACCAGTGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCcatca                                                                                             >  2:1818748/1‑150 (MQ=255)
                                                      accTGGCCACCAATTGCGCAAAACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAACCCCCCCAACGACGACAGCATGCCGGCGTAGGCGATCAATACCGCCCGCGCGGCAAACCCGCGCACcatcaaa                                                                                           >  2:156003/3‑148 (MQ=255)
                                                      aaCTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCatcatc                                                                                           <  2:1597708/150‑1 (MQ=255)
                                                       aCTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATcc                                                                                          <  1:108385/150‑1 (MQ=255)
                                                         tGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCTGTCCATCATCCGa                                                                                        >  1:390527/1‑150 (MQ=255)
                                                         tGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGa                                                                                        >  1:158695/1‑150 (MQ=255)
                                                         tGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGa                                                                                        >  2:711517/1‑150 (MQ=255)
                                                           gCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGAATGc                                                                                                                                             <  1:1014029/95‑1 (MQ=255)
                                                           gCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGAATGc                                                                                                                                             >  2:1014029/1‑95 (MQ=255)
                                                             caccaATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCgg                                                                                    >  2:925056/1‑150 (MQ=255)
                                                              accaATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCGGCATGCCGGCGTAGGCGa                                                                                                                                 >  2:1605445/1‑104 (MQ=255)
                                                              accaATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCGGCATGCCGGCGTAGGCGa                                                                                                                                 <  1:1605445/104‑1 (MQ=255)
                                                              accaATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGt                                                                                   >  2:1595339/1‑150 (MQ=255)
                                                              accaATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGt                                                                                   >  2:349564/1‑150 (MQ=255)
                                                                caATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATc                                                                                                                               <  2:984365/104‑1 (MQ=255)
                                                                caATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATc                                                                                                                               >  1:984365/1‑104 (MQ=255)
                                                                 aaTGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCTATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAg                                                                                >  2:214204/1‑150 (MQ=255)
                                                                  aTGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAgg                                                                               >  2:155990/1‑150 (MQ=255)
                                                                    ggCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGaa                                                                             <  2:1282676/150‑1 (MQ=255)
                                                                    ggCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGaa                                                                             <  1:1119050/150‑1 (MQ=255)
                                                                     gcgcAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAAc                                                                            >  1:701187/1‑150 (MQ=255)
                                                                       gcAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGa                                                                          <  2:476876/150‑1 (MQ=255)
                                                                         aaTACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATc                                                                        <  2:806382/150‑1 (MQ=255)
                                                                           tACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCa                                                                      <  1:1888452/150‑1 (MQ=255)
                                                                                tGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATcc                                                                       >  2:326569/1‑144 (MQ=255)
                                                                                tGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATcc                                                                       <  1:326569/144‑1 (MQ=255)
                                                                                        aCCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCcccaccca                                                         <  1:155990/150‑1 (MQ=255)
                                                                                         cccGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACACGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAg                                                        <  2:512072/150‑1 (MQ=255)
                                                                                          ccGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAgg                                                       >  2:1140551/1‑150 (MQ=255)
                                                                                           cGGCCAGCAACATCAGCCCGAACCTAGTTCG‑CAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGt                                                      <  1:941231/150‑1 (MQ=255)
                                                                                            ggCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGtt                                                     <  2:158695/150‑1 (MQ=255)
                                                                                             gCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTg                                                    <  1:711517/150‑1 (MQ=255)
                                                                                              ccAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGa                                                   <  1:349564/150‑1 (MQ=255)
                                                                                               cAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGGGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGaa                                                  <  1:470905/150‑1 (MQ=255)
                                                                                                aGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAAc                                                 <  2:1266083/150‑1 (MQ=255)
                                                                                                 gCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGCCCCCGTGATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACa                                                <  1:1056948/150‑1 (MQ=255)
                                                                                                   aaCATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGa                                              >  1:1830177/1‑150 (MQ=255)
                                                                                                    aCATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCAt                                                                                                                                               <  2:1198057/52‑1 (MQ=255)
                                                                                                    aCATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCAt                                                                                                                                               >  1:1198057/1‑52 (MQ=255)
                                                                                                      aTCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGaa                                                                             >  1:1896078/1‑116 (MQ=255)
                                                                                                      aTCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGaa                                                                             <  2:1896078/116‑1 (MQ=255)
                                                                                                                   t‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGcgcccg                              >  2:161656/1‑150 (MQ=255)
                                                                                                                      cTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGcgcccgcc                            >  1:1372493/1‑150 (MQ=255)
                                                                                                                       tCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGt                                                      <  1:1738116/123‑1 (MQ=255)
                                                                                                                       tCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGt                                                      >  2:1738116/1‑123 (MQ=255)
                                                                                                                       tCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCt                           >  1:1820978/1‑150 (MQ=255)
                                                                                                                          ccagccccagcccCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGcg                        <  2:1449271/149‑1 (MQ=255)
                                                                                                                          aCAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAAt                                                                                                                            >  2:1769719/1‑50 (MQ=255)
                                                                                                                          aCAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAAt                                                                                                                            <  1:1769719/50‑1 (MQ=255)
                                                                                                                            agccccagccccGCCAACGACGCCAGCATGCCGGCGTAGGCGGTCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCACCCCACCCAGGTTGAACAGAATCTGCAc                                      <  2:878156/134‑1 (MQ=255)
                                                                                                                            agccccagccccGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCACCCCACCCAGGTTGAACAGAATCTGCAc                                      >  1:878156/1‑134 (MQ=255)
                                                                                                                                ccagccccGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGcggccggc                  >  1:1323235/1‑150 (MQ=255)
                                                                                                                                    ccccGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGCGGCCGGCTGaa              <  1:1760172/150‑1 (MQ=255)
                                                                                                                                     cccGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGc                                        <  2:319526/123‑1 (MQ=255)
                                                                                                                                     cccGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGc                                        >  1:319526/1‑123 (MQ=255)
                                                                                                                                         ccAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGCGGCCGGCTGAAGCCtt         <  1:925056/150‑1 (MQ=255)
                                                                                                                                          cAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGGATCTGCACGGCGCCCGCCTGCGGCCGGCTGAAGCCTTg        <  2:28313/150‑1 (MQ=255)
                                                                                                                                           aaCGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGCGGCCGGCTGAAGCCTTgg       >  2:34612/1‑150 (MQ=255)
                                                                                                                                           aaCGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGCGGCCGGCTGAAGCCTTgg       >  2:1431436/1‑150 (MQ=255)
                                                                                                                                            acgacgCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGcggcc                     >  2:248144/1‑135 (MQ=255)
                                                                                                                                            acgacgCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGcggcc                     <  1:248144/135‑1 (MQ=255)
                                                                                                                                              tacgCCGGCATGCTGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGCGGCCGGCTGAAGCCTTGGCCg    >  2:1520062/2‑150 (MQ=255)
                                                                                                                                              gacgCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGCGGCCGGCTGAAGCCTTGGCCg    >  2:334908/1‑150 (MQ=255)
                                                                                                                                                cgCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCATGCGGCCGGCTGAAGCCTTGGCCGAt  <  2:1409212/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CCCGTTGCGCGCACCTGGGTCGAATACAGCGTCGGCGCCAGCGCGTACAGCACCAACTGGCCACCAATGGCGCAATACCCTGCAGTGAACCCGGCCAGCAACATCAGCCCGAAACT‑GCTCGACAGCCCCAGCCCCGCCAACGACGCCAGCATGCCGGCGTAGGCGATCAATACCGCCCGCCCGGCAAACCCGCGGTCCATCATCCGACCGGTCAGGAACGATCCAGCCGCCCCACCCAGGTTGAACAGAATCTGCACGGCGCCCGCCTGCGGCCGGCTGAAGCCTTGGCCGAT  >  NC_002947/3786227‑3786519

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