Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 743,432 A→C 9.1% intergenic (+111/+163) 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,4320AC 9.1% 296.2 / 19.1 99intergenic (+111/+163)PP_0635/PP_0636group II intron‑encoding maturase/cold shock DNA‑binding domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (41/49);  new base C (7/2);  total (48/51)
Fisher's exact test for biased strand distribution p-value = 8.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.12e-01

Only 100 of 107 total aligned reads displayed.

TACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGT  >  NC_002947/743284‑743578
                                                                                                                                                    |                                                                                                                                                  
tACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGaa                                                                                                                                                    <  5:75195/149‑1 (MQ=255)
tACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGaa                                                                                                                                                    >  3:188166/1‑149 (MQ=255)
tACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGaa                                                                                                                                                    <  7:213451/149‑1 (MQ=255)
     gATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAAc                                                                                                                                               <  3:312905/149‑1 (MQ=255)
      aTACGAATAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCACGCTATGTGAAGCAAcc                                                                                                                                              >  2:319631/1‑149 (MQ=255)
        aCGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCtt                                                                                                                                            >  5:134762/1‑149 (MQ=255)
         cGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCTAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCtt                                                                                                                                            <  5:133671/148‑1 (MQ=255)
          gATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCg                                                                                                                                          >  5:158991/1‑149 (MQ=255)
          gATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCg                                                                                                                                          <  6:134762/149‑1 (MQ=255)
          gATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTCCCATGTGTGCCCGCTATGTGAAGCAACCTTCg                                                                                                                                          >  6:195925/1‑149 (MQ=255)
           aTTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCgg                                                                                                                                         <  3:175481/149‑1 (MQ=255)
            ttAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGCGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGt                                                                                                                                        >  6:81077/1‑149 (MQ=255)
            ttAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGCGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGt                                                                                                                                        >  6:75794/1‑149 (MQ=255)
              aaCCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGtt                                                                                                                                       <  7:48816/148‑1 (MQ=255)
                 cGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGTCCGCTATGTAACGCAACCGAAGGTTGctt                                                                                                                                   >  6:353136/1‑147 (MQ=17)
                 cGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGt                                                                                                                                   >  2:63252/1‑149 (MQ=255)
                  ggCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGtt                                                                                                                                  <  5:245526/149‑1 (MQ=255)
                  ggCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGtt                                                                                                                                  <  7:159388/149‑1 (MQ=255)
                        gCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATa                                                                                                                             <  2:47792/148‑1 (MQ=255)
                           gcgTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCg                                                                                                                          >  3:286228/1‑148 (MQ=255)
                                 ccTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGacaca                                                                                                                    >  2:158216/1‑148 (MQ=255)
                                 ccTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGacaca                                                                                                                    >  2:163624/1‑148 (MQ=255)
                                  cTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGacaca                                                                                                                    >  7:232371/1‑147 (MQ=255)
                                  cTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGAcacaca                                                                                                                  <  4:286228/149‑1 (MQ=255)
                                     aaCCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGTCCGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGcacacaa                                                                                                                 <  5:353136/147‑2 (MQ=18)
                                     aaCCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGcg                                                                                                               >  3:172283/1‑149 (MQ=255)
                                     aaCCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGcg                                                                                                               >  2:34591/1‑149 (MQ=255)
                                     aaCCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGcg                                                                                                               >  8:110503/1‑149 (MQ=255)
                                        cgccgTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGcgtcg                                                                                                            >  1:206529/1‑149 (MQ=255)
                                          ccgTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTg                                                                                                          <  7:150867/149‑1 (MQ=255)
                                            gtgtACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGAc                                                                                                        <  1:319631/149‑1 (MQ=255)
                                                 cGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGcc                                                                                                   <  1:63252/149‑1 (MQ=255)
                                                  ggAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCg                                                                                                  >  8:81594/1‑149 (MQ=255)
                                                      ccGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCgggg                                                                                               >  7:222408/1‑148 (MQ=255)
                                                      ccGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGa                                                                                              <  1:158216/149‑1 (MQ=255)
                                                      ccGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGa                                                                                              <  1:163624/149‑1 (MQ=255)
                                                       cGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCaa                                                                                                                                                >  7:114350/1‑98 (MQ=255)
                                                       cGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCaa                                                                                                                                                <  8:114350/98‑1 (MQ=255)
                                                         tACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAac                                                                                           <  1:34591/149‑1 (MQ=255)
                                                           cGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAacta                                                                                         <  4:188166/149‑1 (MQ=255)
                                                            gCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAactac                                                                                        >  7:78112/1‑149 (MQ=255)
                                                             cACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAactact                                                                                       <  6:212174/149‑1 (MQ=255)
                                                               cGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAAACTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTg                                                                                     >  6:449/1‑149 (MQ=255)
                                                                     tgtgAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTcgcg                                                                               <  2:219994/149‑1 (MQ=255)
                                                                     tgtgAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTcgcg                                                                               <  5:195925/149‑1 (MQ=255)
                                                                       cgagaGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCgcggc                                                                             <  5:81077/148‑1 (MQ=255)
                                                                         agagGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCgt                                                                           >  2:58585/1‑149 (MQ=255)
                                                                           agGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCgtg                                                                          <  6:121372/148‑1 (MQ=255)
                                                                            ggACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGtgtg                                                                        <  3:96563/149‑1 (MQ=255)
                                                                            ggACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGtgtg                                                                        <  6:158991/149‑1 (MQ=255)
                                                                               cggcggGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAg                                                                     <  7:2076/149‑1 (MQ=255)
                                                                                   gggTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCaga                                                                 <  8:232371/149‑1 (MQ=255)
                                                                                   gggTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCaga                                                                 <  8:78112/149‑1 (MQ=255)
                                                                                     gtgtGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCg                                                                                                                          <  5:321757/90‑1 (MQ=255)
                                                                                     gtgtGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCg                                                                                                                          >  6:321757/1‑90 (MQ=255)
                                                                                     gtgtGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAgaga                                                               <  4:193240/149‑1 (MQ=255)
                                                                                          aaTCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTg                                                            >  1:250135/1‑147 (MQ=255)
                                                                                           aTCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTa                                                                                                                                 <  8:340128/77‑1 (MQ=255)
                                                                                              cGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAggg                                                       >  5:347202/1‑148 (MQ=255)
                                                                                               gCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAggg                                                       >  6:225134/1‑147 (MQ=255)
                                                                                                   cctACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGt                                                 >  4:76425/1‑149 (MQ=255)
                                                                                                    ctACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCgtgt                                                                         <  5:336100/124‑1 (MQ=255)
                                                                                                    ctACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCgtgt                                                                         >  6:336100/1‑124 (MQ=255)
                                                                                                    ctACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGtt                                                <  2:205174/149‑1 (MQ=255)
                                                                                                    ctACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGtt                                                <  3:76425/149‑1 (MQ=255)
                                                                                                     tACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGtt                                                >  7:35537/1‑148 (MQ=255)
                                                                                                       cTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGa                                              <  2:250135/148‑1 (MQ=255)
                                                                                                        tCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGGTACATAGCGGACACACCGCGTCGTGACTTGCCGGGGAAACTACTGGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGa                                              >  5:198794/1‑147 (MQ=255)
                                                                                                            ttGCGCGAAATAAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCaga                                                                 <  6:30838/124‑1 (MQ=255)
                                                                                                            ttGCGCGAAATAAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCaga                                                                 >  5:30838/1‑124 (MQ=255)
                                                                                                              gcgcGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGc                                       <  7:81594/148‑1 (MQ=255)
                                                                                                               cgcgAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGGTACATAGCGGCCACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGGGGAGCAGAGAGGGTAGGGGAGCGTTGACAGTGGccg                                     >  5:43398/1‑147 (MQ=255)
                                                                                                                  gAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGttt                                  <  1:58585/149‑1 (MQ=255)
                                                                                                                    aaTCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCa                                <  7:110503/149‑1 (MQ=255)
                                                                                                                     aTCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCaa                               <  8:35537/149‑1 (MQ=255)
                                                                                                                      tCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTc                                 >  3:67871/1‑146 (MQ=255)
                                                                                                                      tCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTc                                 <  4:67871/146‑1 (MQ=255)
                                                                                                                         gACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAagcta                           <  4:172283/149‑1 (MQ=255)
                                                                                                                           acgacgctgtgtgTCCGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGCACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAagctagc                         >  1:128514/8‑149 (MQ=21)
                                                                                                                             tGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCtt                       >  3:152665/1‑149 (MQ=255)
                                                                                                                                   gtgtgCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGa                                                     >  8:321109/1‑113 (MQ=255)
                                                                                                                                   gtgtgCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGa                                                     <  7:321109/113‑1 (MQ=255)
                                                                                                                                      tgCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCAc              <  4:304955/149‑1 (MQ=255)
                                                                                                                                        tccGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGCACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGc            >  1:57230/2‑149 (MQ=37)
                                                                                                                                        tccGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGCACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGc            >  4:33177/2‑149 (MQ=37)
                                                                                                                                        tccGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGCACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGc            >  8:251111/2‑149 (MQ=37)
                                                                                                                                         ccGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCgt                                                                           <  4:181203/85‑1 (MQ=255)
                                                                                                                                         ccGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCgt                                                                           >  3:181203/1‑85 (MQ=255)
                                                                                                                                         ccGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAagcta                           <  8:151039/133‑1 (MQ=37)
                                                                                                                                         ccGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAagcta                           >  7:151039/1‑133 (MQ=37)
                                                                                                                                          cGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGa                                                                                              >  2:32770/1‑65 (MQ=255)
                                                                                                                                          cGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGGTGACAGTGGCGGTTGAACGCTAGCGTGGCGGTCACGCtc          >  6:251777/1‑148 (MQ=255)
                                                                                                                                          cGCTATGTGAAGCAACCTTCGGTTGCGTTACATACCGGACACACAGCGTCGTGACTTGCCGGGGa                                                                                              <  1:32770/65‑1 (MQ=255)
                                                                                                                                          cGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGCACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCt           >  6:289071/1‑148 (MQ=37)
                                                                                                                                               tgtgAAGCAACCTTCGGTTTCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGc       <  5:449/147‑1 (MQ=255)
                                                                                                                                                   aaGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGc                    <  5:307310/130‑1 (MQ=255)
                                                                                                                                                   aaGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGc                    >  6:307310/1‑130 (MQ=255)
                                                                                                                                                   aaGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGt  >  1:160962/1‑148 (MQ=255)
                                                                                                                                                    aGCAACCTTCGGTTGCGTTACATAGCGGACTGACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGc            <  6:320859/137‑1 (MQ=255)
                                                                                                                                                    aGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGt  >  3:318571/1‑147 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGT  >  NC_002947/743284‑743578

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