Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,425,657 T→C 8.2% 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,6570TC 8.2% 525.8 / 25.2 183intergenic (+52/‑453)PP_1249/PP_1250DUF4223 domain‑containing protein/group II intron‑encoding maturase
Reads supporting (aligned to +/- strand):  ref base T (83/85);  new base C (8/7);  total (91/92)
Fisher's exact test for biased strand distribution p-value = 7.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.08e-01

Only 100 of 184 total aligned reads displayed.

AAAACTGCACCTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGG  >  NC_002947/1425508‑1425805
                                                                                                                                                     |                                                                                                                                                    
   aCTGCACCTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGCCAt                                                                                                                                                   <  2:51911/150‑1 (MQ=255)
          cTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACAt                                                                                                                                             <  1:2144534/149‑1 (MQ=255)
          cTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATg                                                                                                                                            >  1:358395/1‑150 (MQ=255)
          cTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATg                                                                                                                                            >  1:357941/1‑150 (MQ=255)
            aCGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGc                                                                                                                                          <  1:2324700/150‑1 (MQ=255)
               aCTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGCCATGTGCATGACCgg                                                                                                                                       <  2:2128631/150‑1 (MQ=255)
                   ccTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATGTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGt                                                                                                                                   >  1:2062281/1‑150 (MQ=255)
                       ttCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGTGTCGCGTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGttttt                                                                                                                               <  1:219249/150‑1 (MQ=255)
                       ttCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGTACATGGCCGGGGGttttt                                                                                                                               <  2:1148486/150‑1 (MQ=255)
                       ttCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGCCATGTGCATGACCgggggccttg                                                                                                                               <  1:705108/150‑6 (MQ=255)
                        tCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTATTGCATAGTGCCTGGGGTTTTtg                                                                                                                              <  2:2282853/150‑1 (MQ=255)
                          cACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTtgtg                                                                                                                            >  2:1925599/1‑150 (MQ=255)
                           aCCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTtgtgt                                                                                                                           <  1:2940301/150‑1 (MQ=255)
                           aCCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTtgtgt                                                                                                                           >  2:2314411/1‑150 (MQ=255)
                           aCCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTtgtgt                                                                                                                           <  2:357941/150‑1 (MQ=255)
                           aCCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTtgtgt                                                                                                                           >  1:628337/1‑150 (MQ=255)
                           aCCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTtgt                                                                                                                             >  1:290542/1‑148 (MQ=255)
                             ccTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTgc                                                                                                                         >  1:1867830/1‑150 (MQ=255)
                               ttCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTc                                                                                                                                                     <  2:72453/120‑1 (MQ=255)
                                   gTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGcac                                                                                                                                               >  2:1994288/1‑122 (MQ=255)
                                       ccATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTAATAAGGAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATg                                                                                                               <  1:968036/150‑1 (MQ=255)
                                       ccATCTCCAAGATCATCGGTGGCTGCGGACCGATAGGTAAACTGCCTCAGCAGCAATAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACAGGGCCGGGGGTTTTTGTGTGCGCCAGGCATg                                                                                                               <  1:505734/150‑1 (MQ=255)
                                        cATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGCCATGTGCATGACCgggg                                                                                                                                     <  2:2938482/127‑1 (MQ=255)
                                         aTCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGgc                                                                                                             >  1:1127610/1‑150 (MQ=255)
                                              cAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTtgt                                                                                                                             >  2:1770089/1‑129 (MQ=255)
                                              cAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGtt                                                                                                        >  1:796180/1‑150 (MQ=255)
                                                     atcGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTa                                                                                                 >  1:2247148/1‑150 (MQ=255)
                                                        ggTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATgc                                                                                              <  1:510312/150‑1 (MQ=255)
                                                              tGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTgcgc                                                                                                    >  2:897727/1‑138 (MQ=255)
                                                                 ggCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGc                                                                                     <  1:2658849/150‑1 (MQ=255)
                                                                  gCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCAGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCATTCATTTCGCGTTGCGCGTATCCGCAAcaccgt                                                                                    >  2:503551/1‑145 (MQ=255)
                                                                     cgatcgatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTgg                                                                                 <  2:1867830/150‑1 (MQ=255)
                                                                      gatcgatAAACTGCCTCAGCAGCACTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGt                                                                                >  2:581472/1‑150 (MQ=255)
                                                                       atcgatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGtt                                                                               <  1:1261223/150‑1 (MQ=255)
                                                                        tcgatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATg                                                                                               <  2:1751771/133‑1 (MQ=255)
                                                                          gatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTtgt                                                                                                                             >  2:1541759/1‑101 (MQ=255)
                                                                          gatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTtgt                                                                                                                             <  1:1541759/101‑1 (MQ=255)
                                                                          gatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg                                                                            <  2:2643930/150‑1 (MQ=255)
                                                                          catAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg                                                                            <  1:2432431/149‑1 (MQ=255)
                                                                             aaacgaccTCAGCAGCGGTAATCTTGGCGGTGTCCTCATCGCGACGCTTCCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtggtg                                                                         <  2:300390/144‑1 (MQ=25)
                                                                              aaCTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGtggtgg                                                                        >  2:512893/1‑150 (MQ=255)
                                                                               aCTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAt                                                                                                <  2:1041765/125‑1 (MQ=255)
                                                                               aCTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTAt                                                                                                >  1:1041765/1‑125 (MQ=255)
                                                                                 tGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtggtg                                                                         >  1:1516201/1‑146 (MQ=255)
                                                                                  gCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCaa                                                                    <  2:1316816/150‑1 (MQ=255)
                                                                                    cTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGt                                                                  >  1:2610040/1‑150 (MQ=255)
                                                                                     tCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTgg                                                                                 >  1:1937141/1‑134 (MQ=255)
                                                                                      cagcagcagTAATCTTGGCGGTACCCTGATCGCGCCTGTGGCGTCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTgg                                                                <  1:2212184/150‑1 (MQ=25)
                                                                                      cagcagcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTgg                                                                <  1:2314411/150‑1 (MQ=255)
                                                                                       agcagcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTTGTTGTGGTGGCCAAGTggg                                                               <  2:1790550/150‑1 (MQ=255)
                                                                                       agcagcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTggg                                                               >  1:2107181/1‑150 (MQ=255)
                                                                                        gcagcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg                                                                            <  2:2628654/136‑1 (MQ=255)
                                                                                        gcagcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGt                                                              <  2:2107181/150‑1 (MQ=255)
                                                                                        gcagcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGt                                                              >  2:419804/1‑150 (MQ=255)
                                                                                         cagcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACCTGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTg                                                             >  1:574004/1‑150 (MQ=255)
                                                                                              gTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCtt                                                                                   <  2:2164889/123‑1 (MQ=255)
                                                                                                aaTCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCTAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGa                                                      <  1:1994289/150‑1 (MQ=255)
                                                                                                 aTCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACGTGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtggtg                                                                         >  1:2288699/1‑130 (MQ=255)
                                                                                                 aTCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAg                                                     <  1:879672/150‑1 (MQ=255)
                                                                                                  tCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAgg                                                    >  2:1447781/1‑150 (MQ=255)
                                                                                                    ttGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtggtg                                                                         <  2:1250181/127‑1 (MQ=255)
                                                                                                     tGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTGTGGGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGa                                                 >  1:1749917/1‑150 (MQ=255)
                                                                                                          gTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGGACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTc                                            <  2:2062281/150‑1 (MQ=255)
                                                                                                            aCCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCct                                          <  1:1867212/150‑1 (MQ=255)
                                                                                                            aCCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCGGCATGGCGCGTTGCGCGTATGCGCAACCCGCGTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc                                         <  2:2227705/150‑1 (MQ=255)
                                                                                                              ccTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCtt                                                                                   <  2:39344/107‑1 (MQ=255)
                                                                                                                     gcgACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACAcc                                 <  2:2434157/150‑1 (MQ=255)
                                                                                                                       gACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCaa                                                                    >  2:2830299/1‑113 (MQ=255)
                                                                                                                        aCGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCgg                              <  2:2486989/150‑1 (MQ=255)
                                                                                                                             tGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAgg                         <  1:1305973/150‑1 (MQ=255)
                                                                                                                              gcggcgTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTCTGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc                                         <  2:1985490/133‑1 (MQ=255)
                                                                                                                              gcggcgTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc                                         >  1:1985490/1‑133 (MQ=255)
                                                                                                                                 gcgTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg                                                                            <  2:869471/95‑1 (MQ=255)
                                                                                                                                  cgTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATTGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGc                             <  1:546278/141‑1 (MQ=255)
                                                                                                                                  cgTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATTGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGc                             >  2:546278/1‑141 (MQ=255)
                                                                                                                                   gTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGa                                                 <  1:2665541/120‑1 (MQ=255)
                                                                                                                                   gTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGa                                                 >  2:2665541/1‑120 (MQ=255)
                                                                                                                                      aTCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGaa                     >  1:149056/1‑145 (MQ=255)
                                                                                                                                       tCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTa               >  1:461519/1‑150 (MQ=255)
                                                                                                                                        ccAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGTGCAACCCGCTTGGTTGTGGTGGCCAAGTg                                                                 <  1:176473/99‑1 (MQ=255)
                                                                                                                                         aaaaacccccGGCCATGTGCATGACCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaaa                                               >  2:1511870/6‑116 (MQ=14)
                                                                                                                                             gCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg                                                                            <  2:199464/83‑1 (MQ=255)
                                                                                                                                             gCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg                                                                            >  1:199464/1‑83 (MQ=255)
                                                                                                                                             gCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGt                                                   <  1:1015860/108‑1 (MQ=255)
                                                                                                                                             gCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGa         <  2:1949772/150‑1 (MQ=255)
                                                                                                                                              cccccGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAg        <  2:2991724/150‑1 (MQ=255)
                                                                                                                                              cccccGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAg        <  2:1127610/150‑1 (MQ=255)
                                                                                                                                              cccccGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAg        <  1:1178066/150‑1 (MQ=255)
                                                                                                                                               ccccGGTCATGGACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTcc                                           <  2:341700/114‑1 (MQ=255)
                                                                                                                                               ccccGGTCATGGACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTcc                                           >  1:341700/1‑114 (MQ=255)
                                                                                                                                               ccccGGCCATGTGCATGACCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGa                                                            >  2:2204883/1‑97 (MQ=14)
                                                                                                                                               ccccGGCCATGTGCATGACCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAgg       >  1:3004724/1‑150 (MQ=255)
                                                                                                                                               ccccGGCCATGTGCATGACCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAgg       >  1:560342/1‑150 (MQ=11)
                                                                                                                                                cccGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc                                         >  2:1335634/1‑115 (MQ=255)
                                                                                                                                                cccGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc                                         <  1:1335634/115‑1 (MQ=255)
                                                                                                                                                cccGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAgtgt                 >  1:1890236/1‑139 (MQ=255)
                                                                                                                                                cccGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAgtgt                 <  2:1890236/139‑1 (MQ=255)
                                                                                                                                                cccGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGc      >  1:933022/1‑150 (MQ=255)
                                                                                                                                                  cGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCaa    <  1:400779/150‑1 (MQ=255)
                                                                                                                                                    gTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCACCCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAgg  <  1:2475034/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAAACTGCACCTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGG  >  NC_002947/1425508‑1425805

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