Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,425,663 A→G 8.8% intergenic (+58/‑447) 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,6630AG 8.8% 389.4 / 17.7 133intergenic (+58/‑447)PP_1249/PP_1250DUF4223 domain‑containing protein/group II intron‑encoding maturase
Reads supporting (aligned to +/- strand):  ref base A (67/54);  new base G (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.66e-01

Only 100 of 139 total aligned reads displayed.

ACCTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGC  >  NC_002947/1425516‑1425812
                                                                                                                                                   |                                                                                                                                                           
aCCTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGACCCAAGGCCCCCGGTCATGcacc                                                                                                                                                           >  2:1435305/1‑149 (MQ=255)
       aCTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCgg                                                                                                                                                    <  1:1552772/150‑1 (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)
               ttCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGttttt                                                                                                                                            >  1:640202/1‑150 (MQ=255)
                 ccACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTtgt                                                                                                                                          <  2:2109543/150‑1 (MQ=255)
                     ccTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑gc                                                                                                                                    <  1:1425138/150‑1 (MQ=255)
                     ccTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑gc                                                                                                                                    >  1:1833135/1‑150 (MQ=255)
                      cTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGCCATGTGCATGACCgggg                                                                                                                                                  >  1:2331629/1‑137 (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)
                                cATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATgg                                                                                                                      >  2:2012527/1‑150 (MQ=255)
                                    tcCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGcgcg                                                                                                                  >  2:1687468/1‑150 (MQ=255)
                                      cAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCgggg                                                                                                                                                  <  2:52159/121‑1 (MQ=255)
                                      cAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCgggg                                                                                                                                                  >  1:52159/1‑121 (MQ=255)
                                      cAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGt                                                                                                                                                <  1:906369/123‑1 (MQ=255)
                                      cAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGt                                                                                                                                                >  2:906369/1‑123 (MQ=255)
                                      cAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGtt                                                                                                                >  2:241758/1‑150 (MQ=255)
                                              tcGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTAt                                                                                                        <  2:2204469/150‑1 (MQ=255)
                                                  tGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTTATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCATACgcgc                                                                                                    >  2:2289726/1‑150 (MQ=255)
                                                   ggCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCa                                                                                                   >  1:208099/1‑150 (MQ=255)
                                                     cTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTTATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCATACGCGCAAc                                                                                                 <  1:2289726/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                                                                                    >  1:1392623/1‑150 (MQ=255)
                                                                   atAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCgggg                                                                                                                                                  >  2:63996/1‑92 (MQ=255)
                                                                   atAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtgg                                                                                   >  1:1700886/1‑150 (MQ=255)
                                                                    tAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTCGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtggt                                                                                  >  2:660153/1‑150 (MQ=255)
                                                                      aaCTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGtggtgg                                                                                >  2:1779805/1‑150 (MQ=255)
                                                                         tGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCa                                                                            >  2:2003536/1‑150 (MQ=255)
                                                                             tCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCATACGCGCAAcgc                                                                                               >  2:2079188/1‑126 (MQ=255)
                                                                             tCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCATACGCGCAAcgc                                                                                               <  1:2079188/128‑3 (MQ=255)
                                                                             tCAGCAAAAGTAATCTTGGCGGTACACTGATCGCGACGCTTGCGGCGTGAAACAAGGCCCCCGGTCATGCACATGACCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTg                                                                        >  2:1048264/1‑150 (MQ=25)
                                                                              cagcagcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTgg                                                                       <  2:1866574/150‑1 (MQ=255)
                                                                                  agcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGa                                                                   <  2:1394100/150‑1 (MQ=255)
                                                                                   gcagTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGAc                                                                  >  2:551030/1‑150 (MQ=255)
                                                                                       tAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGt                                                                         <  2:964250/139‑1 (MQ=255)
                                                                                       tAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGt                                                                         >  1:964250/1‑139 (MQ=255)
                                                                                       tAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTgg                                                              >  1:799156/1‑150 (MQ=255)
                                                                                         aTCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTtgtgt                                                                                                                                        >  1:152652/1‑80 (MQ=255)
                                                                                         aTCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTtgtgt                                                                                                                                        <  2:152652/80‑1 (MQ=255)
                                                                                            ttGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCAGATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTAACTTGGAGGTg                                                         >  2:1258220/1‑150 (MQ=255)
                                                                                            ttGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTACGCGTCTGCGAAAGTTACTTGTTTTGGTTTCCTAAGTGTGTGACTTGGAgctc                                                         >  2:624520/1‑147 (MQ=25)
                                                                                               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)
                                                                                                       cTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTa                                              <  1:1687468/150‑1 (MQ=255)
                                                                                                              cgACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCTCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACAccc                                       >  2:741714/1‑150 (MQ=255)
                                                                                                               gACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTAt                                                                                                        >  1:190925/1‑85 (MQ=255)
                                                                                                               gACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTAt                                                                                                        <  2:189924/85‑1 (MQ=255)
                                                                                                                  gCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTa                                                                                                         >  2:1349644/1‑81 (MQ=255)
                                                                                                                  gCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTa                                                                                                         <  1:1349644/81‑1 (MQ=255)
                                                                                                                    ttGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAg                                 >  1:658759/1‑150 (MQ=255)
                                                                                                                        ggcgTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAgggg                              <  1:1254189/149‑1 (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)
                                                                                                                          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                    >  1:800111/1‑150 (MQ=255)
                                                                                                                                   aGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTgc                                                                                                              >  1:2271812/1‑59 (MQ=255)
                                                                                                                                   aGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCg                  <  1:660153/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                                                                     <  2:2296103/96‑1 (MQ=14)
                                                                                                                                      cccccGGCCATGTGCATGACCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGt                                                                     >  1:2296103/1‑96 (MQ=14)
                                                                                                                                       ccccGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGGGGGCCGGAgg              >  1:118288/1‑150 (MQ=255)
                                                                                                                                       ccccGGCCATGTGCATGACCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg                                                                                    >  2:1907343/1‑81 (MQ=14)
                                                                                                                                       ccccGGCCATGTGCATGACCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg                                                                                    >  2:2063086/1‑81 (MQ=14)
                                                                                                                                       ccccGGCCATGTGCATGACCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg                                                                                    <  1:1907343/81‑1 (MQ=14)
                                                                                                                                        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)
                                                                                                                                                   acaTGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAg                                                            >  1:1190112/1‑92 (MQ=255)
                                                                                                                                                   acaTGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAg                                                            <  2:1190112/92‑1 (MQ=255)
                                                                                                                                                   acaTGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTcgc  >  2:808491/1‑150 (MQ=39)
                                                                                                                                                   acaTGGCCGGGGGTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGtggtgg                                                                                <  1:2057635/72‑1 (MQ=255)
                                                                                                                                                   acaTGGCCGGGGGTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGtggtgg                                                                                >  2:2057635/1‑72 (MQ=255)
                                                                                                                                                   |                                                                                                                                                           
ACCTACGACTACCTGTTCCACCCTTCGGTTTCCATCTCCAAGATCATCGGTGGCTGCGGCCCGATCGATAAACTGCCTCAGCAGCAGTAATCTTGGCGGTACCCTGATCGCGACGCTTGCGGCGTGATCCAAGGCCCCCGGTCATGCACATGGCCGGGGGTTTTTGTGT‑‑GCGC‑‑CAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGC  >  NC_002947/1425516‑1425812

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: