Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,425,662 C→T 9.0% intergenic (+57/‑448) 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,6620CT 9.0% 373.3 / 19.8 133intergenic (+57/‑448)PP_1249/PP_1250DUF4223 domain‑containing protein/group II intron‑encoding maturase
Reads supporting (aligned to +/- strand):  ref base C (67/54);  new base T (6/6);  total (73/60)
Fisher's exact test for biased strand distribution p-value = 7.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 134 total aligned reads displayed.

ACCTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTG  >  NC_002947/1425516‑1425808
                                                                                                                                                  |                                                                                                                                                       
aCCTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGACCCAAGGCCCCCGGTCATGcacc                                                                                                                                                      >  2:1435305/1‑149 (MQ=255)
       aCTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCgg                                                                                                                                               <  1:351318/150‑1 (MQ=255)
         tACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCgggg                                                                                                                                             >  2:1040148/1‑150 (MQ=255)
         tACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCgggg                                                                                                                                             >  1:1413739/1‑150 (MQ=255)
         tACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGCCATGTGCATGACCgggg                                                                                                                                             <  1:234333/150‑1 (MQ=255)
          aCCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCggggg                                                                                                                                            >  2:1399723/1‑150 (MQ=255)
            cTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGtt                                                                                                                                          >  2:1043882/1‑150 (MQ=255)
               ttCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGttttt                                                                                                                                       >  1:2155614/1‑150 (MQ=255)
                 ccACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTtgt                                                                                                                                     <  2:2109543/150‑1 (MQ=255)
                     ccTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑gc                                                                                                                               >  1:1833135/1‑150 (MQ=255)
                      cTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGCCATGTGCATGACCgggg                                                                                                                                             <  2:2331629/137‑1 (MQ=255)
                            tttCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGttttt                                                                                                                                       <  1:2194069/137‑1 (MQ=255)
                            tttCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGttttt                                                                                                                                       >  2:2194069/1‑137 (MQ=255)
                               ccATCTCCAAGATCATCGGTGGTTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGCCATGTGCATGACCggg                                                                                                                                              >  1:201744/1‑127 (MQ=255)
                               ccATCTCCAAGATCATCGGTGGTTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGCCATGTGCATGACCggg                                                                                                                                              <  2:201744/127‑1 (MQ=255)
                                cATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATgg                                                                                                                 >  2:2012527/1‑150 (MQ=255)
                                    tcCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CCGGCATGGcgcg                                                                                                             >  1:772033/1‑150 (MQ=255)
                                    tcCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGcgcg                                                                                                             >  2:1687468/1‑150 (MQ=255)
                                    tcCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGcgcg                                                                                                             >  1:1328264/1‑150 (MQ=255)
                                      cAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCgggg                                                                                                                                             >  1:52159/1‑121 (MQ=255)
                                      cAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGt                                                                                                                                           >  2:906369/1‑123 (MQ=255)
                                              tcGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTAt                                                                                                   <  2:2204469/150‑1 (MQ=255)
                                                  tGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTTATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCATACgcgc                                                                                               >  2:2289726/1‑150 (MQ=255)
                                                     cTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTTATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCATACGCGCAAc                                                                                            <  1:2289726/150‑1 (MQ=255)
                                                     cTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAAc                                                                                            <  1:1399723/150‑1 (MQ=255)
                                                     cTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCATACGCGCAAc                                                                                            <  1:2354189/150‑1 (MQ=255)
                                                             cgatcgatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGATGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGCGGTTTTGGGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTgg                                                                                    <  1:436670/150‑1 (MQ=255)
                                                             cgatcgatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTgg                                                                                    >  1:1866574/1‑150 (MQ=255)
                                                               atcgatAAACTGCCTCAGCAGCCGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGGCGGGGGGTTTTGTGT‑‑GCGC‑‑CGGGCATGGCGCGTTGCGCGGATGCGCAACCCGCGTGGtt                                                                                  >  2:308089/1‑150 (MQ=38)
                                                               atcgatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGtt                                                                                  <  2:885920/150‑1 (MQ=255)
                                                                  gatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg                                                                               <  2:208099/150‑1 (MQ=255)
                                                                  gatAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg                                                                               >  1:1392623/1‑150 (MQ=255)
                                                                   atAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCgggg                                                                                                                                             >  2:63996/1‑92 (MQ=255)
                                                                   atAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtgg                                                                              >  1:1700886/1‑150 (MQ=255)
                                                                   atAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAATCCCCCGGTCATGCACATGGCCgggg                                                                                                                                             <  1:63996/92‑1 (MQ=255)
                                                                    tAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTCGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtggt                                                                             >  2:660153/1‑150 (MQ=255)
                                                                     aaaCTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtggtg                                                                            <  1:2285057/150‑1 (MQ=255)
                                                                     aaaCTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtggtg                                                                            <  2:1055322/150‑1 (MQ=255)
                                                                      aaCTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGtggtgg                                                                           >  2:1779805/1‑150 (MQ=255)
                                                                          gCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCaa                                                                       <  2:356507/150‑1 (MQ=255)
                                                                             tCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCATACGCGCAAcgc                                                                                          <  1:2079188/128‑3 (MQ=255)
                                                                                  agcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGa                                                               <  2:1394100/150‑1 (MQ=255)
                                                                                   gcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGAc                                                              >  2:551030/1‑150 (MQ=255)
                                                                                    cagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGATTGGTTGTGGTGGCCAAGTGGGTGACt                                                             >  2:2372920/1‑150 (MQ=255)
                                                                                      gTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGTCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTg                                                           <  1:882072/150‑1 (MQ=255)
                                                                                      gTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTg                                                           >  1:2340695/1‑150 (MQ=255)
                                                                                       tAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGt                                                                     <  2:964250/139‑1 (MQ=255)
                                                                                       tAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGt                                                                     >  1:964250/1‑139 (MQ=255)
                                                                                         aTCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTtgtgt                                                                                                                                   <  2:152652/80‑1 (MQ=255)
                                                                                         aTCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTtgtgt                                                                                                                                   >  1:152652/1‑80 (MQ=255)
                                                                                         aTCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTCCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAg                                                        <  2:1328264/150‑1 (MQ=255)
                                                                                               gCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGaaa                                                  <  2:1413739/150‑1 (MQ=255)
                                                                                                  gTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCAGGCACATGCCCGGGGACCTTTGTGTACGCGCAACCGGCAATGGCC‑‑‑ACGCAAGGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGCCACAGGGGTGAAAGTc                                               >  1:544293/1‑150 (MQ=12)
                                                                                                    aCCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCct                                             <  1:1043882/150‑1 (MQ=255)
                                                                                                     cccTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGc                                                                          >  1:1054311/1‑120 (MQ=255)
                                                                                                     cccTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGc                                                                          <  2:1054311/120‑1 (MQ=255)
                                                                                                     cccTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctc                                            <  1:1278093/150‑1 (MQ=255)
                                                                                                      ccTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCggg                                                                                                                                              <  2:211558/56‑1 (MQ=255)
                                                                                                      ccTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCggg                                                                                                                                              >  1:211558/1‑56 (MQ=255)
                                                                                                      ccTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctct                                           >  1:468330/1‑150 (MQ=255)
                                                                                                      ccTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctct                                           >  1:1674285/1‑150 (MQ=255)
                                                                                                      ccTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTGGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCctct                                           >  1:1203041/1‑150 (MQ=255)
                                                                                                       cTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTa                                          <  1:1687468/150‑1 (MQ=255)
                                                                                                            cgcgACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTacacac                                     <  2:1833135/150‑1 (MQ=255)
                                                                                                               gACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTAt                                                                                                   <  2:190925/85‑1 (MQ=255)
                                                                                                               gACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTAt                                                                                                   <  2:189924/85‑1 (MQ=255)
                                                                                                               gACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTAt                                                                                                   >  1:190925/1‑85 (MQ=255)
                                                                                                                  gCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTa                                                                                                    >  2:1349644/1‑81 (MQ=255)
                                                                                                                    ttGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAg                             >  1:658759/1‑150 (MQ=255)
                                                                                                                        ggcgTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAgggg                          >  2:1254189/1‑149 (MQ=255)
                                                                                                                        ggcgTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGa                         <  2:2340695/150‑1 (MQ=255)
                                                                                                                         gcgTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGaa                        <  1:2003536/150‑1 (MQ=255)
                                                                                                                         gcgTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGaa                        <  2:544293/150‑1 (MQ=255)
                                                                                                                          cgTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGT‑‑CAGGCATGGCGCGTTGCGCGTATGCGCACCCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAg                       <  1:417628/150‑1 (MQ=255)
                                                                                                                          cgTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTacacac                                     <  2:1477334/136‑1 (MQ=255)
                                                                                                                          cgTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTacacac                                     >  1:1477334/1‑136 (MQ=255)
                                                                                                                             gATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAAAATGTGTGGGTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAgtgt                    >  1:1860661/1‑150 (MQ=17)
                                                                                                                              aTCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGtt                   >  2:780052/1‑150 (MQ=255)
                                                                                                                               tCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTa                  >  1:1063828/1‑150 (MQ=255)
                                                                                                                                 cAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGc                <  2:640202/150‑1 (MQ=255)
                                                                                                                                 cAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGc                >  1:800111/1‑150 (MQ=255)
                                                                                                                                   aGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCg              <  1:660153/150‑1 (MQ=255)
                                                                                                                                   aGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCg              <  2:772033/150‑1 (MQ=255)
                                                                                                                                     gCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGa            >  2:2295992/1‑150 (MQ=255)
                                                                                                                                      cccccGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAg           <  1:1080709/150‑1 (MQ=255)
                                                                                                                                      cccccGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAg           >  1:317971/1‑150 (MQ=255)
                                                                                                                                      cccccGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAg           <  2:317971/150‑1 (MQ=255)
                                                                                                                                      cccccGGCCATGTGCATGACCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGt                                                                 >  1:2296103/1‑96 (MQ=14)
                                                                                                                                      cccccGGCCATGTGCATGACCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGt                                                                 <  2:2296103/96‑1 (MQ=14)
                                                                                                                                       ccccGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGGGGGCCGGAgg          >  1:118288/1‑150 (MQ=255)
                                                                                                                                       ccccGGCCATGTGCATGACCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg                                                                               <  1:1907343/81‑1 (MQ=14)
                                                                                                                                       ccccGGCCATGTGCATGACCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg                                                                               >  2:2063086/1‑81 (MQ=14)
                                                                                                                                       ccccGGCCATGTGCATGACCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg                                                                               <  1:2063086/81‑1 (MQ=14)
                                                                                                                                       ccccGGCCATGTGCATGACCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg                                                                               >  2:1907343/1‑81 (MQ=14)
                                                                                                                                       ccccGGCCATGTGCATGACCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAgg          >  1:1435225/1‑150 (MQ=11)
                                                                                                                                        cccGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGa                                                         <  1:1339134/102‑1 (MQ=255)
                                                                                                                                        cccGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGa                                                         >  2:1339134/1‑102 (MQ=255)
                                                                                                                                        cccGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGc         >  1:2400270/1‑150 (MQ=255)
                                                                                                                                          cGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCaa       >  2:2258340/1‑150 (MQ=255)
                                                                                                                                               aTGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGgtg  >  2:531006/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                       
ACCTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTG  >  NC_002947/1425516‑1425808

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: