Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,458,795 A→G 5.9% intergenic (+43/‑225) pcaJ → / → trkH‑I 3‑oxoadipate CoA‑transferase subunit A/potassium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,458,7950AG 5.9% 203.2 / 2.3 70intergenic (+43/‑225)pcaJ/trkH‑I3‑oxoadipate CoA‑transferase subunit A/potassium transporter
Reads supporting (aligned to +/- strand):  ref base A (34/31);  new base G (2/2);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.90e-01

GGCCGTACTGGAAGTGACACCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGATAGCGGCGTAAAAAATTGCCGGCCAC  >  NC_002947/4458647‑4458944
                                                                                                                                                    |                                                                                                                                                     
ggCCGTACTGGAAGTGACACCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCAcaatcaa                                                                                                                                                      <  1:782528/150‑1 (MQ=255)
          gAAGTGACACCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGca                                                                                                                                            >  1:1699344/1‑150 (MQ=255)
                acacCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACTATTGATTGTGCCGca                                                                                                                                            <  2:667666/144‑1 (MQ=255)
                acacCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACTATTGATTGTGCCGca                                                                                                                                            >  1:667666/1‑144 (MQ=255)
                  acCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGcgcc                                                                                                                                    >  1:1169910/1‑150 (MQ=255)
                      gAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGa                                                                                                                                <  2:1169910/150‑1 (MQ=255)
                           gCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCAGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAgg                                                                                                                           >  1:1995745/1‑150 (MQ=255)
                            cTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAggg                                                                                                                          <  2:1926847/150‑1 (MQ=255)
                             tGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGa                                                                                                                         >  1:2297851/1‑150 (MQ=255)
                             tGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATAAATAGTGCCGCACAAGCGCCGCGAGCAGGGa                                                                                                                         <  1:2178346/150‑1 (MQ=255)
                               aaaGTGGTCGAAATCTGCGCGGACATCGACTTTGAAGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGc                                                                                                                       <  2:1219639/150‑1 (MQ=255)
                                 aGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCAc                                                                                                                     <  1:1984011/150‑1 (MQ=255)
                                 aGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCAc                                                                                                                     >  2:362443/1‑150 (MQ=255)
                                     gTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCgg                                                                                                                 <  2:282787/150‑1 (MQ=255)
                                       cGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTg                                                                                                               >  1:2337885‑M2/1‑149 (MQ=255)
                                        gAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCATATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTgc                                                                                                              >  1:24261‑M2/1‑148 (MQ=255)
                                             cTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCaccgcgtgctccctgctcgcggcgcttgtg                                                                                                         <  1:1782070‑M1/150‑31 (MQ=255)
                                              tGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGAGCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCAGACCCTTTAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTgaga                                                                                                        <  1:1009712/150‑1 (MQ=255)
                                                   ggACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTAcc                                                                                                   <  1:758629/150‑1 (MQ=255)
                                                     aCATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTAc                                                                                                    <  1:517566/147‑1 (MQ=255)
                                                     aCATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTAc                                                                                                    >  2:517566/1‑147 (MQ=255)
                                                       aTCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGa                                                                                                                         <  1:2051730/124‑1 (MQ=255)
                                                       aTCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGa                                                                                                                         >  2:2051730/1‑124 (MQ=255)
                                                          gACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGcc                                                                                                                                               >  1:1831160/1‑99 (MQ=255)
                                                          gACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGcc                                                                                                                                               <  2:1831160/99‑1 (MQ=255)
                                                               tGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCt                                                                                       >  2:728107/1‑150 (MQ=255)
                                                                 aCGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTc                                                                                     <  1:362443/150‑1 (MQ=255)
                                                                    aGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGaa                                                                                  <  2:2335067/150‑1 (MQ=255)
                                                                       tGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGa                                                                                                                         >  1:923866/1‑108 (MQ=255)
                                                                       tGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGa                                                                                                                         <  2:923866/108‑1 (MQ=255)
                                                                       tGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCACCGCGTGCTCCCTGCTCGCGGCgcttgtgcggcactattgattgtgccgcagagg                                                                               >  1:422658‑M2/1‑117 (MQ=255)
                                                                          aGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGAt                                                                            <  1:234141/150‑1 (MQ=255)
                                                                              aCTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAgg                                                                                                                           >  2:2253949/1‑99 (MQ=255)
                                                                              aCTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAgg                                                                                                                           <  1:2253949/99‑1 (MQ=255)
                                                                                tCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAAct                                                                      >  2:990534/1‑150 (MQ=255)
                                                                                 cAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACTATTGATTGTGCCGcagagg                                                                                                                                        >  2:2284936/1‑79 (MQ=25)
                                                                                 cAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACTATTGATTGTGCCGcagagg                                                                                                                                        <  1:2284936/83‑5 (MQ=25)
                                                                                  aGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATg                                                                           <  2:769517/143‑1 (MQ=255)
                                                                                  aGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATg                                                                           >  1:769517/1‑143 (MQ=255)
                                                                                  aGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAActct                                                                    >  2:2285040/1‑150 (MQ=255)
                                                                                   gTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACtctc                                                                   <  1:2285040/150‑1 (MQ=255)
                                                                                        gTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGgaga                                                              <  1:990534/150‑1 (MQ=255)
                                                                                        gTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGgaga                                                              <  1:728107/150‑1 (MQ=255)
                                                                                            cGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCCTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGGGAGGta                                                          >  1:2002387/1‑149 (MQ=255)
                                                                                             gCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCa                                                         >  1:661855/1‑150 (MQ=255)
                                                                                             gCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCAATGCATCTGCTGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCa                                                         >  2:1835651/1‑150 (MQ=255)
                                                                                              cTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGaa                                                                                  >  1:787077/1‑124 (MQ=255)
                                                                                              cTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGaa                                                                                  <  2:787077/124‑1 (MQ=255)
                                                                                               tGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCg                                                                              <  1:900895/127‑1 (MQ=255)
                                                                                               tGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCg                                                                              >  2:900895/1‑127 (MQ=255)
                                                                                                  tCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCa                                                    >  1:1787883/1‑150 (MQ=255)
                                                                                                      gTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGa                                                >  1:2094552/1‑150 (MQ=255)
                                                                                                         aTCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCt                                                                                                  <  2:2169358/97‑1 (MQ=255)
                                                                                                         aTCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCt                                                                                                  >  1:2169358/1‑97 (MQ=255)
                                                                                                          tCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCAGCTTCGAAGCCGATGGAACTCGCGGAGAGGTCAACCCAATGAGCTa                                            >  2:1424597/1‑150 (MQ=255)
                                                                                                              ggAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGc                                                                                >  2:502410/1‑110 (MQ=255)
                                                                                                              ggAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGc                                                                                <  1:502410/110‑1 (MQ=255)
                                                                                                               gAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAg                                       >  2:55011/1‑150 (MQ=255)
                                                                                                                aaGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGATTGAGACTACCTTTACAGCCGCTTCGAAGCAGATGCAACTATCGGAGAGGTCAACCCAATGAGCTATCGAGc                                      <  1:55011/150‑1 (MQ=255)
                                                                                                                 aGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCg                                     >  1:2309424/1‑150 (MQ=255)
                                                                                                                    cTTAGGGCTTGCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTGCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGa                                  <  1:883321/150‑1 (MQ=255)
                                                                                                                          gCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGa                            <  2:2094552/150‑1 (MQ=255)
                                                                                                                               ccATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCaccgcgtgctccctgctcgcggcgcttgtgcggca                                                                                                    >  1:899920‑M1/1‑38 (MQ=255)
                                                                                                                               ccATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCaccgcgtgctccctgctcgcggcgcttgtgcggca                                                                                                    <  2:899920‑M1/73‑36 (MQ=255)
                                                                                                                               ccATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCg                                     >  1:656356/1‑136 (MQ=255)
                                                                                                                               ccATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCg                                     <  2:656356/136‑1 (MQ=255)
                                                                                                                                 aTGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGATAgcggc                     <  1:1490691/150‑1 (MQ=255)
                                                                                                                                 aTGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTCGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGATAgcggc                     <  2:1699344/150‑1 (MQ=255)
                                                                                                                                    ccTCTGCGGCACAATCAATAGTGCCGCACAAGCgccgcgagcagggagcacgcggtgcttgtgcggcactattgattgtgccgcagaggcatggaaagccctaagtcttcctggatcacttgatcagcggcacgccactgagtttctgca                  >  1:422658‑M1/1‑33 (MQ=255)
                                                                                                                                      tctGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGATAGCGGCGTaaa                >  2:340049/1‑150 (MQ=255)
                                                                                                                                              acaATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGGGAGGTCAACCTAATGAGCTATCGAGCGGGAGCCAGATAGCGGCGTAAAAAATTGcc        >  1:1848484/1‑150 (MQ=255)
                                                                                                                                                    aatAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGATAGCGGCGTAAAAAATTGCCGGCCAc  >  1:1608019/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCCGTACTGGAAGTGACACCTGAAGGGCTGAAAGTGGTCGAAATCTGCGCGGACATCGACTTTGACGAGCTGCAGAAACTCAGTGGCGTGCCGCTGATCAAGTGATCCAGGAAGACTTAGGGCTTTCCATGCCTCTGCGGCACAATCAATAGTGCCGCACAAGCGCCGCGAGCAGGGAGCACGCGGTGCTTGAGACTACCTTTACAGCCGCTTCGAAGCCGATGGAACTCTCGGAGAGGTCAACCCAATGAGCTATCGAGCGGGAGCCAGATAGCGGCGTAAAAAATTGCCGGCCAC  >  NC_002947/4458647‑4458944

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