Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,926,267 C→T 6.3% A90T (GCC→ACC)  peaA ← quinohemoprotein amine dehydrogenase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,926,2670CT 6.3% 380.4 / 9.4 127A90T (GCC→ACC) peaAquinohemoprotein amine dehydrogenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (61/57);  new base T (4/4);  total (66/61)
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.62e-01

Only 100 of 129 total aligned reads displayed.

GTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGCTCGGCCGCCATGGCCTGGGTGCTCAGCAGCG  >  NC_002947/3926118‑3926416
                                                                                                                                                      |                                                                                                                                                     
gTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGc                                                                                                                                                       <  2:1795167/150‑1 (MQ=255)
   cTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGAt                                                                                                                                                    >  1:1134025/1‑150 (MQ=255)
    tCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATg                                                                                                                                                   <  1:2843005/150‑1 (MQ=255)
      gACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGGGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCGCCCGCATGCAGGTGAt                                                                                                                                                    <  1:312331/147‑1 (MQ=255)
      gACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCGCCCGCATGCAGGTGAt                                                                                                                                                    >  2:137946/1‑147 (MQ=255)
      gACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCGCCCGCATGCAGGTGAt                                                                                                                                                    <  1:137946/147‑1 (MQ=255)
              tGAGCCGAAGCTCCATGGCGTTGCGAACGCACTGGGTTTCGCAGGCCGCGAGTATCGGTTGGCCGGCTAGGGATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATcagcc                                                                                                                                         <  2:1439776/150‑1 (MQ=255)
                 gCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATcagccagc                                                                                                                                      >  2:1552813/1‑150 (MQ=255)
                 gCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATcagccagc                                                                                                                                      <  2:1590561/150‑1 (MQ=255)
                          ccATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGgtg                                                                                                                             <  1:3046852/150‑1 (MQ=255)
                           gaTGGCGTGCGCTACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGgtgt                                                                                                                            >  1:1164223/2‑150 (MQ=255)
                           cATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGgtgt                                                                                                                            >  1:1969173/1‑150 (MQ=255)
                            aTGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTc                                                                                                                           <  2:1277542/150‑1 (MQ=255)
                            aTGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTc                                                                                                                           <  1:1636890/150‑1 (MQ=255)
                             tGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCt                                                                                                                          >  1:399733/1‑150 (MQ=255)
                                    gcgcACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCtg                                                                                                                   <  2:25321/150‑1 (MQ=255)
                                      gcACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCCGCCAGCCTTCCGGTGTCTTTCGCtggt                                                                                                                 <  2:1134025/150‑1 (MQ=255)
                                      gcACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCtggt                                                                                                                 <  2:1397173/150‑1 (MQ=255)
                                       cACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCtggtg                                                                                                                >  2:1285767/1‑150 (MQ=255)
                                       cACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCACATCAGCCAGCCTTCCGGTGTCTTGCGCtggtg                                                                                                                <  2:2553644/150‑1 (MQ=255)
                                           gCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACGAGCGTGCGACGATCGTCGTTGCTGATCTGAAGGCCATGCATCACCTGCATGCGAGCGATGCGCATCAGCCCGCCTTCCAGTGTCTTGCGCTGGTGGCt                                                                                                             >  2:232333/1‑150 (MQ=255)
                                              aTCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCgg                                                                                                                               >  2:2127120/1‑129 (MQ=255)
                                              aTCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCgg                                                                                                                               <  1:2127120/129‑1 (MQ=255)
                                               tCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATcagcca                                                                                                                                        <  2:719564/119‑1 (MQ=255)
                                               tCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATcagcca                                                                                                                                        >  1:719564/1‑119 (MQ=255)
                                                     tCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCta                                                                                                   <  1:2516417/150‑1 (MQ=255)
                                                     tCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCta                                                                                                   <  2:2393271/150‑1 (MQ=255)
                                                        cTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAg                                                                                                <  1:1325758/150‑1 (MQ=255)
                                                         tGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAgg                                                                                               >  2:2270500/1‑150 (MQ=255)
                                                         tGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAgg                                                                                               >  1:1287819/1‑150 (MQ=255)
                                                         tGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAgg                                                                                               >  1:1484357/1‑150 (MQ=255)
                                                          gggCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGgt                                                                                              <  1:1401830/150‑1 (MQ=255)
                                                           ggCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATTCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCtt                                                                                                                         >  1:1750727/1‑122 (MQ=255)
                                                            gcgcCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGgtgt                                                                                            <  2:1751356/150‑1 (MQ=255)
                                                              gcCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCAACACCTGCATGCGGGAGATGCTAATCAGCCCGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCg                                                                                          >  1:985772/1‑150 (MQ=255)
                                                                cAGGCCCTGTTTGTCGGCCAGGTATTTGACAAGCGTGCGACGATAGTCGTCGCTGATCTGACGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGtt                                                                                        >  2:1122331/1‑150 (MQ=255)
                                                                 aGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTg                                                                                       <  2:399733/150‑1 (MQ=255)
                                                                 aGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTg                                                                                       >  1:777288/1‑150 (MQ=255)
                                                                  ggCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGc                                                                                      <  1:1598871/150‑1 (MQ=255)
                                                                   gCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGcc                                                                                     >  1:2251153/1‑150 (MQ=255)
                                                                   gCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGcc                                                                                     >  1:697550/1‑150 (MQ=255)
                                                                        gTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCgg                                                                                <  2:2361745/150‑1 (MQ=255)
                                                                         tttGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGtt                                                                                        <  2:483476/141‑1 (MQ=255)
                                                                         tttGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGtt                                                                                        >  1:483476/1‑141 (MQ=255)
                                                                           tGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCGTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGa                                                                             <  2:1986999/150‑1 (MQ=255)
                                                                            gTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGat                                                                            >  1:1294381/1‑150 (MQ=255)
                                                                              cGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGatat                                                                          <  1:2127193/150‑1 (MQ=255)
                                                                                gCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATgg                                                                        <  2:3025556/150‑1 (MQ=255)
                                                                                gCCAGGTATTTGACCAGCGTGCGACGACCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCtgg                                                                                                                  >  2:1701223/1‑108 (MQ=255)
                                                                                gCCAGGTATTTGACCAGCGTGCGACGACCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCtgg                                                                                                                  <  1:1701223/108‑1 (MQ=255)
                                                                                 ccAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGTGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGgtg                                                                                             >  2:2414264/1‑128 (MQ=255)
                                                                                 ccAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGTGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGgtg                                                                                             <  1:2414264/128‑1 (MQ=255)
                                                                                           gACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCa                                                             >  1:2119384/1‑150 (MQ=255)
                                                                                            aCCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGatat                                                                          >  1:2945742/1‑136 (MQ=255)
                                                                                            aCCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGatat                                                                          <  2:2945742/136‑1 (MQ=255)
                                                                                            aCCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCAc                                                            >  1:2051092/1‑150 (MQ=255)
                                                                                            aCCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCAc                                                            <  1:2270500/150‑1 (MQ=255)
                                                                                              cAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACtt                                                          >  1:2998081/1‑150 (MQ=255)
                                                                                                gCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGt                                                        >  1:865078/1‑150 (MQ=255)
                                                                                                gCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGt                                                        >  2:3030140/1‑150 (MQ=255)
                                                                                                   tGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCt                                                     <  2:2051092/150‑1 (MQ=255)
                                                                                                      gacgaTCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCa                                                  >  2:2467767/1‑150 (MQ=255)
                                                                                                         gaTCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCg                                                                                                                                >  2:1404956/1‑69 (MQ=255)
                                                                                                         gaTCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCa                                                                      <  1:1355882/127‑1 (MQ=255)
                                                                                                         gaTCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCa                                                                      >  2:1355882/1‑127 (MQ=255)
                                                                                                         gaTCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACa                                               <  2:1294381/150‑1 (MQ=255)
                                                                                                         gaTCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACa                                               <  2:2998081/150‑1 (MQ=255)
                                                                                                                gtcgCTGATCTGCAGGCCATGCATCACCTGCATGCGGGCGCTGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATGGGTGTCGTTGTTTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCgg                                        >  2:1337726/1‑150 (MQ=255)
                                                                                                                gtcgCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCgg                                        >  2:2887010/1‑150 (MQ=255)
                                                                                                                   gCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGCTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGcc                                     >  1:32448/1‑150 (MQ=255)
                                                                                                                    cTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGccc                                    <  1:3030140/150‑1 (MQ=255)
                                                                                                                      gATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTg                                  >  1:2714449/1‑150 (MQ=255)
                                                                                                                         cTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGCTc                               >  2:2922127/1‑150 (MQ=255)
                                                                                                                          tGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGCTCg                              >  2:2225086/1‑150 (MQ=255)
                                                                                                                           gAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCgg                                                                                                      <  1:666400/77‑1 (MQ=255)
                                                                                                                            aaGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGCTCGgc                            >  1:2213705/1‑150 (MQ=255)
                                                                                                                              ggCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGCTCGgccg                          >  1:1553506/1‑150 (MQ=255)
                                                                                                                                ccATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTc                                                      <  1:1277241/120‑1 (MQ=255)
                                                                                                                                ccATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTc                                                      >  2:1277241/1‑120 (MQ=255)
                                                                                                                                atgagCATCGCCCGCATGCAGGTGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCgg                                                                                <  1:2461592/90‑1 (MQ=37)
                                                                                                                                 cATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAg                                              >  2:1469539/1‑127 (MQ=255)
                                                                                                                                 cATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAg                                              <  1:1469539/127‑1 (MQ=255)
                                                                                                                                 cATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGGCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGCTCggaagacc                       >  2:2174838/1‑144 (MQ=255)
                                                                                                                                  aTGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCgg                                                                                <  1:307533/92‑1 (MQ=255)
                                                                                                                                   tGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGCTCGGCCGCCATg                     <  1:1032539/150‑1 (MQ=255)
                                                                                                                                       tCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGCTCGGCCGCCATGGCCt                 >  1:2581240/1‑150 (MQ=255)
                                                                                                                                         aCCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGtt                                                       <  2:703743/110‑1 (MQ=255)
                                                                                                                                         aCCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGCTCGGCCGCCATGGCCTgg               <  2:2251153/150‑1 (MQ=255)
                                                                                                                                         aCCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGCTCGGCCGCCATGGCCTgg               <  2:304692/150‑1 (MQ=255)
                                                                                                                                         aCCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTCCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGCTCGGCCGCCATGGCCTgg               <  2:2213705/150‑1 (MQ=255)
                                                                                                                                          ccTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTc                                                      >  2:2313031/1‑110 (MQ=255)
                                                                                                                                           cTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTc                                                      >  1:1787020/1‑109 (MQ=255)
                                                                                                                                           cTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTc                                                      <  2:1787020/109‑1 (MQ=255)
                                                                                                                                           cTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGCTCGGCCGCCATGGCCTGGGt             >  2:400682/1‑150 (MQ=255)
                                                                                                                                                tGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGCTCGGCCGCCATGGCCTGGGTGCTca        <  2:1164223/150‑1 (MQ=255)
                                                                                                                                                tGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGCTCGGCCGCCATGGCCTGGGTGCTca        <  2:2581240/150‑1 (MQ=255)
                                                                                                                                                   gggCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGCTCGGCCGCCATGGCCTGGGTGCTcagca     <  1:1122331/150‑1 (MQ=255)
                                                                                                                                                   gggCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGCTCGGCCGCCATGGCCTGGGTGCTcagca     >  2:2086837/1‑150 (MQ=255)
                                                                                                                                                    ggCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGCTCGGCCGCCATGGCCTGGGTGCTcagcag    >  1:1929721/1‑150 (MQ=255)
                                                                                                                                                      cGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGCTCGGCCGCCATGGCCTGGGTGCGCagcagca  <  2:777288/150‑2 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGTGGCTGATGCGGCTATAGGTGTCGTTGCCTTCGGGGATATGGCACCCCATGCACTTGTTCTGCAACAGGCTCGGGCCCTGCTCGGCCGCCATGGCCTGGGTGCTCAGCAGCG  >  NC_002947/3926118‑3926416

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: