Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,425,650 T→C 6.6% intergenic (+52/‑453) PP_1249 → / → PP_1250 DUF4223 domain‑containing protein/group II intron‑encoding maturase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,425,6500TC 6.6% 262.6 / 5.4 90intergenic (+52/‑453)PP_1249/PP_1250DUF4223 domain‑containing protein/group II intron‑encoding maturase
Reads supporting (aligned to +/- strand):  ref base T (42/40);  new base C (3/3);  total (47/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

TGCACCTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGG  >  NC_002947/1425506‑1425798
                                                                                                                                                |                                                                                                                                                     
tGCACCTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATg                                                                                                                                                   >  6:292912/1‑149 (MQ=255)
  cACCTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATg                                                                                                                                                   <  6:336294/147‑1 (MQ=255)
     cTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACAt                                                                                                                                              >  5:141049/1‑149 (MQ=255)
     cTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACAt                                                                                                                                              >  8:57483/1‑149 (MQ=255)
     cTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACAt                                                                                                                                              >  8:41636/1‑149 (MQ=255)
      tACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATg                                                                                                                                             >  3:14395/1‑149 (MQ=255)
        cGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGc                                                                                                                                           <  2:28535/149‑1 (MQ=255)
          aCTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCg                                                                                                                                         <  2:192147/149‑1 (MQ=255)
          aCTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGCCATGTGCATGACCg                                                                                                                                         <  4:229529/149‑1 (MQ=255)
               cTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGt                                                                                                                                    >  4:136083/1‑149 (MQ=255)
               cTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGt                                                                                                                                    >  4:130511/1‑149 (MQ=255)
               cTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGGCATGCACATGGCCGGGGGt                                                                                                                                    >  4:51715/1‑149 (MQ=255)
                    cgtcgcTTCGGTTTCTATCTCGAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCTGCTGGTGTAATGTTGGCGGTAGCCTGATCGCGACGCTTGTGTCGTGTTCCAAGGCCGCCTGTCATGCTGATGGGCGGGGGttttt                                                                                                                                <  4:160789/143‑1 (MQ=255)
                    ccACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGttttt                                                                                                                                <  7:198066/148‑1 (MQ=255)
                    ccACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGttttt                                                                                                                                <  4:8152/148‑1 (MQ=255)
                        ccTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTtgtgtg                                                                                                                           >  5:82561/1‑149 (MQ=255)
                          ttCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTgcg                                                                                                                         >  6:47353/1‑149 (MQ=255)
                             ggTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGcc                                                                                                                       <  7:312783/148‑1 (MQ=255)
                                ttCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAgg                                                                                                                    >  4:154505/1‑148 (MQ=255)
                                    aTCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATgg                                                                                                               >  4:48845/1‑149 (MQ=255)
                                     tctcCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGc                                                                                                                                           >  7:192016/1‑120 (MQ=255)
                                     tctcCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGc                                                                                                                                           <  8:192016/120‑1 (MQ=255)
                                      ctcCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGgcg                                                                                                             >  2:132718/1‑149 (MQ=255)
                                       tcCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGgcgc                                                                                                            >  3:20109/1‑149 (MQ=255)
                                       tcCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGgcgc                                                                                                            >  8:291834/1‑149 (MQ=255)
                                         cAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑Gt                                                                                                         >  5:41768/1‑149 (MQ=255)
                                         cAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑Gt                                                                                                         >  4:177860/1‑149 (MQ=255)
                                          aaGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCCGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCATGCCTGgcgcaca                                                                                                         <  3:268325/149‑4 (MQ=255)
                                          aaGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑Gtt                                                                                                        >  6:84523/1‑149 (MQ=255)
                                          aaGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGGCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑Gtt                                                                                                        >  4:217351/1‑149 (MQ=255)
                                          aaGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGCCATGTGCATGACCgggg                                                                                                                                      >  6:312922/1‑120 (MQ=255)
                                          aaGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGCCATGTGCATGACCgggg                                                                                                                                      <  5:312922/120‑1 (MQ=255)
                                            gATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATTCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTgc                                                                                                      <  5:288913/149‑1 (MQ=255)
                                            gATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCa                                                                                                                  <  3:70192/138‑1 (MQ=255)
                                            gATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCa                                                                                                                  >  4:70192/1‑138 (MQ=255)
                                                    gTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATgc                                                                                              >  1:137935/1‑149 (MQ=255)
                                                          gCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAAcc                                                                                        >  8:297202/1‑149 (MQ=255)
                                                           cGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCAGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAAc                                                                                         <  4:18785/147‑1 (MQ=255)
                                                                cgatcgatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGttttt                                                                                                                                <  8:248766/104‑1 (MQ=255)
                                                                cgatcgatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTg                                                                                  <  7:41636/149‑1 (MQ=255)
                                                                 gatcgatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTgg                                                                                 <  1:132718/149‑1 (MQ=255)
                                                                     gatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTtgt                                                                             <  3:51715/149‑1 (MQ=255)
                                                                     gatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTtgt                                                                             >  4:280858/1‑149 (MQ=255)
                                                                       tAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTgtgg                                                                           <  3:280858/149‑1 (MQ=255)
                                                                         aaCTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTgtggtg                                                                         <  5:257193/149‑1 (MQ=255)
                                                                         aaCTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCAAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTgtggtg                                                                         <  6:141049/149‑1 (MQ=255)
                                                                          aCTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGtggtgg                                                                        <  4:4526/149‑1 (MQ=255)
                                                                          aCTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGtggtgg                                                                        >  3:333821/1‑149 (MQ=255)
                                                                              ccTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCaa                                                                    >  7:333403/1‑149 (MQ=255)
                                                                              ccTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCaa                                                                    <  5:292912/149‑1 (MQ=255)
                                                                               cTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAg                                                                   >  4:279463/1‑149 (MQ=255)
                                                                                tCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGcc                                                                      <  5:339455/145‑1 (MQ=255)
                                                                                tCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGcc                                                                      >  6:339455/1‑145 (MQ=255)
                                                                                  agcagcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTgg                                                                >  3:151388/1‑149 (MQ=255)
                                                                                       cagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGAc                                                           >  2:147639/1‑149 (MQ=255)
                                                                                       cagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGAc                                                           >  7:343794/1‑149 (MQ=255)
                                                                                                 ggCGGTACCCTTATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAg                                                     <  2:326787/145‑1 (MQ=255)
                                                                                                 ggCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAg                                                     <  3:154505/145‑1 (MQ=255)
                                                                                                 ggCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAg                                                     <  2:321676/145‑1 (MQ=255)
                                                                                                   cGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaaa                                               <  5:84523/149‑1 (MQ=255)
                                                                                                    ggTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAg                                              <  3:177860/149‑1 (MQ=255)
                                                                                                     gTACCCTGATCGCGACGCTTGCGGCGTGATGCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTgg                                                       >  5:73856/1‑139 (MQ=255)
                                                                                                     gTACCCTGATCGCGACGCTTGCGGCGTGATGCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTgg                                                       <  6:73856/139‑1 (MQ=255)
                                                                                                        cccTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCct                                          >  7:329430/1‑149 (MQ=255)
                                                                                                         ccTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc                                         <  3:279463/149‑1 (MQ=255)
                                                                                                         ccTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc                                         >  2:343059/1‑149 (MQ=255)
                                                                                                         ccTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc                                         <  4:20109/149‑1 (MQ=255)
                                                                                                         ccTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc                                         <  8:333403/149‑1 (MQ=255)
                                                                                                          cTGATAGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc                                         >  1:44877/1‑148 (MQ=255)
                                                                                                           tGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGAGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTa                                       >  5:88473/1‑149 (MQ=255)
                                                                                                             aTCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTaca                                     <  2:44877/149‑1 (MQ=255)
                                                                                                             aTCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTa                                       <  5:47353/147‑1 (MQ=255)
                                                                                                                gcgACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTacacac                                  >  5:352688/1‑149 (MQ=255)
                                                                                                                       ttGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGc                             >  7:8272/1‑147 (MQ=255)
                                                                                                                       ttGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGc                             >  7:137368/1‑147 (MQ=255)
                                                                                                                         gcggcgTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCCGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCaa                           <  7:291834/147‑1 (MQ=255)
                                                                                                                          cggcgTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAgg                         <  3:48845/148‑1 (MQ=255)
                                                                                                                               tGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAgt                   <  7:57483/149‑1 (MQ=255)
                                                                                                                                 aTCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAgtgt                 <  7:297202/149‑1 (MQ=255)
                                                                                                                                   ccAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGGAAGGGGAAgtgt                 >  5:285217/1‑147 (MQ=255)
                                                                                                                                    cAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTa               <  2:137935/148‑1 (MQ=255)
                                                                                                                                       ggCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCg           >  6:376181/1‑149 (MQ=255)
                                                                                                                                       ggCCCCCGGTCATGCACATGGCCGGGGGTTTTCGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCg           >  6:96262/1‑149 (MQ=255)
                                                                                                                                         cccccGGCCATGTGCATGACCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCgg          >  2:184076/1‑149 (MQ=14)
                                                                                                                                          ccccGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAg        <  4:151388/149‑1 (MQ=255)
                                                                                                                                           cccGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAgg       >  4:49786/1‑149 (MQ=255)
                                                                                                                                           cccGGCCATGTGCATGACCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCa                                                                     >  7:254782/1‑87 (MQ=14)
                                                                                                                                           cccGGCCATGTGCATGACCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCa                                                                     <  8:254782/87‑1 (MQ=14)
                                                                                                                                              ggTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCaa    <  6:82561/149‑1 (MQ=255)
                                                                                                                                               gTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAg   >  6:198763/1‑149 (MQ=255)
                                                                                                                                                tCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAgg  <  4:333821/149‑1 (MQ=255)
                                                                                                                                                tCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAgg  <  4:14395/149‑1 (MQ=255)
                                                                                                                                                tCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAgg  <  3:136083/149‑1 (MQ=255)
                                                                                                                                                tCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAgg  <  3:130511/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
TGCACCTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGC‑GTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGG  >  NC_002947/1425506‑1425798

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: ‑