Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,349,034 T→G 5.6% intergenic (‑399/+75) PP_3820 ← / ← galU group II intron‑encoding maturase/UTP‑glucose‑1‑phosphate uridylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,349,0340TG 5.6% 255.1 / 2.2 81intergenic (‑399/+75)PP_3820/galUgroup II intron‑encoding maturase/UTP‑glucose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (40/36);  new base G (2/3);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGATGTAGCCTTC  >  NC_002947/4348890‑4349174
                                                                                                                                                |                                                                                                                                               
gCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTAcg                                                                                                                                             >  8:224748/1‑149 (MQ=32)
      aTTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAc                                                                                                                                       <  2:110800/149‑1 (MQ=32)
      aTTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCGTAAGCGCAAc                                                                                                                                       <  7:185486/149‑1 (MQ=17)
      aTTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGCTTGCGCTTACGCGCAAc                                                                                                                                       >  4:207693/1‑149 (MQ=17)
         ttGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAc‑‑‑gc                                                                                                                                  <  8:122782/148‑1 (MQ=32)
             cTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACAACCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCAATTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCAt                                                                                                                             >  4:154827/1‑149 (MQ=17)
               cGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGc                                                                                                                           >  1:192736/1‑149 (MQ=32)
                gggCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGCGTTGCGCGTAAGCGCAAC‑‑‑ccgcttgg                                                                                                                            <  6:66109/147‑5 (MQ=11)
                       cttcAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGca                                                                                                                   >  5:320201/1‑149 (MQ=32)
                        ttcAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGca                                                                                                                   >  6:309933/1‑148 (MQ=32)
                          cAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGcaca                                                                                                                 >  6:359870/1‑148 (MQ=32)
                          cAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGcaca                                                                                                                 <  4:218456/148‑1 (MQ=32)
                           aGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCAGGCATGgcgcgttg                                                                                                                <  5:3201/148‑5 (MQ=12)
                            gATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGcacaca                                                                                                               >  7:108867/1‑148 (MQ=32)
                             aTTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACaa                                                                                                              >  3:352630/1‑148 (MQ=32)
                               tCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACaaaa                                                                                                            >  5:14583/1‑148 (MQ=32)
                                         ccGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCAcccc                                                                                                  >  3:346545/1‑148 (MQ=35)
                                           gcgACAACCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCAATTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGaa                                                                                               <  3:154827/149‑1 (MQ=255)
                                               cacCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGAt                                                                                           >  7:345805/1‑149 (MQ=255)
                                                acCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATg                                                                                          <  3:181317/149‑1 (MQ=255)
                                                 cccTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATgg                                                                                         <  8:29747/149‑1 (MQ=255)
                                                  ccTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGc                                                                                        >  1:250227/1‑149 (MQ=255)
                                                  ccTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGc                                                                                        >  4:21666/1‑149 (MQ=255)
                                                   cTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCt                                                                                       <  4:147808/149‑1 (MQ=255)
                                                   cTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCt                                                                                       <  7:224748/149‑1 (MQ=255)
                                                      acagccGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTg                                                                                    <  2:326151/145‑1 (MQ=255)
                                                       ccTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGa                                                                                   <  8:132153/149‑1 (MQ=255)
                                                        cTCCGGCTAGCACTTACCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAg                                                                                  <  7:110984/149‑1 (MQ=255)
                                                         tCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCCCTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAgg                                                                                 >  3:58680/1‑149 (MQ=255)
                                                         tCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAgg                                                                                 >  6:17078/1‑149 (MQ=255)
                                                           cGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGc                                                                                <  4:346545/148‑1 (MQ=255)
                                                               tAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTCGGGGGGGGGGCAACGGAGGGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGc                                                                                >  1:119937/1‑147 (MQ=17)
                                                                 gCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAgg                                                                         >  4:341479/1‑149 (MQ=255)
                                                                   aCTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGc                                                                       >  5:195190/1‑149 (MQ=255)
                                                                        cccTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTAt                                                                   <  8:345805/148‑1 (MQ=255)
                                                                         ccTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCa                                                                 <  6:320201/149‑1 (MQ=255)
                                                                           ttGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAg                                                                >  5:380710/1‑148 (MQ=255)
                                                                                  gtgtgtAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCtt                                                        <  5:359870/149‑1 (MQ=255)
                                                                                  gtgtgtAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCtt                                                        >  2:193055/1‑149 (MQ=255)
                                                                                     tgtAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGc                                                      <  3:341479/148‑1 (MQ=255)
                                                                                       tAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGc                                                      <  1:204037/146‑1 (MQ=255)
                                                                                       tAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGc                                                      >  2:204037/1‑146 (MQ=255)
                                                                                       tAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCa                                                    <  7:317939/148‑1 (MQ=255)
                                                                                         gagGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTc                                                 >  6:57768/1‑149 (MQ=255)
                                                                                         gagGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAACCCTTGCCAGTc                                                 <  8:108867/149‑1 (MQ=255)
                                                                                           ggACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCtt                                               <  5:17078/149‑1 (MQ=255)
                                                                                           ggACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCtt                                               >  3:15083/1‑149 (MQ=255)
                                                                                               tttCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTc                                                 >  3:22482/1‑143 (MQ=255)
                                                                                               tttCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTc                                                 <  4:22482/143‑1 (MQ=255)
                                                                                               tttCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTgtag                                           >  5:233033/1‑149 (MQ=255)
                                                                                               tttCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGAAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTgtag                                           >  1:102090/1‑149 (MQ=255)
                                                                                                 tCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTgtag                                           <  6:380710/147‑1 (MQ=255)
                                                                                                 tCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGAAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTgtag                                           <  5:309933/147‑1 (MQ=255)
                                                                                                   aCCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGtagtag                                        <  5:57768/148‑1 (MQ=255)
                                                                                                    ccTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTg                                              <  6:251410/141‑1 (MQ=255)
                                                                                                    ccTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTg                                              >  5:251410/1‑141 (MQ=255)
                                                                                                          aGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGaa                                <  2:192736/149‑1 (MQ=255)
                                                                                                            tCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAg                               >  8:349782/1‑148 (MQ=255)
                                                                                                            tCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAg                               >  4:328177/1‑148 (MQ=255)
                                                                                                                ccACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTc                                                 <  4:329290/126‑1 (MQ=255)
                                                                                                                ccACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTc                                                 >  3:329290/1‑126 (MQ=255)
                                                                                                                  aCTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCTGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGt                        >  1:72449/1‑149 (MQ=255)
                                                                                                                  aCTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGt                        >  2:29268/1‑149 (MQ=255)
                                                                                                                      ggCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTg                      <  2:119937/147‑1 (MQ=255)
                                                                                                                       gCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTc                   <  6:195190/149‑1 (MQ=255)
                                                                                                                         caccacAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTcgc                 >  8:117658/1‑149 (MQ=255)
                                                                                                                            cacaACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGaa                          >  1:56381/1‑137 (MQ=255)
                                                                                                                            cacaACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTc              >  6:124057/1‑149 (MQ=255)
                                                                                                                               aaCCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGAt           >  7:141780/1‑149 (MQ=255)
                                                                                                                               aaCCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGAt           >  1:244160/1‑149 (MQ=255)
                                                                                                                                 ccAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTAt                                                                   >  4:269223/1‑91 (MQ=255)
                                                                                                                                 ccAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTAt                                                                   <  3:269223/91‑1 (MQ=255)
                                                                                                                                 ccAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGATGt         <  4:15083/149‑1 (MQ=255)
                                                                                                                                 ccAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGATGt         <  7:349782/149‑1 (MQ=255)
                                                                                                                                  cAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGATGTa        <  7:117658/149‑1 (MQ=255)
                                                                                                                                  atggcgcGTTGCGCGTAAGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGaa                                <  2:98235/122‑1 (MQ=255)
                                                                                                                                  atggcgcGTTGCGCGTAAGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGaa                                >  1:98235/4‑125 (MQ=255)
                                                                                                                                   aaGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGATTATCAGTAAGCCTTGCGAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGATGTAg       >  4:231651/1‑149 (MQ=255)
                                                                                                                                    aGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGATGTAGc      <  4:79065/149‑1 (MQ=255)
                                                                                                                                    aGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGATGTAGc      >  7:289765/1‑149 (MQ=255)
                                                                                                                                        ggTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGATGTAGCCTTc  <  8:141780/149‑1 (MQ=255)
                                                                                                                                                ttACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAg                         >  4:312212/1‑118 (MQ=255)
                                                                                                                                                ttACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAg                         <  3:312212/118‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GCCAGCATTTTGCCTCGGGCTTCCTTCAGATTCGCGGTCACCCGCGACACCCTTGCCTCCGGCTAGCACTTCCCCTTGCCGGGTGTGTAGAGGACTTTCACCTCCAAGTCACCCACTTGGCCACCACAACCAAGCGGGTTGCGCTTACGCGCAAC‑‑‑GCGCCATGCCTGGCGCACACAAAAAAGCCACCCCGAAGGATGGCTTTGAGGCCAGCAGGGCTTATCAGTAAGCCTTGCCAGTCTTGTAGTAGTTCTCGAAGCAGAAGTTGGTCGCGTCGATGTAGCCTTC  >  NC_002947/4348890‑4349174

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: