Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 743,433 G→A 9.2% intergenic (+109/+165) PP_0635 → / ← PP_0636 group II intron‑encoding maturase/cold shock DNA‑binding domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947743,4330GA 9.2% 447.7 / 25.7 153intergenic (+109/+165)PP_0635/PP_0636group II intron‑encoding maturase/cold shock DNA‑binding domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (62/76);  new base A (6/8);  total (69/84)
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.99e-01

Only 100 of 159 total aligned reads displayed.

CGATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTG  >  NC_002947/743284‑743582
                                                                                                                                                     |                                                                                                                                                     
cGATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTAtgtg                                                                                                                                                       <  2:5222871/150‑1 (MQ=255)
 gATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGa                                                                                                                                                      >  1:205531/1‑150 (MQ=255)
  aTACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGaa                                                                                                                                                     <  2:3806683/150‑1 (MQ=255)
  aTACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGaa                                                                                                                                                     <  2:447810/150‑1 (MQ=255)
        gATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAAcc                                                                                                                                               >  1:898084/1‑150 (MQ=255)
        gATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAAcc                                                                                                                                               >  2:639013/1‑150 (MQ=255)
        gATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAAcc                                                                                                                                               >  2:5393550/1‑150 (MQ=255)
            cGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCg                                                                                                                                           <  1:5441436/150‑1 (MQ=255)
            cGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCg                                                                                                                                           <  1:5312852/150‑1 (MQ=255)
            cGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCg                                                                                                                                           <  1:4693876/150‑1 (MQ=255)
             gATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTCCTCGATTGCGCGAAATCAGACCTGCCGTGTGTGCCCGCTATGTGAAGCAACCTTCgg                                                                                                                                          >  2:3329837/1‑150 (MQ=255)
             gATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCgg                                                                                                                                          >  1:4461648/1‑150 (MQ=255)
             gATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCgg                                                                                                                                          <  1:3227944/150‑1 (MQ=255)
               ttAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGtt                                                                                                                                        <  1:2541041/150‑1 (MQ=255)
                  aCCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGaa                                                                                                                                                     >  1:4764499/1‑134 (MQ=255)
                  aCCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGaa                                                                                                                                                     <  2:4764499/134‑1 (MQ=255)
                   ccGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGt                                                                                                                                    <  1:444399/150‑1 (MQ=255)
                    cGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGtt                                                                                                                                   <  2:810437/150‑1 (MQ=255)
                    cGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGtt                                                                                                                                   <  2:1684985/150‑1 (MQ=255)
                    cGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTAACGCAACCGAAGGTTGcttc                                                                                                                                   <  1:3739329/150‑4 (MQ=18)
                     ggCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTTCTCGATTGCGGGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGttt                                                                                                                                  <  2:3002282/150‑2 (MQ=255)
                     ggCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTa                                                                                                                                  >  2:2106701/1‑150 (MQ=255)
                     ggCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGGGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCGTCGGTTGCGTTa                                                                                                                                  >  2:2588215/1‑150 (MQ=255)
                           gCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGc                                                                                                                            <  2:2373974/150‑1 (MQ=255)
                           gCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACAGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGTCCGCTATGTAACGCAACCGAAGGTTGCTTCACATAGc                                                                                                                            <  2:4703926/150‑1 (MQ=18)
                              gcgTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGa                                                                                                                         >  1:3054448/1‑150 (MQ=255)
                              gcgTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGa                                                                                                                         <  1:3580194/150‑1 (MQ=255)
                                    ccTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGc                                                                                                                            <  1:1306064/141‑1 (MQ=255)
                                     cTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAg                                                                                                                  >  2:2186353/1‑150 (MQ=255)
                                      tGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGc                                                                                                                 >  1:4034806/1‑150 (MQ=255)
                                       gtaCCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCGTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACTTAGCGGACACACAGcg                                                                                                                >  2:3779450/3‑150 (MQ=255)
                                       gAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATCGCGGACACACAGcg                                                                                                                >  2:594277/1‑150 (MQ=255)
                                           cgccgTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACGGcgtcgt                                                                                                            >  2:2705610/1‑150 (MQ=255)
                                                     ngAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCgg                                                                                                  >  1:1595017/2‑150 (MQ=255)
                                                     ggAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCgg                                                                                                  >  1:4850083/1‑150 (MQ=255)
                                                     ggAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCgg                                                                                                  >  2:3928638/1‑150 (MQ=255)
                                                     ggAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCgg                                                                                                  >  2:2920318/1‑150 (MQ=255)
                                                     ggAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCgg                                                                                                  >  1:4164379/1‑150 (MQ=255)
                                                     ggAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCgg                                                                                                  >  1:92465/1‑150 (MQ=255)
                                                       aaCCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCgggg                                                                                                <  2:3668877/150‑1 (MQ=255)
                                                           gTACGCACGGTGGTGTGAGAGGTCGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAac                                                                                            <  2:365261/150‑1 (MQ=255)
                                                           gTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAac                                                                                            <  1:1009128/150‑1 (MQ=255)
                                                                 aCGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGACACTACTTg                                                                                      >  2:5450226/1‑150 (MQ=255)
                                                                 aCGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTg                                                                                      <  1:4342955/150‑1 (MQ=255)
                                                                     tggtgTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTcg                                                                                  <  1:276727/150‑1 (MQ=255)
                                                                      ggtgTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTcgc                                                                                 >  2:4179283/1‑150 (MQ=255)
                                                                          tgAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACa                                                                                                                                <  1:4591092/99‑1 (MQ=255)
                                                                          tgAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACa                                                                                                                                >  2:4591092/1‑99 (MQ=255)
                                                                          tgAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCgcggcg                                                                             >  2:4466200/1‑150 (MQ=255)
                                                                             gagGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCgtgt                                                                          <  1:3732187/150‑1 (MQ=255)
                                                                              agGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGtgtg                                                                         <  2:898084/150‑1 (MQ=255)
                                                                              agGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACGTGGTCGCGGCGtgtg                                                                         <  2:4850083/150‑1 (MQ=255)
                                                                               ggACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTgg                                                                        >  2:228567/1‑150 (MQ=255)
                                                                                 aCGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAg                                                                      <  2:3201429/150‑1 (MQ=255)
                                                                                  cggcggGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGc                                                                     <  1:4309824/150‑1 (MQ=255)
                                                                                  cggcggGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGc                                                                     <  1:5185946/150‑1 (MQ=255)
                                                                                   ggcggGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCa                                                                    <  2:1030620/150‑1 (MQ=255)
                                                                                     cggGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCaga                                                                  <  2:1612971/150‑1 (MQ=255)
                                                                                     cggGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCaga                                                                  <  2:2711004/150‑1 (MQ=255)
                                                                                            gAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACGACTTGGTCGCGGCGTGTGGAGCAGAGACTGTa                                                           >  2:392356/1‑150 (MQ=255)
                                                                                             aaTCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAg                                                          <  2:5066175/150‑1 (MQ=255)
                                                                                               tCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAggg                                                        <  1:846354/150‑1 (MQ=255)
                                                                                               tCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAggg                                                        <  1:3541114/150‑1 (MQ=255)
                                                                                               tCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCATTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAggg                                                        <  2:2500576/150‑1 (MQ=255)
                                                                                                      cctACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGtt                                                 <  1:639013/150‑1 (MQ=255)
                                                                                                        tACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTTGAGCAGAGACTGTAGGGGAGCGTTGa                                               >  2:925400/1‑150 (MQ=255)
                                                                                                        tACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGa                                               >  2:1027739/1‑150 (MQ=255)
                                                                                                        tACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGa                                               >  1:3721740/1‑150 (MQ=255)
                                                                                                         aCTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGAc                                              <  2:5354554/150‑1 (MQ=255)
                                                                                                          cTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACa                                             >  2:5139761/1‑150 (MQ=255)
                                                                                                          cTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACa                                             >  2:1249018/1‑150 (MQ=255)
                                                                                                          cTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACa                                             >  2:1248900/1‑150 (MQ=255)
                                                                                                                 gcgcGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGctg                                      <  2:1186511/150‑3 (MQ=255)
                                                                                                                  cgcgAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGt                                     >  2:1406278/1‑150 (MQ=255)
                                                                                                                  cgcgAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGt                                     <  2:2041064/150‑1 (MQ=255)
                                                                                                                     gAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAATCAACCTTAGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCAGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACGGTAGGGGAGCGTTGACAGTGGCGGTTTc                                  >  2:1842136/1‑150 (MQ=255)
                                                                                                                       aaTCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCaa                                <  2:4332220/150‑1 (MQ=255)
                                                                                                                       aaTCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCaa                                <  2:5302510/150‑1 (MQ=255)
                                                                                                                         tCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAagc                              <  2:245818/150‑1 (MQ=255)
                                                                                                                           aGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTctttggc                            <  2:3194677/150‑7 (MQ=255)
                                                                                                                           aGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTctttggc                            <  2:467605/150‑7 (MQ=255)
                                                                                                                                  ccTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGc                     <  1:345642/150‑1 (MQ=255)
                                                                                                                                   cTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCg                    <  1:844971/150‑1 (MQ=255)
                                                                                                                                   cTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCg                    <  1:3589583/150‑1 (MQ=255)
                                                                                                                                     tgtgtgCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGt                  <  1:1510950/150‑1 (MQ=255)
                                                                                                                                         tgtccGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGCACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACg              >  1:3451724/4‑150 (MQ=37)
                                                                                                                                         tgtccGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGCACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACg              >  2:5169287/4‑150 (MQ=37)
                                                                                                                                         tgCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACg              >  1:5314100/1‑150 (MQ=255)
                                                                                                                                         tgCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACg              <  2:1198098/150‑1 (MQ=255)
                                                                                                                                          gCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACCGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGc             >  2:3870730/1‑150 (MQ=255)
                                                                                                                                          gCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGc                     >  1:1227494/1‑142 (MQ=255)
                                                                                                                                           cccGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCt            >  1:2370851/1‑150 (MQ=255)
                                                                                                                                            ccGCTATGTAACGCAACCGAGGGTTGCTTCACATAGCGGGCACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACt                                                              <  1:5369948/99‑1 (MQ=21)
                                                                                                                                            ccGCTATGTAACGCAACCGAGGGTTGCTTCACATAGCGGGCACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACt                                                              >  2:5369948/1‑99 (MQ=21)
                                                                                                                                             cGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGCACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCaa                                <  2:1764090/128‑1 (MQ=37)
                                                                                                                                             cGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGCACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCaa                                >  1:1764090/1‑128 (MQ=37)
                                                                                                                                              gCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACtt                                                                                       <  1:1009743/72‑1 (MQ=14)
                                                                                                                                              gCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACtt                                                                                       >  2:1009743/1‑72 (MQ=14)
                                                                                                                                                    tgAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGt   <  1:2106701/150‑1 (MQ=255)
                                                                                                                                                     gAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTg  <  1:4497056/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTG  >  NC_002947/743284‑743582

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