Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,155,334 C→T 7.9% G68G (GGG→GGA) ‡ PP_3658 ← aromatic compound MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,155,3340CT 7.9% 255.1 / 9.9 89G68G (GGG→GGA) ‡PP_3658aromatic compound MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (40/42);  new base T (3/4);  total (43/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

AGCGCATGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATAAGTCC  >  NC_002947/4155189‑4155480
                                                                                                                                                 |                                                                                                                                                  
aGCGCATGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGa                                                                                                                                                >  2:1781322/1‑150 (MQ=255)
 gcgcATGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGaa                                                                                                                                               >  1:1246861/1‑150 (MQ=255)
  cgcATGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAAc                                                                                                                                              >  1:1462377/1‑150 (MQ=255)
      tGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATg                                                                                                                                          <  1:76650/150‑1 (MQ=255)
          aCTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTTCGGGACAATGCTgg                                                                                                                                      <  2:1446283/150‑2 (MQ=255)
          aCTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTTCGGGACAATGCTg                                                                                                                                       <  2:70608/149‑1 (MQ=255)
          aCTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTTCGGGACAATGCTg                                                                                                                                       >  1:70608/1‑149 (MQ=255)
            tGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAACGACATAGCGCCCAGCATTGTCCCGAACAATGCTGagc                                                                                                                                    <  1:2268638/150‑1 (MQ=255)
            tGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGagc                                                                                                                                    <  1:38728/150‑1 (MQ=255)
            tGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGagc                                                                                                                                    >  2:2041322/1‑150 (MQ=255)
              aaGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGagcca                                                                                                                                  <  1:2112330/150‑1 (MQ=255)
               aGGTCAGCGTGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGagccag                                                                                                                                 >  2:408562/1‑150 (MQ=255)
                  tCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAgccagcca                                                                                                                              <  2:1246861/150‑1 (MQ=255)
                  tCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAgccagcca                                                                                                                              <  1:642870/150‑1 (MQ=255)
                    aGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCAtt                                                                                                                            >  1:1841137/1‑150 (MQ=255)
                        ggCTGTCAGTGAACCCTGACGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCTCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTAccc                                                                                                                        <  2:1165588/150‑1 (MQ=255)
                             tCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGgcg                                                                                                                   >  1:1984371/1‑150 (MQ=255)
                                gTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTg                                                                                                                >  1:684882/1‑150 (MQ=255)
                                     cccTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGc                                                                                                           >  2:330312/1‑150 (MQ=255)
                                        tGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTc                                                                                                        >  1:367377/1‑150 (MQ=255)
                                               gCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGc                                                                                                 >  2:426044/1‑150 (MQ=255)
                                               gCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATACCGCCTAGCATTGCCCCGACCAATGCTGCGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGc                                                                                                 <  1:1548274/150‑1 (MQ=255)
                                                     aaaaCGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCTCAGCATTGTTCGGGACAATGCTggg                                                                                                                                     >  2:2334640/1‑106 (MQ=38)
                                                     aaaaCGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCTCAGCATTGTTCGGGACAATGCTggg                                                                                                                                     <  1:2334640/108‑3 (MQ=38)
                                                      aaaCGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGTGTGCGTCAATTTCCATATCgg                                                                                          >  2:559295/1‑150 (MQ=255)
                                                         cGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCa                                                                                       <  1:408562/150‑1 (MQ=255)
                                                            tGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCatga                                                                                    >  1:1284773/1‑150 (MQ=255)
                                                             gAACAGCTCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCatgat                                                                                   <  1:1611080/150‑1 (MQ=255)
                                                                 aGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGAtggctg                                                                              >  2:1561964/1‑150 (MQ=255)
                                                                  gcgcTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGAtggctg                                                                              <  2:1984371/150‑1 (MQ=255)
                                                                   cgcTACACAGATCGCGATCATTTTCGGCAGGCAACCCCCATTTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGAtggctgg                                                                             <  2:1398291/150‑1 (MQ=255)
                                                                    gcTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGagccagc                                                                                                                                >  2:152053/1‑98 (MQ=255)
                                                                    gcTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGagccagc                                                                                                                                <  1:152053/98‑1 (MQ=255)
                                                                      tACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTggcag                                                                          <  1:1781322/150‑1 (MQ=255)
                                                                        cacaGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTggcaggg                                                                        <  1:1189617/150‑1 (MQ=255)
                                                                         acaGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTggcagggc                                                                       >  1:2103678/1‑150 (MQ=255)
                                                                             aTCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGcc                                                                   >  2:1533266/1‑150 (MQ=255)
                                                                                  gATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCTCAGCATTGTTCGGGACAATGCTggg                                                                                                                                     <  1:2161026/79‑3 (MQ=25)
                                                                                  gATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCTCAGCATTGTTCGGGACAATGCTggg                                                                                                                                     >  2:2161026/1‑77 (MQ=25)
                                                                                  gATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCccac                                                              >  2:1143121/1‑150 (MQ=255)
                                                                                   aTCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCccacc                                                             >  2:882381/1‑150 (MQ=255)
                                                                                    tCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGgc                                                                                                                    <  1:646430/94‑1 (MQ=255)
                                                                                    tCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGgc                                                                                                                    >  2:646430/1‑94 (MQ=255)
                                                                                      aTTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCccacca                                                            <  2:1511956/148‑1 (MQ=255)
                                                                                      aTTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCccacca                                                            >  1:1511956/1‑148 (MQ=255)
                                                                                      aTTTTCGGCCGGCCAATTCGATATGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACt                                                          <  2:351714/150‑1 (MQ=255)
                                                                                             gccggccAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGt                                                   >  1:1686861/1‑150 (MQ=255)
                                                                                                 gccAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTa                                               >  2:2010918/1‑150 (MQ=255)
                                                                                                         gATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCatga                                                                                    >  1:56362/1‑105 (MQ=255)
                                                                                                         gATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCatga                                                                                    <  2:56362/105‑1 (MQ=255)
                                                                                                         gATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGa                                       >  2:2122188/1‑150 (MQ=255)
                                                                                                          aTCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTTCCAAGTCGTAGCCATCGAt                                      >  1:2035779/1‑150 (MQ=255)
                                                                                                          aTCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGAt                                      >  1:2220795/1‑150 (MQ=255)
                                                                                                            cTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATaa                                    <  2:367377/150‑1 (MQ=255)
                                                                                                              gCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTggcagggca                                                                      >  2:1179235/1‑114 (MQ=255)
                                                                                                              gCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTggcagggca                                                                      <  1:1179235/114‑1 (MQ=255)
                                                                                                               ccAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGa                                       >  2:1533316/1‑144 (MQ=255)
                                                                                                               ccAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGa                                       <  1:1533316/144‑1 (MQ=255)
                                                                                                                  aTGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGaa                              <  1:882381/150‑1 (MQ=255)
                                                                                                                   tgtgCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAAt                             >  2:659774/1‑150 (MQ=255)
                                                                                                                   tgtgCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAAt                             <  1:2010918/150‑1 (MQ=255)
                                                                                                                   ggtgCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCCCCACTGCCAAGTCGTAGCCATCGATAATCAGAAt                             <  2:463301/149‑1 (MQ=255)
                                                                                                                      gCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCaa                          <  2:2220795/150‑1 (MQ=255)
                                                                                                                      gCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCaa                          <  1:426044/150‑1 (MQ=255)
                                                                                                                      gCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCaa                          <  2:600497/150‑1 (MQ=255)
                                                                                                                       ccTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCaa                          >  1:2229973/1‑149 (MQ=255)
                                                                                                                       ccTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCaa                          <  2:2229973/149‑1 (MQ=255)
                                                                                                                         tAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCATCACTGCCACGTCGTAGCCATCGATAATCAGAACCAacac                       <  2:1841137/150‑1 (MQ=255)
                                                                                                                         tAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCaa                                                     <  1:815299/120‑1 (MQ=255)
                                                                                                                         tAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCaa                                                     >  2:815299/1‑120 (MQ=255)
                                                                                                                         tAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGaa                              <  1:365187/143‑1 (MQ=255)
                                                                                                                         tAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGaa                              >  2:365187/1‑143 (MQ=255)
                                                                                                                         tAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAacac                       >  2:2357581/1‑150 (MQ=255)
                                                                                                                         tAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAacac                       <  1:1561964/150‑1 (MQ=255)
                                                                                                                         tAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAacac                       >  2:169260/1‑150 (MQ=255)
                                                                                                                                  aGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCagaa              <  1:1533266/150‑1 (MQ=255)
                                                                                                                                    cgcCCCGGATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATACCGGCCATGATGGCCGGCAGGGCAGCCCCCATCCCTGCGAAGTCGTAGCTATCGATAATCAGAATCAACACGCACCagaata            <  1:559295/150‑1 (MQ=255)
                                                                                                                                     gcCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCACTCAGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCagaatag           >  2:329072/1‑150 (MQ=255)
                                                                                                                                     gcCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCagaatag           <  2:2103678/150‑1 (MQ=255)
                                                                                                                                       ccAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCagaatagaa         <  2:206301/150‑1 (MQ=255)
                                                                                                                                       ccAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCagaatagaa         >  1:206301/1‑150 (MQ=255)
                                                                                                                                        cAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCagaatagaat        >  1:280857/1‑150 (MQ=255)
                                                                                                                                         aGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGt                                                   >  2:584761/1‑106 (MQ=255)
                                                                                                                                         aGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGt                                                   <  1:584761/106‑1 (MQ=255)
                                                                                                                                          gCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATaa      <  2:2017345/150‑1 (MQ=255)
                                                                                                                                          gCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATaa      <  2:1462377/150‑1 (MQ=255)
                                                                                                                                             ttGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGtt                                                                                                             <  1:16570/44‑1 (MQ=255)
                                                                                                                                             ttGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGtt                                                                                                             >  2:16570/1‑44 (MQ=255)
                                                                                                                                              tGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATAAGTcc  <  1:621621/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AGCGCATGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATAAGTCC  >  NC_002947/4155189‑4155480

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