Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,486,224 C→G 6.9% A258G (GCC→GGC)  spuC‑I → polyamine:pyruvate transaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,486,2240CG 6.9% 225.0 / 8.0 72A258G (GCC→GGC) spuC‑Ipolyamine:pyruvate transaminase
Reads supporting (aligned to +/- strand):  ref base C (32/35);  new base G (3/2);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.73e-01

GAGAAGAAGATCCTCGAAGTCGGCGTGGACAACGTTGCCGCCTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATCCCCATGGGCGGGGTGATCGTGCGTGACG  >  NC_002947/2486076‑2486364
                                                                                                                                                    |                                                                                                                                            
gagaAGAAGATCCTCGAAGTCGGCGTGGACAACGTTGCCGCCTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGcc                                                                                                                                             >  2:812828/1‑150 (MQ=255)
   aagaagaTCCTCGAAGTCGGCGTGGACAACGTTGCCGCCTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCcgac                                                                                                                                          >  2:374757/1‑150 (MQ=255)
        gaTCCTCGAAGTCGGCGTGGACAACGTTGCCGCCTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGt                                                                                                                                     >  2:571755/1‑150 (MQ=255)
          tCCTCGAAGTCGGCGTGGACAACGTTGCCGCCTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGa                                                                                                                                   <  1:374757/150‑1 (MQ=255)
           ccTCGAAGTCGGCGTGGACAACGTTGCCGCCTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGAt                                                                                                                                  <  1:1462645/150‑1 (MQ=255)
                        gTGGACAACGTTGCCGCCTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGc                                                                                                                     <  1:446473/150‑1 (MQ=255)
                         tGGACAACGTTGCCGCCTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGcc                                                                                                                    <  1:590913/150‑1 (MQ=255)
                             cAACGTTGCCGCCTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTAc                                                                                                                >  2:1813767/1‑150 (MQ=255)
                              aaCGTTGCCGCCTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTAcc                                                                                                               >  1:1922480/1‑150 (MQ=255)
                                 gTTGCCGCCTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGt                                                                                                            <  2:278742/150‑1 (MQ=255)
                                  ttGCCGCCTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTg                                                                                                           >  2:1651805/1‑150 (MQ=255)
                                     ccgccTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGt                                                                                                        >  1:921942/1‑150 (MQ=255)
                                     ccgccTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGt                                                                                                        <  1:180611/150‑1 (MQ=255)
                                     ccgccTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGt                                                                                                        <  1:1651805/150‑1 (MQ=255)
                                       gccTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTgg                                                                                                      >  2:830156/1‑150 (MQ=255)
                                       gccTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATGTGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTgg                                                                                                      <  1:512760/150‑1 (MQ=255)
                                         cTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCgg                                                                                                             >  2:1649208/1‑141 (MQ=255)
                                         cTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCgg                                                                                                             <  1:1649208/141‑1 (MQ=255)
                                         cTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGt                                                                                                            <  1:950632/142‑1 (MQ=255)
                                         cTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGt                                                                                                            >  2:950632/1‑142 (MQ=255)
                                         cTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGtt                                                                                                    >  1:1546867/1‑150 (MQ=255)
                                             aTCGCCGGGCCGATCCAGGGCGCCGGTGGCGGGATCCTTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGTCCAGGTACGACATCATCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGt                                                                                                <  2:890030/150‑1 (MQ=255)
                                              tCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGtt                                                                                                                         <  1:1661784/124‑1 (MQ=255)
                                              tCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGtt                                                                                                                         >  2:1661784/1‑124 (MQ=255)
                                                  cgagccgaTCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGa                                                                                           <  2:1546867/150‑1 (MQ=255)
                                                    agccgaTCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGatt                                                                                         >  1:1159660/1‑150 (MQ=255)
                                                               ggCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTAGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTc                                                                              >  2:502745/1‑150 (MQ=255)
                                                                gcgcCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCa                                                                             >  2:1829743/1‑150 (MQ=255)
                                                                 cgcCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCaa                                                                            >  1:533924/1‑150 (MQ=255)
                                                                 cgcCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCaa                                                                            <  2:79748/150‑1 (MQ=255)
                                                                        ggCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGcctgac                                                                     <  1:1813767/150‑1 (MQ=255)
                                                                          cGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGcctgacct                                                                   >  2:437621/1‑150 (MQ=255)
                                                                           gTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGGCGACGAAgagg                                                                                                                                   <  1:1066300/85‑4 (MQ=255)
                                                                           gTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGGCGACGAAgagg                                                                                                                                   >  2:1066300/1‑82 (MQ=255)
                                                                           gTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTCCGGTACCGGTTATTACGACCTCAAGcctgacctg                                                                  >  2:1699709/1‑150 (MQ=255)
                                                                            tGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGcctgacctga                                                                 <  1:872820/150‑1 (MQ=255)
                                                                            tGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGcctg                                                                       <  1:1356548/144‑1 (MQ=255)
                                                                            tGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGcctg                                                                       >  2:1356548/1‑144 (MQ=255)
                                                                             gATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCtgat                                                                <  2:1159660/150‑1 (MQ=255)
                                                                             gATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCtgat                                                                <  1:812828/150‑1 (MQ=255)
                                                                               tCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCtgatga                                                              <  2:1922480/150‑1 (MQ=255)
                                                                                    ccgccgCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATc                                                         <  1:1829743/150‑1 (MQ=255)
                                                                                        cgcAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCa                                                     >  2:547645/1‑150 (MQ=255)
                                                                                          cAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAg                                                   <  1:437621/150‑1 (MQ=255)
                                                                                             aCCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGcct                                                                        >  2:1636472/1‑126 (MQ=255)
                                                                                             aCCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGcct                                                                        <  1:1636472/126‑1 (MQ=255)
                                                                                                 aTTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTCCGACCTCAAGCCTGACCTGAGTACCATCGCCAAGGGCCTGa                                            <  1:830156/150‑1 (MQ=255)
                                                                                                       cGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCg                                      <  2:533924/150‑1 (MQ=255)
                                                                                                                 ggAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTAGCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACAATCGCCCAGGGCCTGACCTCCGGTTGCATccc                            >  2:132778/1‑150 (MQ=255)
                                                                                                                    aaTCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGCCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTTACCTAATGACCATCGCCAAGCGCCTGACCTCCGGTTACATCCCCAt                         <  2:921424/150‑1 (MQ=255)
                                                                                                                     aTCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGc                                                                          >  1:1349700/1‑100 (MQ=255)
                                                                                                                     aTCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGc                                                                          <  2:1349700/100‑1 (MQ=255)
                                                                                                                      tCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTAcc                                                                                             >  2:1856945/1‑80 (MQ=255)
                                                                                                                      tCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTAcc                                                                                             <  1:1856945/80‑1 (MQ=255)
                                                                                                                      tCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATCCCCATgg                       >  2:696912/1‑150 (MQ=255)
                                                                                                                            ccAGGTACGACATCCTCTTAGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATCCCCATGGGCgggg                 <  1:502745/150‑1 (MQ=255)
                                                                                                                             cAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCAGCGCCAAGGGCCTGACCTCCGGTTACATCCCCATGGGCGGGGt                >  2:1603311/1‑150 (MQ=255)
                                                                                                                              aGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATCCCCATGGGCGGGTTg               <  1:571755/150‑1 (MQ=255)
                                                                                                                                gTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGAcc                                          >  2:1906976/1‑121 (MQ=255)
                                                                                                                                gTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGAcc                                          <  1:1906976/121‑1 (MQ=255)
                                                                                                                                gTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATCCCCATGGGCGGGGTGa              >  2:1622626/1‑149 (MQ=255)
                                                                                                                                gTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATCCCCATGGGCGGGGTGa              <  1:1622626/149‑1 (MQ=255)
                                                                                                                                   cGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGcctgacctga                                                                 <  1:786430/95‑1 (MQ=255)
                                                                                                                                   cGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGcctgacctga                                                                 >  2:786430/1‑95 (MQ=255)
                                                                                                                                       aTCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATCCCCATGGGCGGGGTGATcgtgc        <  1:1909602/148‑1 (MQ=255)
                                                                                                                                       aTCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATCCCCATGGGCGGGGTGATcgtgc        >  2:1909602/1‑148 (MQ=255)
                                                                                                                                         ccTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATCCCCATGGGCGGGGTGATCGTGCGTGa    >  1:204760/1‑150 (MQ=255)
                                                                                                                                          cactTCGTCGGCGACGAAGAGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCCAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATCCCCATGGGAAAGCTGGTCGTACGTGAc   >  2:415681/3‑150 (MQ=18)
                                                                                                                                           tctTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGAAGACCATCGCCAAGGGCCTGACCTCCGGTTACATCCCCATGGGCGGGGTGATCGTGCGTGACg  >  2:850226/1‑150 (MQ=255)
                                                                                                                                            ctTCGTCGGCGACGAAGAGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCg                                      <  2:1133283/113‑1 (MQ=255)
                                                                                                                                            ctTCGTCGGCGACGAAGAGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCg                                      >  1:1133283/1‑113 (MQ=255)
                                                                                                                                                 tcgCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATCCCCATGGGCgggg                 <  2:513925/129‑1 (MQ=255)
                                                                                                                                                 tcgCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATCCCCATGGGCgggg                 >  1:513925/1‑129 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
GAGAAGAAGATCCTCGAAGTCGGCGTGGACAACGTTGCCGCCTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATCCCCATGGGCGGGGTGATCGTGCGTGACG  >  NC_002947/2486076‑2486364

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