Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 743,460 A→G 8.5% intergenic (+139/+135) 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,4600AG 8.5% 298.7 / 12.0 106intergenic (+139/+135)PP_0635/PP_0636group II intron‑encoding maturase/cold shock DNA‑binding domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (48/48);  new base G (5/4);  total (54/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

CCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGCTAGCGACG  >  NC_002947/743317‑743603
                                                                                                                                               |                                                                                                                                                  
ccTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGacaca                                                                                                                                                >  5:229054/1‑148 (MQ=255)
 cTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGacaca                                                                                                                                                >  7:237742/1‑147 (MQ=255)
      ccgccgTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGcgtc                                                                                                                                         <  3:125779/149‑1 (MQ=255)
       cgccgTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGTCCGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGcacacaaggcag                                                                                                                                        >  6:29686/1‑143 (MQ=18)
       cgccgTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGcgtcg                                                                                                                                        >  7:57603/1‑149 (MQ=255)
             gtACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGAGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACtt                                                                                                                                  >  4:371822/1‑149 (MQ=255)
              cacaccaCCGTGCGTACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTg                                                                                                                                 >  7:312584/7‑149 (MQ=255)
                cGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGTCCGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGcacacaaggcag                                                                                                                                        <  5:61397/140‑7 (MQ=21)
                cGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGTCCGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGcacacaaggcag                                                                                                                                        >  6:61397/1‑134 (MQ=21)
                  gAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCgg                                                                                                                             <  5:400395/149‑1 (MQ=255)
                    aCCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCgggg                                                                                                                           <  5:152612/149‑1 (MQ=255)
                     ccGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGa                                                                                                                          <  6:46500/149‑1 (MQ=255)
                      cGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGaa                                                                                                                         <  6:17958/149‑1 (MQ=255)
                       gTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGaaa                                                                                                                        <  8:57603/149‑1 (MQ=255)
                        tACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAac                                                                                                                       >  6:6140/1‑149 (MQ=255)
                        tACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAac                                                                                                                       >  3:395072/1‑149 (MQ=255)
                           gCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAactac                                                                                                                    >  7:165548/1‑149 (MQ=255)
                           gCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAactac                                                                                                                    >  2:133656/1‑149 (MQ=255)
                            cACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAactact                                                                                                                   <  5:301726/149‑1 (MQ=255)
                             aCGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACtt                                                                                                                  >  5:322323/1‑149 (MQ=255)
                              cGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTg                                                                                                                 <  7:95224/149‑1 (MQ=255)
                              cGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGAGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTg                                                                                                                 <  3:371822/149‑1 (MQ=255)
                               ggtggtGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCgg                                                                                                                             <  6:196686/136‑1 (MQ=255)
                               ggtggtGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCgg                                                                                                                             >  5:196686/1‑136 (MQ=255)
                                gtggtgTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGt                                                                                                               <  8:312584/149‑1 (MQ=255)
                                    tgtgAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTcgcg                                                                                                           <  5:278548/149‑1 (MQ=255)
                                      tgAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCgcggc                                                                                                         <  2:292769/149‑1 (MQ=255)
                                        agagGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCgt                                                                                                       >  2:370768/1‑149 (MQ=255)
                                        agagGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCgt                                                                                                       >  3:404395/1‑149 (MQ=255)
                                                 cggGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCag                                                                                              <  4:395072/149‑1 (MQ=255)
                                                  gggTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCaga                                                                                             <  5:6140/149‑1 (MQ=255)
                                                    gtgtGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCGga‑‑‑ga                                                                                        >  4:324504/1‑149 (MQ=255)
                                                    gtgtGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAga‑‑‑ga                                                                                        >  2:393496/1‑149 (MQ=255)
                                                           tCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAgg                                                                                 >  6:174109/1‑149 (MQ=255)
                                                             cGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAggg                                                                                <  5:307758/148‑1 (MQ=255)
                                                             cGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGNAggg                                                                                <  5:125462/148‑1 (MQ=255)
                                                              gCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCGGA‑‑‑GACTGTAGGGGa                                                                              <  3:324504/149‑1 (MQ=255)
                                                              gCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGAAACACAGCGTCGTGACTTGACGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAggg                                                                                >  7:77282/1‑147 (MQ=255)
                                                                ctcctACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACtt                                                                                                                                  >  6:171406/1‑98 (MQ=255)
                                                                ctcctACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACtt                                                                                                                                  <  5:171406/98‑1 (MQ=255)
                                                                 tcctACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGa                                                                              <  6:162592/146‑1 (MQ=255)
                                                                 tcctACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGa                                                                              >  5:162592/1‑146 (MQ=255)
                                                                      cTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGa                                                                       <  8:237742/148‑1 (MQ=255)
                                                                       tCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACtt                                                                                                                                  <  5:136822/91‑1 (MQ=255)
                                                                       tCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACtt                                                                                                                                  >  6:136822/1‑91 (MQ=255)
                                                                       tCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGa                                                                       >  8:299057/1‑147 (MQ=255)
                                                                         gATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGtgtg                                                                                                    <  7:305528/119‑1 (MQ=255)
                                                                         gATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGtgtg                                                                                                    >  8:305528/1‑119 (MQ=255)
                                                                         gATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAga‑‑‑ga                                                                                        <  6:361550/128‑1 (MQ=255)
                                                                         gATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGt                                                                   >  7:307838/1‑149 (MQ=255)
                                                                             gcgcGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCg                                                               >  4:236359/1‑149 (MQ=255)
                                                                              cgcgAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACTCAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCg                                                               <  8:225581/148‑1 (MQ=255)
                                                                              cgcgAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCgg                                                              >  5:383200/1‑149 (MQ=255)
                                                                              cgcgAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCg                                                               <  1:393496/148‑1 (MQ=255)
                                                                               gcgAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGCCACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGCCAGTGGCGGt                                                             >  6:307925/1‑149 (MQ=255)
                                                                               gcgAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGt                                                             >  7:254960/1‑149 (MQ=255)
                                                                               gcgAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGt                                                             >  5:297352/1‑149 (MQ=255)
                                                                                      cAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCCGGGAAACTACTTGGTCGCGGCGTGTGGAGCaga                                                                                             >  8:51145/1‑113 (MQ=255)
                                                                                        gACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAagcta                                                    <  4:404395/149‑1 (MQ=255)
                                                                                        gACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAagcta                                                    <  6:297352/149‑1 (MQ=255)
                                                                                              ccTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTgg                                              <  3:236359/149‑1 (MQ=255)
                                                                                                  gtgtgCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCg                                            <  5:174109/147‑1 (MQ=255)
                                                                                                     tgCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCAc                                       <  6:74368/149‑1 (MQ=255)
                                                                                                       tccGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGCACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCaa                                                        <  1:269523/129‑1 (MQ=25)
                                                                                                       tccGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGCACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCaa                                                        >  2:374940/2‑130 (MQ=25)
                                                                                                       tccGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGCACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCaa                                                        <  1:374940/129‑1 (MQ=25)
                                                                                                       tccGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGCACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCaa                                                        >  2:269523/2‑130 (MQ=25)
                                                                                                       tccGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGCACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCa                                        <  4:44483/145‑1 (MQ=37)
                                                                                                       tccGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGCACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCa                                        >  3:44483/2‑146 (MQ=37)
                                                                                                       cccGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGc                                     <  6:383200/149‑1 (MQ=255)
                                                                                                          gCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAgg                                                                                 >  4:238664/1‑102 (MQ=255)
                                                                                                          gCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAgg                                                                                 <  3:238664/102‑1 (MQ=255)
                                                                                                              tgtgAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGcc                               <  7:193945/148‑1 (MQ=255)
                                                                                                                tgAAGCAACCTTCGGTTGCGTTACATAGCGGACACAAAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAg                            >  1:67957/1‑149 (MQ=255)
                                                                                                                 gAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGt                           >  4:89983/1‑149 (MQ=255)
                                                                                                                  aaGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCgg                                                              >  2:144979/1‑113 (MQ=255)
                                                                                                                  aaGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCgg                                                              <  1:144979/113‑1 (MQ=255)
                                                                                                                  aaGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGa                             <  5:147149/146‑1 (MQ=255)
                                                                                                                  aaGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGt                           <  1:133656/148‑1 (MQ=255)
                                                                                                                  aaGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGt                           >  4:264126/1‑148 (MQ=255)
                                                                                                                    gCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAg                        <  7:299057/149‑1 (MQ=255)
                                                                                                                    gCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAg                        <  6:93909/149‑1 (MQ=255)
                                                                                                                      aaCCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGt                                                                                    >  6:108933/1‑87 (MQ=255)
                                                                                                                      aaCCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAgg                       >  2:169772/1‑148 (MQ=255)
                                                                                                                      aaCCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGc                      <  3:89983/149‑1 (MQ=255)
                                                                                                                             ggTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGa                <  8:307838/148‑1 (MQ=255)
                                                                                                                             ggTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGAt               >  6:5662/1‑149 (MQ=255)
                                                                                                                              gTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCTGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATc              <  6:229054/149‑1 (MQ=255)
                                                                                                                              gTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATc              <  3:264126/149‑1 (MQ=255)
                                                                                                                                   gTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGAGGAGTGTAGGGGAGCGTTGAAATTGGCGGTTTCAAGCTAGGTTGGCGGTCCCGCTAAGCCGAGTGAGGCATTCGATACt            >  1:127135/1‑149 (MQ=255)
                                                                                                                                   gTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGc          >  1:93321/1‑148 (MQ=255)
                                                                                                                                   gTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGc          >  1:15652/1‑148 (MQ=255)
                                                                                                                                    ttACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGc          >  6:307115/1‑147 (MQ=255)
                                                                                                                                    ttACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGAGTTGACAGTGGGGTTTGAAGGTTGCGTGGCGGTCACGCTTAGCCGAGTGACGCCTTCGAAGCTGCTac       >  8:219888/1‑148 (MQ=255)
                                                                                                                                      aCATAGCGCACACACAGCCTCGTGACTTGCCGGGGAAACTACTTTGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTACGGGACCGTTGACAGTCGCGGTTTCAATCTAGCATGCCGGTCACGCAAAGCCGACTCAGCCATACGATCCTCCTAGc      >  4:204906/1‑149 (MQ=255)
                                                                                                                                          aGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGa                             >  7:71426/1‑122 (MQ=255)
                                                                                                                                          aGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGa                             <  8:71426/122‑1 (MQ=255)
                                                                                                                                           gCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGCTAGCGACg  >  1:351574/1‑148 (MQ=255)
                                                                                                                                               acacacAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGc          >  2:320541/1‑136 (MQ=255)
                                                                                                                                               acacacAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGc          <  1:320541/136‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
CCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGA‑‑‑GACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGGTCACGCTAAGCCGAGTGAGGCATTCGATCCTGCTAGCGACG  >  NC_002947/743317‑743603

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: