Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,450,603 T→C 6.1% K249K (AAA→AAG tynA ← tyramine oxidase, copper‑requiring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,450,6030TC 6.1% 432.0 / 13.7 133K249K (AAA→AAGtynAtyramine oxidase, copper‑requiring
Reads supporting (aligned to +/- strand):  ref base T (75/49);  new base C (4/4);  total (80/53)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.76e-01

Only 100 of 139 total aligned reads displayed.

TCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATGGGTTGCCAGGAGAGCAGTTTGTTGTTTTG  >  NC_000913/1450454‑1450752
                                                                                                                                                     |                                                                                                                                                     
tCAACGACCGCCACCAGTTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACt                                                                                                                                                       <  1:1718087/150‑1 (MQ=255)
 cAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACtt                                                                                                                                                      >  1:1540061/1‑150 (MQ=255)
  aaCGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACttt                                                                                                                                                     >  1:2224680/1‑150 (MQ=255)
     gACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACtttttt                                                                                                                                                  <  1:2366812/150‑1 (MQ=255)
     gACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACtttttt                                                                                                                                                  <  1:2305157/150‑1 (MQ=255)
      aCCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTc                                                                                                                                                 >  1:1055592/1‑150 (MQ=255)
      aCCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTc                                                                                                                                                 <  1:2331210/150‑1 (MQ=255)
      aCCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTc                                                                                                                                                 <  1:2293319/150‑1 (MQ=255)
         gCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACATTTTTCGCa                                                                                                                                              >  2:7945/1‑150 (MQ=255)
            accaGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCa                                                                                                                                           >  2:1542418/1‑150 (MQ=255)
            accaGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCa                                                                                                                                           <  1:416266/150‑1 (MQ=255)
            accaGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCa                                                                                                                                           <  2:2001056/150‑1 (MQ=255)
            accaGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCa                                                                                                                                           <  1:1597108/150‑1 (MQ=255)
               aGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAACCACTTTTTTCGCATCAGTa                                                                                                                                        >  1:1080007/1‑150 (MQ=255)
                ggTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACCTTTTTCGCATCAGTaa                                                                                                                                       <  2:2806803/150‑1 (MQ=255)
                ggTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACCTTTTTCGCATCAGTaa                                                                                                                                       <  2:1362958/150‑1 (MQ=255)
                ggTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACCTTTTTCGCATCAGTaa                                                                                                                                       <  2:2806798/150‑1 (MQ=255)
                    ttCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACCTTTTTCGCATCAGTAATAcc                                                                                                                                   >  1:897144/1‑150 (MQ=255)
                    ttCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACCTTTTTCGCATCAGTAATAcc                                                                                                                                   >  1:897247/1‑150 (MQ=255)
                      cGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACcgc                                                                                                                                 <  1:2341031/150‑1 (MQ=255)
                           ggATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTttc                                                                                                                            >  1:2153478/1‑150 (MQ=255)
                             aTGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTttctt                                                                                                                          <  1:1160789/150‑1 (MQ=255)
                              tgtgCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTttcttc                                                                                                                         >  2:2691239/1‑150 (MQ=255)
                                  cccAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACgg                                                                                                                     >  1:699489/1‑150 (MQ=255)
                                   ccAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGc                                                                                                                    >  1:2087418/1‑150 (MQ=255)
                                   ccAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGc                                                                                                                    >  2:606135/1‑150 (MQ=255)
                                    cAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCa                                                                                                                   >  2:2672118/1‑150 (MQ=255)
                                    cAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACCTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCa                                                                                                                   >  1:1251452/1‑150 (MQ=255)
                                       tagtTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCg                                                                                                                >  1:1005650/1‑150 (MQ=255)
                                                tCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAAttct                                                                                                       <  1:2111309/150‑1 (MQ=255)
                                                tCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAAttct                                                                                                       <  1:467707/150‑1 (MQ=255)
                                                tCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAAttct                                                                                                       >  1:2689202/1‑150 (MQ=255)
                                                        aTCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACtgtt                                                                                               >  1:656947/1‑150 (MQ=255)
                                                        aTCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACtgtt                                                                                               >  2:2284181/1‑150 (MQ=255)
                                                        aTCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACtgtt                                                                                               >  2:2147479/1‑150 (MQ=255)
                                                        aTCAAGATAGCTGATGACTATGAGCACCCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTCATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCTCtgtt                                                                                               >  2:2147478/1‑150 (MQ=255)
                                                          cAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACtgttgt                                                                                             <  2:1755103/150‑1 (MQ=255)
                                                             gATAGCTGATGACTTTGAGCAACCGGGCATCTTGATTCAGGCCAACTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTtaa                                                                                          <  2:2263241/150‑1 (MQ=255)
                                                               tAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTtaata                                                                                        <  1:1615144/150‑1 (MQ=255)
                                                                  cTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATg                                                                                     >  2:2502457/1‑150 (MQ=255)
                                                                  cTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATg                                                                                     >  2:2489252/1‑150 (MQ=255)
                                                                    gatgaCTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGtt                                                                                   <  2:342242/150‑1 (MQ=255)
                                                                    gatgaCTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGtt                                                                                   <  2:342250/150‑1 (MQ=255)
                                                                    gatgaCTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGtt                                                                                   <  2:448180/150‑1 (MQ=255)
                                                                      tgaCTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCt                                                                                 >  1:2855844/1‑150 (MQ=255)
                                                                      tgaCTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCt                                                                                 >  1:1186100/1‑150 (MQ=255)
                                                                      tgaCTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCt                                                                                 >  1:1759035/1‑150 (MQ=255)
                                                                         cTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGca                                                                              >  2:1734654/1‑150 (MQ=255)
                                                                          tttGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGcac                                                                             >  2:2868333/1‑150 (MQ=255)
                                                                            tGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCacac                                                                           <  2:47609/150‑1 (MQ=255)
                                                                              aGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTg                                                                         >  2:2434832/1‑150 (MQ=255)
                                                                              aGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTg                                                                         <  1:1952849/150‑1 (MQ=255)
                                                                               gCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTgg                                                                        >  1:697794/1‑150 (MQ=255)
                                                                                cAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGc                                                                       >  1:2004452/1‑150 (MQ=255)
                                                                                cAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGc                                                                       <  1:1802873/150‑1 (MQ=255)
                                                                                      ggCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAatc                                                                 >  2:1272569/1‑150 (MQ=255)
                                                                                       gCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAatca                                                                <  2:697794/150‑1 (MQ=255)
                                                                                         aTCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAatcatc                                                              >  1:2515231/1‑150 (MQ=255)
                                                                                            ttGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAg                                                           <  2:1817148/150‑1 (MQ=255)
                                                                                                   aGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACAcc                                                    <  1:1581719/150‑1 (MQ=255)
                                                                                                     gCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCAt                                                  >  2:490775/1‑150 (MQ=255)
                                                                                                      ccATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATa                                                 >  2:703062/1‑150 (MQ=255)
                                                                                                       cATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATAc                                                <  1:475860/150‑1 (MQ=255)
                                                                                                          cTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGt                                             <  2:2646188/150‑1 (MQ=255)
                                                                                                                cATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGt                                       <  2:1484563/150‑1 (MQ=255)
                                                                                                                      aaTAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTaa                                 >  1:182603/1‑150 (MQ=255)
                                                                                                                      aaTAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTaa                                 >  1:2581274/1‑150 (MQ=255)
                                                                                                                      aaTAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACGGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTaa                                 >  1:761069/1‑150 (MQ=255)
                                                                                                                      aaTAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAAAACCGCGTTTCTTCACGGCAGCGGCAAAATCTTCACTGTTGTTAATAATGTTATGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTat                                 >  1:761136/1‑149 (MQ=255)
                                                                                                                       aTAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCCGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAAt                                <  1:172575/150‑1 (MQ=255)
                                                                                                                        tAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATg                               <  2:2824869/150‑1 (MQ=255)
                                                                                                                          aCCTACGGTCAGCGGCGTGGTACTCCCTTTTTTCGCCTCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATCCCGTGGGCGTCTTTAAgggg                             <  2:2671250/150‑1 (MQ=255)
                                                                                                                          aCCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATggg                             <  2:2329145/150‑1 (MQ=255)
                                                                                                                          aCCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATggg                             <  2:579036/150‑1 (MQ=255)
                                                                                                                           ccTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATGGGt                            <  2:1005650/150‑1 (MQ=255)
                                                                                                                              aCGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATGGGTTGc                         >  1:1521713/1‑150 (MQ=255)
                                                                                                                              aCGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATGGGTTGc                         >  2:128175/1‑150 (MQ=255)
                                                                                                                              aCGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTAGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATGGGTTGc                         >  1:2568237/1‑150 (MQ=255)
                                                                                                                              aCGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAAAACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCAATGTTGTTAATAATGTTATGCACACTGGCGAAATCAACCAGCAACACCAAACCGAGGGCGACTATAAAGGGATGc                         >  1:100771/1‑150 (MQ=255)
                                                                                                                               cGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATGGGTTGcc                        <  1:1403246/150‑1 (MQ=255)
                                                                                                                                ggTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCAGTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACTCCATACCGTGGGCGTGTTTAATGGGTTGcct                       >  2:2413622/1‑149 (MQ=255)
                                                                                                                                ggTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATGGGTTGccc                       >  2:224029/1‑149 (MQ=255)
                                                                                                                                ggTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATGGGTTGccc                       >  2:1886580/1‑149 (MQ=255)
                                                                                                                                ggTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATGGGTTGCCa                       >  2:224030/1‑150 (MQ=255)
                                                                                                                                ggTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATGGGTTGCCa                       >  2:2167376/1‑150 (MQ=255)
                                                                                                                                 gTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATGGGTTGCCAg                      >  2:1252769/1‑150 (MQ=255)
                                                                                                                                     gcggcgTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATGGGTTGCCAGgaga                  <  2:2909793/150‑1 (MQ=255)
                                                                                                                                        gcgTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATGGGTTGCCAGGAGAGCa               >  1:964076/1‑150 (MQ=255)
                                                                                                                                           tGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATGGGTTGCCAGGAGAGCAGtt            <  2:1055592/150‑1 (MQ=255)
                                                                                                                                             gTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATGGGTTGCCAGGAGAGCAGTttg          <  2:935941/150‑1 (MQ=255)
                                                                                                                                             gTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATGGGTTGCCAGGAGAGCAGTttg          >  1:297372/1‑150 (MQ=255)
                                                                                                                                              tAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATGGGTTGCCAGGAGAGCAGTttgt         >  2:1296787/1‑150 (MQ=255)
                                                                                                                                              tAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATGGGTTGCCAGGAGAGCAGTttgt         >  2:2137997/1‑150 (MQ=255)
                                                                                                                                               aaTCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATGGGTTGCCAGGAGAGCAGTttgtt        <  2:1595951/150‑1 (MQ=255)
                                                                                                                                               aaTCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATGGGTTGCCAGGAGAGCAGTttgtt        <  1:2434832/150‑1 (MQ=255)
                                                                                                                                                 tCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATGGGTTGCCAGGAGAGCTGTTtgttgt      >  2:759602/1‑150 (MQ=255)
                                                                                                                                                 tCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATGGGTTGCCAGGAGAGCAGTTtgttgt      >  2:760094/1‑150 (MQ=255)
                                                                                                                                                    cTTTTTTTGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATGGGTTGCCAGGAGAGCAGTTTGTTGtttt   >  2:2174849/1‑150 (MQ=255)
                                                                                                                                                    cTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATGGGTTGCCAGGAGAGCAGTTTGTTGtttt   >  2:970143/1‑150 (MQ=255)
                                                                                                                                                     ttttttCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATGGGTTGCCAGGAGAGCAGTTTGTTGTTTTg  <  2:2496026/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCAACGACCGCCACCAGGTTTTCGATGGGATGTGCCCAGTAGTTGCCATCACCGACATCAAGATAGCTGATGACTTTGAGCAACCGGGCATCTTGTTTCAGGCCATCTTTACCATCGAAATAACCTACGGTCAGCGGCGTGGTAATCACTTTTTTCGCATCAGTAATACCGCGTTTCTTCACGGCAGCGGCAAATTCTTCACTGTTGTTAATAATGTTCTGCACACTGGCGAAATCATCCAGCAACACCATACCGTGGGCGTCTTTAATGGGTTGCCAGGAGAGCAGTTTGTTGTTTTG  >  NC_000913/1450454‑1450752

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: