Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,012,381 G→A 5.2% Q562* (CAG→TAG)  PP_1793 ← group 2 family glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,012,3810GA 5.2% 600.6 / 13.5 190Q562* (CAG→TAG) PP_1793group 2 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (94/86);  new base A (7/3);  total (101/89)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 190 total aligned reads displayed.

GATGGCTATCCCGGTGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATCCCGCCTTGATGC  >  NC_002947/2012234‑2012529
                                                                                                                                                   |                                                                                                                                                    
gATGGCTATCCCGGTGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGCATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTa                                                                                                                                                    >  1:3095097/1‑150 (MQ=255)
 aTGGCTATCCCGGTGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAg                                                                                                                                                   <  1:1404753/150‑1 (MQ=255)
         cccGGTGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGACGCGCTACAGCGCGcccg                                                                                                                                           <  1:1454477/150‑2 (MQ=255)
          ccGGTGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCg                                                                                                                                          >  2:4055079/1‑150 (MQ=255)
            ggTGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCggg                                                                                                                                        >  1:1487023/1‑150 (MQ=255)
            ggTGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCggg                                                                                                                                        >  1:1565154/1‑150 (MQ=255)
            ggTGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCggg                                                                                                                                        >  1:1565225/1‑150 (MQ=255)
              tGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCa                                                                                                                                      <  2:5082061/150‑1 (MQ=255)
                 cAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGAGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTa                                                                                                                                                    <  2:2521061/133‑1 (MQ=255)
                  aGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCtt                                                                                                                                  >  1:1418388/1‑150 (MQ=255)
                   gACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCttt                                                                                                                                 <  2:2544938/150‑1 (MQ=255)
                   gACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCttt                                                                                                                                 <  2:1470017/150‑1 (MQ=255)
                   gACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCttt                                                                                                                                 >  2:3310867/1‑150 (MQ=255)
                      aaGAGAGAGGGTAGAATCAGCCCCTGAGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATAGCTCTGTTGTGAGAAGCGGGACAGTCGGGCGCGCTGTAGCGCGTCACGGGCAGCTTTGAt                                                                                                                              >  1:4827817/1‑150 (MQ=255)
                        gagagagGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCg                                                                                                                            <  1:4953326/150‑1 (MQ=255)
                         agagagGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCgg                                                                                                                           <  2:543611/150‑1 (MQ=255)
                           agagGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGt                                                                                                                         >  1:3980042/1‑150 (MQ=255)
                           agagGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGt                                                                                                                         >  1:1215393/1‑150 (MQ=255)
                               ggTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCg                                                                                                                     >  2:3017997/1‑150 (MQ=255)
                                gTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTAATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGc                                                                                                                    >  1:3484294/1‑150 (MQ=255)
                                 tAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCa                                                                                                                   <  1:3265070/150‑1 (MQ=255)
                                   gAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAgg                                                                                                                 <  2:2389051/150‑1 (MQ=255)
                                     aTCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCa                                                                                                               >  1:5033311/1‑150 (MQ=255)
                                        aGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCaaaa                                                                                                            >  2:5009816/1‑150 (MQ=255)
                                            ccTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCa                                                                                                        <  2:1907342/150‑1 (MQ=255)
                                                 gCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACtata                                                                                                   <  1:3498392/150‑1 (MQ=255)
                                                  ccGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAg                                                                                                  >  1:1043734/1‑150 (MQ=255)
                                                  ccGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAg                                                                                                  <  2:2867389/150‑1 (MQ=255)
                                                  ccGTCATGAATGTGTAGCGTATTACCTGCCTCAATCCGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCCCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCAcgaaag                                                                                                  >  1:1043418/1‑145 (MQ=255)
                                                   cGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGc                                                                                                 <  2:5414475/150‑1 (MQ=255)
                                                     tCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGGTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGt                                                                                               >  2:4502040/1‑150 (MQ=255)
                                                     tCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTGTGTTGTGCGAAGCGGGCCTGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGt                                                                                               >  2:4501983/1‑150 (MQ=255)
                                                       aTGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGc                                                                                             <  1:4531424/150‑1 (MQ=255)
                                                        tGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGcc                                                                                            <  1:23276/150‑1 (MQ=255)
                                                        tGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGcc                                                                                            >  1:4265497/1‑150 (MQ=255)
                                                             gtgtAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAgc                                                                                       <  2:4682472/150‑1 (MQ=255)
                                                              tgtAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAgcg                                                                                      >  1:2221839/1‑150 (MQ=255)
                                                                tAGCGTTTTACCTGCCTCCATCGGCCTTTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTTCCGAAGCGcc                                                                                    <  2:1043418/150‑1 (MQ=255)
                                                                 aGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCt                                                                                   <  2:1404296/150‑1 (MQ=255)
                                                                 aGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCt                                                                                   <  2:1201759/150‑1 (MQ=255)
                                                                  gCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTg                                                                                  <  1:185647/150‑1 (MQ=255)
                                                                  gCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTg                                                                                  >  2:3428295/1‑150 (MQ=255)
                                                                         aCCTGCCTCCGTCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAgcg                                                                           >  1:3511406/1‑150 (MQ=255)
                                                                             gCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAgcggcgg                                                                       >  1:2708799/1‑150 (MQ=255)
                                                                             gCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAgcggcgg                                                                       >  1:1299293/1‑150 (MQ=255)
                                                                                tCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACtata                                                                                                   >  1:3774833/1‑119 (MQ=255)
                                                                                tCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGgcgc                                                                    >  2:1216845/1‑150 (MQ=255)
                                                                                 ccctCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCa                                                                   <  2:5299726/147‑1 (MQ=255)
                                                                                   aTCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAg                                                                 >  1:3493391/1‑150 (MQ=255)
                                                                                      ggCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATa                                                              <  1:4737910/150‑1 (MQ=255)
                                                                                        ccATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGc                                                            >  1:706818/1‑150 (MQ=255)
                                                                                          aTTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCgt                                                          >  1:2613334/1‑150 (MQ=255)
                                                                                          aTTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCgt                                                          >  1:2841272/1‑150 (MQ=255)
                                                                                           ttGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCgtt                                                         <  2:2123554/150‑1 (MQ=255)
                                                                                           ttGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCgtt                                                         <  2:940770/150‑1 (MQ=255)
                                                                                           ttGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCgtt                                                         <  2:4859237/150‑1 (MQ=255)
                                                                                              cGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCgttgtt                                                      <  1:4352557/150‑1 (MQ=255)
                                                                                              cGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCgttgtt                                                      <  1:1734883/150‑1 (MQ=255)
                                                                                               gAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAgc                                                                                       >  1:626728/1‑116 (MQ=255)
                                                                                                aGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTtgttgt                                                    <  1:3691755/150‑1 (MQ=255)
                                                                                                     ttttGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCg                                               >  1:4117247/1‑150 (MQ=255)
                                                                                                     ttttGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCg                                               >  2:2723868/1‑150 (MQ=255)
                                                                                                         gTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGAcgc                                           <  2:2293577/150‑1 (MQ=255)
                                                                                                          ttCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGAcgcg                                          <  2:4884374/150‑1 (MQ=255)
                                                                                                            cATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCa                                        <  1:28768/150‑1 (MQ=255)
                                                                                                             aTCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAg                                       >  2:3679637/1‑150 (MQ=255)
                                                                                                             aTCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAg                                       >  1:5346779/1‑150 (MQ=255)
                                                                                                             aTCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAg                                       >  2:5149198/1‑150 (MQ=255)
                                                                                                              tCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAgt                                      <  1:2088283/150‑1 (MQ=255)
                                                                                                               cGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAgtg                                     <  1:944329/150‑1 (MQ=255)
                                                                                                                gCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAgtgt                                    <  2:111052/150‑1 (MQ=255)
                                                                                                                  tctGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGttt                                  >  2:960533/1‑150 (MQ=255)
                                                                                                                      ttgtGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTg                                                                                  <  2:3182409/98‑1 (MQ=255)
                                                                                                                      ttgtGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTc                              <  1:3017997/150‑1 (MQ=255)
                                                                                                                       tgtgCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCTAAGCGACTTGATAGCGGCGGCGCAAGATATCGTTGTTGTGTTCGACGCGCAGTGTTTTTTcc                             <  2:721982/150‑1 (MQ=255)
                                                                                                                         tgCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTgcacaacagagcgatgaacaaaaagcctcgcaatggccgatggaggcaggtaatacgctacacattcatgacggcacaggggctgattctaccctctctctt                           <  2:1298740‑M2/150‑103 (MQ=255)
                                                                                                                         tgCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATGGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGc                           >  2:5183740/1‑150 (MQ=255)
                                                                                                                         tgCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATGGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGc                           >  2:5183789/1‑150 (MQ=255)
                                                                                                                          gCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGcc                          <  1:4463327/150‑1 (MQ=255)
                                                                                                                           cGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTgcacaacagagcgatgaacaaaaagcctcgcaatggccgatggaggcaggtaatacgctacacattcatgacggcacaggggctgattctaccctctctcttgt                         <  1:1510022‑M2/150‑105 (MQ=255)
                                                                                                                           cGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGccc                         <  2:1785085/150‑1 (MQ=255)
                                                                                                                           cGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGccc                         >  2:5421762/1‑150 (MQ=255)
                                                                                                                           cGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGccc                         <  1:2344681/150‑1 (MQ=255)
                                                                                                                               gCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCgaa                     >  1:3171725/1‑150 (MQ=255)
                                                                                                                                cGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCgaaa                    <  2:4120510/150‑1 (MQ=255)
                                                                                                                                  ggCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCgaaaga                  >  2:2277490/1‑150 (MQ=255)
                                                                                                                                   gCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCgaaagaa                 >  1:1355664/1‑150 (MQ=255)
                                                                                                                                      cccgggacgcgcTACAGCGCGCCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATc              >  2:1367234/8‑150 (MQ=255)
                                                                                                                                       gTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATcc             <  1:2462126/150‑1 (MQ=255)
                                                                                                                                       gTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATcc             >  2:1983988/1‑150 (MQ=255)
                                                                                                                                        cgggacgcgcTACAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATccc            >  1:2733066/6‑150 (MQ=255)
                                                                                                                                        cgggacgcgcTACAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATccc            >  1:5261556/6‑150 (MQ=255)
                                                                                                                                          ggacgcgcTACAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATCGCGGCGGCGCAAGATAGCGTGGTTGGGTGCGCCGCGCCGTGTTTTTTCCGCCCCGAAAGAAATCCCGc          >  2:1833172/4‑150 (MQ=255)
                                                                                                                                             cgcgcTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATCCCGCCtt       <  1:986874/150‑1 (MQ=255)
                                                                                                                                               cgcTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGTCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTTCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATCCCGCCTTGa     <  1:3202437/150‑1 (MQ=255)
                                                                                                                                                gcTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATCCCGCCTTGAt    <  1:3903359/150‑1 (MQ=255)
                                                                                                                                                gcTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATCCCGCCTTGAt    <  2:615081/150‑1 (MQ=255)
                                                                                                                                                gcTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATCCCGCCTTGAt    >  2:3371078/1‑150 (MQ=255)
                                                                                                                                                  tGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATCCCGCCTTGATgc  <  2:3980042/150‑1 (MQ=255)
                                                                                                                                                  tGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATCCCGCCTTGATgc  >  2:3228929/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GATGGCTATCCCGGTGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATCCCGCCTTGATGC  >  NC_002947/2012234‑2012529

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