Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,486,224 C→G 6.4% 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.4% 234.1 / 4.7 78A258G (GCC→GGC) spuC‑Ipolyamine:pyruvate transaminase
Reads supporting (aligned to +/- strand):  ref base C (37/35);  new base G (3/2);  total (40/38)
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.37e-01

GAAGAAGATCCTCGAAGTCGGCGTGGACAACGTTGCCGCCTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATCCCCATGGGCGGGGTGATCGTGCGTGACGAG  >  NC_002947/2486078‑2486366
                                                                                                                                                  |                                                                                                                                              
gaagaagaTCCTCGAAGTTAGCCTGGACAACGTTGCCGCCTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCcga                                                                                                                                             <  1:97553/150‑1 (MQ=255)
      gaTCCTCGAAGTCGGCGTGGACAACGTTGCCGCCTTCATCGCCGAGCCTATCCAGGGCGCTGGTGGCGTGATCATTCCGCCGAAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGGCGACGAAga                                                                                                                                       <  1:289023/150‑2 (MQ=255)
      gaTCCTCGAAGTCGGCGTGGACAACGTTGCCGCCTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGt                                                                                                                                       <  1:1491871/150‑1 (MQ=255)
       aTCCTCGAAGTCGGCGTGGACAACGTTGCCGCCTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTg                                                                                                                                      >  1:1074694/1‑150 (MQ=255)
           tCGAAGTCGGCGTGGACAACGTTGCCGCCTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCt                                                                                                                                  >  1:1622138/1‑150 (MQ=255)
             gAAGTCGGCGTGGACAACGTTGCCGCCTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATAATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATNTGc                                                                                                                                <  2:1193657/150‑1 (MQ=255)
                gTCGGCGTGGACAACGTTGCCGCTTTCATCGCCGAGCCGCTCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGGTCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCggg                                                                                                                             <  2:127307/150‑1 (MQ=255)
                 tCGGCGTGGACAACGTTGCCGCCTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCCTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAAGCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGt                                                                                                                            >  2:1000636/1‑150 (MQ=255)
                         gACAACGTTGCCGCCTTCATCGCTGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGc                                                                                                                                >  2:944493/1‑138 (MQ=255)
                         gACAACGTTGCCGCCTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGc                                                                                                                                <  1:944493/138‑1 (MQ=255)
                                    cgccTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTg                                                                                                         >  2:170100/1‑150 (MQ=255)
                                    cgccTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTg                                                                                                         >  1:32309/1‑150 (MQ=255)
                                        ttCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTc                                                                                                     >  1:1846593/1‑150 (MQ=255)
                                           aTCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTAGCGGTGAGTGGTTCGGt                                                                                                  >  1:1625725/1‑150 (MQ=255)
                                             cGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTAc                                                                                                <  2:1352187/150‑1 (MQ=255)
                                             cGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTAc                                                                                                <  1:1195935/150‑1 (MQ=255)
                                              gccgagccgaTCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTAcc                                                                                               >  1:6041/1‑150 (MQ=255)
                                                ctcggcgaTCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCg                                                                                                                         >  2:541138/6‑122 (MQ=255)
                                                ctcggcgaTCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCg                                                                                                                         <  1:541138/117‑1 (MQ=255)
                                                 gagccgaTCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGCGGTTCGGTACCGat                                                                                            >  2:639897/1‑150 (MQ=255)
                                                  agccgaTCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGatt                                                                                           <  2:113772/150‑1 (MQ=255)
                                                       aTCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGa                                                                                                            <  2:1016085/128‑1 (MQ=255)
                                                       aTCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGa                                                                                                            >  1:1016085/1‑128 (MQ=255)
                                                       aTCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTAc                                                                                      <  1:170100/150‑1 (MQ=255)
                                                       aTCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTAc                                                                                      <  2:1458098/150‑1 (MQ=255)
                                                         ccAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGa                                                                                    <  2:32309/150‑1 (MQ=255)
                                                         ccAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGa                                                                                    >  2:1167133/1‑150 (MQ=255)
                                                          cAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGAc                                                                                   <  2:1622138/150‑1 (MQ=255)
                                                           aGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGAcc                                                                                  >  2:1301347/1‑150 (MQ=255)
                                                           aGGGCGCCGGTGGCGTGAGCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGAcc                                                                                  <  2:1376347/150‑1 (MQ=255)
                                                               cgcCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCaa                                                                              <  2:488846/150‑1 (MQ=255)
                                                               cgcCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCaa                                                                              <  2:1074694/150‑1 (MQ=255)
                                                                  cGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCANGcc                                                                           <  2:1200862/150‑1 (MQ=255)
                                                                   ggTGGCGTGATCATTCCGCCACAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGcct                                                                          >  2:1824161/1‑150 (MQ=255)
                                                                    gTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGt                                                                                                                                       <  1:440177/88‑1 (MQ=255)
                                                                    gTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGt                                                                                                                                       >  2:440177/1‑88 (MQ=255)
                                                                     tGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGcctga                                                                        <  2:1113425/150‑1 (MQ=255)
                                                                         gTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGcctgacctg                                                                    >  1:46781/1‑150 (MQ=255)
                                                                            aTCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGa                                                                                    <  1:857441/131‑1 (MQ=255)
                                                                            aTCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGa                                                                                    >  2:857441/1‑131 (MQ=255)
                                                                             tCATTCCGCCACAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCtgatga                                                                >  1:442851/1‑150 (MQ=255)
                                                                                    gccgcAGAGCTATTGGCAGAAGAGCAAGGAAATCCTGGCCAGGTACGACAGCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGc                                                         <  2:305322/150‑1 (MQ=255)
                                                                                        cAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAg                                                     <  1:1167133/150‑1 (MQ=255)
                                                                                         aGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCTGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAgg                                                    >  2:1015444/1‑150 (MQ=255)
                                                                                             cTATTGGCCGAAGATCAACGAAATCCTGGTCAGGTACGACATCCCCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGCCCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCt                                                <  1:579962/150‑1 (MQ=255)
                                                                                                ttGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGAc                                             >  2:792008/1‑150 (MQ=255)
                                                                                                   gtcGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGaa                                                                                                                                         >  2:879342/3‑55 (MQ=255)
                                                                                                   gtcGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGaa                                                                                                                                         >  2:881275/3‑55 (MQ=255)
                                                                                                   gtcGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGaa                                                                                                                                         <  1:881275/53‑1 (MQ=255)
                                                                                                   gtcGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGaa                                                                                                                                         <  1:879342/53‑1 (MQ=255)
                                                                                                   gCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTAc                                                                                      <  2:57159/106‑1 (MQ=255)
                                                                                                   gCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTAc                                                                                      >  1:57159/1‑106 (MQ=255)
                                                                                                     cGAAGTCCAAGGAAATGCTGGGCAGGTACGGGAGCGTCTTCGAGGCCGAGGAAGGGATCTGGGGGTTTGGCCGTACGGGTGAGTGGTTCGGTAGCGGTGAGTAGGACCTCGAGCCTGCGCTGGTGACCGTCGCCAAGGGCCTGACCTCCg                                        <  2:565628/150‑1 (MQ=255)
                                                                                                      gAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCgg                                       <  1:792008/150‑1 (MQ=255)
                                                                                                       aaGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGt                                      <  1:1419628/150‑1 (MQ=255)
                                                                                                          aTCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTAc                                   >  2:866342/1‑150 (MQ=255)
                                                                                                             aaGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATc                                >  1:1428031/1‑150 (MQ=255)
                                                                                                              aGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATcc                               >  2:1174002/1‑150 (MQ=255)
                                                                                                               ggAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCaa                                                                              >  1:537692/1‑102 (MQ=255)
                                                                                                               ggAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCaa                                                                              <  2:537692/102‑1 (MQ=255)
                                                                                                               ggAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATccc                              >  2:1285506/1‑150 (MQ=255)
                                                                                                               ggAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATccc                              >  2:91385/1‑150 (MQ=255)
                                                                                                                gAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATcccc                             >  1:802148/1‑150 (MQ=255)
                                                                                                                gAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATcccc                             >  2:1913903/1‑150 (MQ=255)
                                                                                                                gAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATcccc                             <  2:172004/150‑1 (MQ=255)
                                                                                                                  aaTCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATCCCCAt                           <  2:6041/150‑1 (MQ=255)
                                                                                                                        ggCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCgg                                                                                                                              >  2:1108469/1‑45 (MQ=255)
                                                                                                                        ggCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCgg                                                                                                                              <  1:1108469/45‑1 (MQ=255)
                                                                                                                          ccAGGTACGACATCCTCTTCGTCGACGCCGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGCTAATTACGACCTCAAGCCTGAGTTCTCCCGTGTTTGCAGGTGCCTGACCTCCGGTTACATCCCCATGGGCgggg                   <  1:1773698/150‑1 (MQ=18)
                                                                                                                          ccAGGTACGACATCCTCTTCGGCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATCCCCATGGGCgggg                   >  2:1459903/1‑150 (MQ=255)
                                                                                                                                    agataactTCGTCGGCGACGAAGAGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATCCCCATGGGCggg                    >  2:1168359/7‑139 (MQ=255)
                                                                                                                                    agaccactTCGTCGGCGACGAAGAGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATCCCCATGGGCggg                    <  1:1168359/133‑1 (MQ=255)
                                                                                                                                        ctctTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATCCCCATGGGCGGGGTGATCGTGCGTGAc     >  2:1023675/1‑150 (MQ=255)
                                                                                                                                        ctctTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATCCCCATGGGCGGGGTGATCGTGCGTGAc     <  1:625623/150‑1 (MQ=255)
                                                                                                                                         tctTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATCCCCATGGGCGGGGTGATCGTGCGTGACg    >  1:884086/1‑150 (MQ=255)
                                                                                                                                          ctTCGTCGGCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATCCCCATGGGCGGGGTGATCGTGCGTGACGa   >  1:1409698/1‑150 (MQ=255)
                                                                                                                                          ctTCGTCGGCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATCCCCATGGGCGGGGTGATCGTGCGTGACGa   >  2:1824396/1‑150 (MQ=255)
                                                                                                                                           ttCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATCCCCATGGGCGGGGTGATCGTGCGGGACGAg  >  1:1717348/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
GAAGAAGATCCTCGAAGTCGGCGTGGACAACGTTGCCGCCTTCATCGCCGAGCCGATCCAGGGCGCCGGTGGCGTGATCATTCCGCCGCAGACCTATTGGCCGAAGATCAAGGAAATCCTGGCCAGGTACGACATCCTCTTCGTCGCCGACGAAGTGATCTGCGGGTTCGGCCGTACCGGTGAGTGGTTCGGTACCGATTATTACGACCTCAAGCCTGACCTGATGACCATCGCCAAGGGCCTGACCTCCGGTTACATCCCCATGGGCGGGGTGATCGTGCGTGACGAG  >  NC_002947/2486078‑2486366

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