Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,909,797 C→T 5.7% V353M (GTG→ATG)  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,909,7970CT 5.7% 213.0 / 2.3 70V353M (GTG→ATG) PP_1710phthalate permease family MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (36/30);  new base T (2/2);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

CGGTGTGACGATGCCGCCGATGTTCGCCGCCAGGCTGGTAATGCCGCCGGTCAGGCCGATCAGCTCTTTTGGCGCCACTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCAAACGACAGCGAGTAGCCCTTGCGGATCAGCAGGTCGGAGACGATGCCCGCCAGGAG  >  NC_002947/1909657‑1909944
                                                                                                                                            |                                                                                                                                                   
cGGTGTGACGATGCCGCCGATGTTCGCCGCCAGGCTGGTAATGCCGCCGGTCAGGCCGATCAGCTCTTTTGGCGCCACTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATgca                                                                                                                                            <  1:1028082/150‑1 (MQ=255)
      gACGATGCCGCCGATGTTCGCCGCCAGGCTGGTAATGCCGCCGGTCAGGCCGATCAGCTCTTTTGGCGCCACTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATgcacagcac                                                                                                                                      >  2:1073569/1‑150 (MQ=255)
         gATGCCGCCGATGTTCGCCGCCAGGCTGGTCATGCCGCCGGTCAGGCCGATCAGCTCTTTTGGCGCCACTTCCGACACGGCAGCCCACGACGACGAGCCACGTCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACAtt                                                                                                                                   >  1:1573426/1‑150 (MQ=255)
           tGCCGCCGATGTTCGCCGCCAGGCTGGTAATGCCGCCGGTCAGGCCGATCAGCTCTTTTGGCGCCACTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTgg                                                                                                                                 >  2:1402542/1‑150 (MQ=255)
              cgccgATGTTCGCCGCCAGGCTGGTAATGCCGCCGGTCAGGCCGATCAGCTCTTTTGGCGCCACTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGt                                                                                                                              >  1:326179/1‑150 (MQ=255)
              cgccgATGTTCGCCGCCAGGCTGGTAATGCCGCCGGTCAGGCCGATCAGCTCTTTTGGCGCCACTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGt                                                                                                                              >  1:78450/1‑150 (MQ=255)
               gccgATGTTCGCCGCCAGGCTGGTAATGCCGCCGGTCAGGCCGATCAGCTCTTTTGGCGCCACTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTc                                                                                                                             >  1:655349/1‑150 (MQ=255)
                ccgATGTTCGCCGCCAGGCTGGTAATGCCGCCGGTCAGGCCGATCAGCTCTTTTGGCGCCACTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCg                                                                                                                            <  1:1402542/150‑1 (MQ=255)
                 cgATGTTCGCCGCCAGGCTGGTAATGCCGCCGGTCAGGCCGATCAGCTCTTTTGGCGCCACTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGAACAGCACATTGGAGTCgg                                                                                                                           <  2:1922886/150‑1 (MQ=255)
                   aTGTTCGCCGCCAGGCTGGTAATGCCGCCGGTCAGGCCGATCAGCTCTTTTGGCGCCACTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTg                                                                                                                         <  2:78450/150‑1 (MQ=255)
                    tGTTCGCCGCCAGGCTGGTAATGCCGCCGGTCAGGCCGATCAGCTCTTTTGGCGCCACTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGa                                                                                                                        <  1:1144119/150‑1 (MQ=255)
                     gTTCGCCGCCAGGCTGGTAATGCCGCCGGTCAGGCCGATCAGCTCTTTTGGCGCCACTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGaa                                                                                                                       <  1:1616674/150‑1 (MQ=255)
                              cAGGCTGGTAATGCCGCCGGTCAGGCCGATCAGCTCTTTTGGCGCCACTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAg                                                                                                              >  2:1913734/1‑150 (MQ=255)
                                     gTAATGCCGCCGGTTAGGCCGATCAGCTCTTTTGGCGCCACTTCCGACACGGCAGCCTACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATc                                                                                                       <  2:1808392/150‑1 (MQ=255)
                                     gTAATGCCGCCGGTCAGGCCGATCAGCTCTTTTGGCGCCACTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATc                                                                                                       <  1:1073569/150‑1 (MQ=255)
                                      tAATGCCGCCGGTCAGGCCGATCAGCTCTTTTGGCGCCACTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTg                                                                                                                                  <  2:239319/122‑1 (MQ=255)
                                      tAATGCCGCCGGTCAGGCCGATCAGCTCTTTTGGCGCCACTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTg                                                                                                                                  >  1:239319/1‑122 (MQ=255)
                                             gccgGTCAGGCCGATCAGCTCTTTTGGCGCCCCAGGCGACGCGGCAGCCCACGACGACGAGGCTATGCCCTGGGAAAAAAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATccc                                                                                               >  1:1911086/1‑150 (MQ=255)
                                             gccgGTCAGGCCGATCAGCTCTTTTGGCGCCACTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATccc                                                                                               >  2:680527/1‑150 (MQ=255)
                                                    aGGCCGATCAGCTCTTTTGGCGCCACTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATgcacag                                                                                                                                         <  1:581741/101‑1 (MQ=255)
                                                    aGGCCGATCAGCTCTTTTGGCGCCACTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATgcacag                                                                                                                                         >  2:581741/1‑101 (MQ=255)
                                                       ccGATCAGCTCTTTTGGCGCCACTTCTTACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAg                                                                                     <  2:1416029/150‑1 (MQ=255)
                                                       ccGATCAGCTCTTTTGGCGCCACTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAg                                                                                     <  2:655349/150‑1 (MQ=255)
                                                          aTCAGCTCTTTTGGCGCCACTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGccc                                                                                  >  2:721543/1‑150 (MQ=255)
                                                                tctTTTGGCGCCACTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCa                                                                                             <  2:1060388/133‑1 (MQ=255)
                                                                tctTTTGGCGCCACTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCa                                                                                             >  1:1060388/1‑133 (MQ=255)
                                                                       gcgcCACTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGGCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCa                                                                     >  2:1091397/1‑150 (MQ=255)
                                                                        cgcCACTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAg                                                                    >  1:742437/1‑150 (MQ=255)
                                                                             cTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATc                                                                                                       >  1:1466178/1‑110 (MQ=255)
                                                                             cTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATc                                                                                                       <  2:1466178/110‑1 (MQ=255)
                                                                                  gACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCa                                                          >  1:936443/1‑150 (MQ=255)
                                                                                       ggCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGATCGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCAAACGa                                                     <  1:1513843/150‑1 (MQ=255)
                                                                                        gCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCAAACGAc                                                    <  1:721543/150‑1 (MQ=255)
                                                                                         cAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCg                                                              <  1:187885/139‑1 (MQ=255)
                                                                                         cAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCg                                                              >  2:187885/1‑139 (MQ=255)
                                                                                            cccACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCAAACGACAGCg                                                >  1:1310960/1‑150 (MQ=255)
                                                                                             ccACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATc                                                                                        >  1:1535491/1‑109 (MQ=255)
                                                                                             ccACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATc                                                                                        >  1:1535705/1‑109 (MQ=255)
                                                                                             ccACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATc                                                                                        <  2:1535705/109‑1 (MQ=255)
                                                                                             ccACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATc                                                                                        <  2:1535491/109‑1 (MQ=255)
                                                                                             ccACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCAAACGACAGCGa                                               >  1:1553638/1‑150 (MQ=255)
                                                                                                    gacgaGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCCGGCAAACGACAGCGAGGAAccc                                        >  2:1079744/1‑150 (MQ=255)
                                                                                                      cgaGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCAAACGACAGCGAGTAGCCCtt                                      >  2:619776/1‑150 (MQ=255)
                                                                                                       gaGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCATGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCAAACGACAGCGAg                                              <  1:24213/141‑1 (MQ=255)
                                                                                                       gaGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCATGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCAAACGACAGCGAg                                              >  2:24213/1‑141 (MQ=255)
                                                                                                       gaGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTGCACCAGCAAGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCAAACGACAGCGAGTAGCCCTTg                                     >  1:315393/1‑150 (MQ=255)
                                                                                                         ggCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCAAACGACAGCGAGTAGCCCTTGCg                                   <  2:1312907/150‑1 (MQ=255)
                                                                                                         ggCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCAAACGACAGCGAGTAGCCCTTGCg                                   >  2:1917410/1‑150 (MQ=255)
                                                                                                           cTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCAc                                                                                                                 <  1:445373/70‑1 (MQ=255)
                                                                                                           cTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCAc                                                                                                                 >  2:445373/1‑70 (MQ=255)
                                                                                                            tATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCAAACGACAGCGAGTAGCCCTTGCGGAt                                <  1:1917410/150‑1 (MQ=255)
                                                                                                              tGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACCGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCAAACGACAGCGAGTAGCCCTTGCGGATca                              <  2:1325019/150‑1 (MQ=255)
                                                                                                                cccTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGTTGAAGTTCACCAGCACGATCGACAt                                                                                                  >  2:821780/1‑80 (MQ=255)
                                                                                                                cccTGCGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGTTGAAGTTCACCAGCACGATCGACAt                                                                                                  <  1:821780/80‑1 (MQ=255)
                                                                                                                 ccTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCAAACGACAGCGAGTAGCCCTTGCGGATcagca                           >  2:1311659/1‑150 (MQ=255)
                                                                                                                 ccTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCAAACGACAGCGAGTAGCCCTTGCGGATcagca                           >  1:673223/1‑150 (MQ=255)
                                                                                                                       cAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTg                                                                                                                                  >  1:1868193/1‑41 (MQ=255)
                                                                                                                       cAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTg                                                                                                                                  <  2:1868193/41‑1 (MQ=255)
                                                                                                                        aaaGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCAAACGACAGCGAGTAGCCCTTGCGGATCAGCAGGTCGga                    <  2:1553638/150‑1 (MQ=255)
                                                                                                                          aGAACGCCAGGGTCAGCATGGCGATGCACAGCACATTGGAGTCGGTGAAGt                                                                                                                     <  2:293420/51‑1 (MQ=255)
                                                                                                                          aGAACGCCAGGGTCAGCATGGCGATGCACAGCACATTGGAGTCGGTGAAGt                                                                                                                     >  1:293420/1‑51 (MQ=255)
                                                                                                                              cGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCAAACGACAGCGAGTAGCCCTTGCGGATCAGCAGGTCGGAGACGAt              >  1:713005/1‑150 (MQ=255)
                                                                                                                               gCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCAAACGACAGCGAGTAGCCCTTGCGGATCAGCAGGTCGGAGACGATg             >  1:763412/1‑150 (MQ=255)
                                                                                                                                ccAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCAAACGACAGCGAGTAGCCCTTGCGGATCAGCAGGTCGGAGACGATGc            >  1:1248309/1‑150 (MQ=255)
                                                                                                                                 cAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCAAACGa                                                     >  2:871777/1‑108 (MQ=255)
                                                                                                                                 cAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCAAACGa                                                     <  1:871777/108‑1 (MQ=255)
                                                                                                                                         gCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGAACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCAAACGACAGCGAGTAGCCCTTGCGGATCAGCAGGTCGGAGACGATGCCCGCCagga   <  2:326179/150‑1 (MQ=255)
                                                                                                                                          cACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCAAACGACAGCGAGTAGCCCTTGCGGATCAGCAGGTCGGAGACg                <  2:1284388/136‑1 (MQ=255)
                                                                                                                                          cACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCAAACGACAGCGAGTAGCCCTTGCGGATCAGCAGGTCGGAGACg                >  1:1284388/1‑136 (MQ=255)
                                                                                                                                          cACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCAAACGACAGCGAGTAGCCCTTGCGGATCAGCAGGTCGGAGACGATGCCCGCCaggag  <  1:1913734/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
CGGTGTGACGATGCCGCCGATGTTCGCCGCCAGGCTGGTAATGCCGCCGGTCAGGCCGATCAGCTCTTTTGGCGCCACTTCCGACACGGCAGCCCACGACGACGAGGCTATGCCCTGGGCAAAGAACGCCAGGGTCAGCACGGCGATGCACAGCACATTGGAGTCGGTGAAGTTCACCAGCACGATCGACATCCCCAGCATCGAGCCCACCACCAGCGGCAGCTTGCGGGCAAACGACAGCGAGTAGCCCTTGCGGATCAGCAGGTCGGAGACGATGCCCGCCAGGAG  >  NC_002947/1909657‑1909944

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