Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,932,096 C→G 5.9% intergenic (+37/+19) purT → / ← eda phosphoribosylglycinamide formyltransferase 2/KHG/KDPG aldolase; 2‑dehydro‑3‑deoxy‑phosphogluconate/4‑hydroxy‑2‑ oxoglutarate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,932,0960CG 5.9% 733.2 / 27.8 202intergenic (+37/+19)purT/edaphosphoribosylglycinamide formyltransferase 2/KHG/KDPG aldolase; 2‑dehydro‑3‑deoxy‑phosphogluconate/4‑hydroxy‑2‑ oxoglutarate aldolase
Reads supporting (aligned to +/- strand):  ref base C (113/76);  new base G (7/5);  total (121/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 203 total aligned reads displayed.

CGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAG  >  NC_000913/1931947‑1932244
                                                                                                                                                      |                                                                                                                                                     
      ttGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCgg                                                                                                                                                 <  1:3180928/150‑1 (MQ=255)
       tGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCGAGGATCggg                                                                                                                                                <  1:2268936/150‑1 (MQ=255)
       tGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCGAGGATCggg                                                                                                                                                <  1:2435195/150‑1 (MQ=255)
           ggCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCA‑tt                                                                                                                                           <  2:10047/150‑1 (MQ=255)
           ggCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑tt                                                                                                                                           >  2:1728883/1‑150 (MQ=255)
            gCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑ttt                                                                                                                                          <  1:2572218/150‑1 (MQ=255)
            gCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑ttt                                                                                                                                          >  1:843865/1‑150 (MQ=255)
                ccGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGAc                                                                                                                                      <  2:227785/150‑1 (MQ=255)
                      tcTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTAc                                                                                                                                >  2:1151939/1‑150 (MQ=255)
                       cTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AACCAAACGCTCCCACAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACa                                                                                                                               <  2:2761977/150‑1 (MQ=255)
                       cTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACa                                                                                                                               <  2:2174868/150‑1 (MQ=255)
                         gggggTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGc                                                                                                                             <  2:1400869/150‑1 (MQ=255)
                          ggggTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCt                                                                                                                            >  1:1727893/1‑150 (MQ=255)
                                gCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGACCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAgcgc                                                                                                                      >  2:12832/1‑150 (MQ=255)
                                 cACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGcc                                                                                                                     >  2:270960/1‑150 (MQ=255)
                                       cTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTAc                                                                                                               >  2:819430/1‑150 (MQ=255)
                                        tACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACa                                                                                                              >  1:2967194/1‑150 (MQ=255)
                                        tACTGCAGAGAGTGTTGTTGACGCCATTGAAAGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATTCCCGATCCTCGAACGGGCA‑TGTTGACTTTTACATCTTAGCGCCTTCTACa                                                                                                              >  2:2235892/1‑150 (MQ=255)
                                         aCTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAg                                                                                                             <  2:2879824/150‑1 (MQ=255)
                                           tGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACCGcg                                                                                                           >  2:2543947/1‑149 (MQ=255)
                                           tGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCt                                                                                                           >  2:2543963/1‑150 (MQ=255)
                                             cAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTc                                                                                                         >  2:850356/1‑150 (MQ=255)
                                               gagagTGTTGTTGACGCCATTGAACGCGCGCAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCAc                                                                                                       <  1:2150705/150‑1 (MQ=255)
                                                  agTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCAcgcg                                                                                                    >  2:2487388/1‑150 (MQ=255)
                                                   gtgtTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACgcgc                                                                                                   >  1:1175204/1‑150 (MQ=255)
                                                     gttgttGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCa                                                                                                 <  1:746095/150‑1 (MQ=255)
                                                        gttgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCt                                                                                              >  1:3191068/1‑150 (MQ=255)
                                                          tgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTa                                                                                            <  2:2102720/150‑1 (MQ=255)
                                                           gACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAg                                                                                           >  2:2314686/1‑150 (MQ=255)
                                                            aCGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGt                                                                                          >  1:364357/1‑150 (MQ=255)
                                                             cGCCATTGAACGCGCGAAGCACGCCGCCGGACTGGTATAAGTACAGGGTTTAACCCGGC‑AAAATAACGCTTCAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTa                                                                                         <  1:3086570/150‑1 (MQ=255)
                                                             cGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTa                                                                                         <  1:2621650/150‑1 (MQ=255)
                                                              gCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTaa                                                                                        >  1:2826221/1‑150 (MQ=255)
                                                                  ttGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCg                                                                                    >  1:124909/1‑150 (MQ=255)
                                                                   tGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAACACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCgg                                                                                   <  2:2491282/150‑1 (MQ=255)
                                                                   tGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCgg                                                                                   <  2:1231712/150‑1 (MQ=255)
                                                                     aaCGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑GAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTc                                                                                 >  1:554401/1‑150 (MQ=255)
                                                                     aaCGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTc                                                                                 <  1:3014177/150‑1 (MQ=255)
                                                                     aaCGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTc                                                                                 <  1:1559555/150‑1 (MQ=255)
                                                                        gcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTa                                                                              >  2:2510862/1‑150 (MQ=255)
                                                                        gcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTa                                                                              >  1:573760/1‑150 (MQ=255)
                                                                        gcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCATa                                                                              >  1:2457543/1‑150 (MQ=255)
                                                                        gcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTGAGTAATGCGGTCGTa                                                                              >  1:2343905/1‑150 (MQ=255)
                                                                         cgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTaa                                                                             <  1:725264/150‑1 (MQ=255)
                                                                          gcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAAt                                                                            <  2:1307246/150‑1 (MQ=255)
                                                                            gAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCg                                                                          >  2:807266/1‑150 (MQ=255)
                                                                            gAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCg                                                                          >  1:1018122/1‑150 (MQ=255)
                                                                              aGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGcc                                                                        >  2:1554129/1‑150 (MQ=255)
                                                                                   gccgccGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCtt                                                                   >  1:1187644/1‑150 (MQ=255)
                                                                                      gccgGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCa                                                                <  2:3134602/150‑1 (MQ=255)
                                                                                      gccgGACAGGAAAAACTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TATTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCa                                                                <  1:3219029/150‑1 (MQ=255)
                                                                                       ccgGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAg                                                               <  1:1320106/150‑1 (MQ=255)
                                                                                       ccgGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAg                                                               >  2:2468848/1‑150 (MQ=255)
                                                                                       ccgGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAg                                                               <  1:2744573/150‑1 (MQ=255)
                                                                                       ccgGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAg                                                               <  2:3076440/150‑1 (MQ=255)
                                                                                       ccgGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAg                                                               <  1:1873226/150‑1 (MQ=255)
                                                                                        cgGACAGGTAATAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAgc                                                              <  2:74533/150‑1 (MQ=255)
                                                                                        cgGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAgc                                                              <  2:3192409/150‑1 (MQ=255)
                                                                                         ggACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAgcg                                                             >  2:1923475/1‑150 (MQ=255)
                                                                                           aCAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCa                                                           >  2:2516697/1‑150 (MQ=255)
                                                                                               gTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCtc                                                       >  2:1975786/1‑149 (MQ=255)
                                                                                               gTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTg                                                       >  2:1456609/1‑150 (MQ=255)
                                                                                               gTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTg                                                       >  2:1256394/1‑150 (MQ=255)
                                                                                                tAAATGTACCGGGTTACACCCGGC‑AAAAAAACGCTACCAACATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCGTCCAGCGCATCTGc                                                      <  2:2641333/150‑1 (MQ=255)
                                                                                                tAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGTTCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGc                                                      <  1:3131102/150‑1 (MQ=255)
                                                                                                tAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTACTCGCCCGCTTCCAGCGCATCTGc                                                      >  2:2607523/1‑150 (MQ=255)
                                                                                                tAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGc                                                      <  2:1175204/150‑1 (MQ=255)
                                                                                                tAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCAACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCAGCTGc                                                      >  2:2043508/1‑150 (MQ=255)
                                                                                                  aaaGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCg                                                    >  1:1117140/1‑150 (MQ=255)
                                                                                                          gggTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAAccagc                                            <  2:1204245/150‑1 (MQ=255)
                                                                                                          gggTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAAccagc                                            <  2:3250119/150‑1 (MQ=255)
                                                                                                              tAAACCCGGCAAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAAccagccag                                         >  1:1371197/1‑150 (MQ=255)
                                                                                                              tAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAgg                                        >  2:1595175/1‑150 (MQ=255)
                                                                                                                 aCCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAac                                     <  1:2130324/150‑1 (MQ=255)
                                                                                                                  cccGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAacc                                    <  1:2648767/150‑1 (MQ=255)
                                                                                                                   ccGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAacca                                   >  1:1265954/1‑150 (MQ=255)
                                                                                                                   ccGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAacca                                   >  2:454911/1‑150 (MQ=255)
                                                                                                                    cGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAaccac                                  >  1:1995379/1‑150 (MQ=255)
                                                                                                                       c‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGa                               <  1:2705053/150‑1 (MQ=255)
                                                                                                                       c‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGa                               <  1:2516566/150‑1 (MQ=255)
                                                                                                                         aaaaaaaCGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTAGAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGAt                              >  1:2486408/1‑150 (MQ=255)
                                                                                                                              aaCGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGcaca                         <  2:2227746/150‑1 (MQ=255)
                                                                                                                                      aaaaaTGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCtt                 >  1:3198768/1‑150 (MQ=255)
                                                                                                                                      aaaaaTGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCtt                 <  2:2791743/150‑1 (MQ=255)
                                                                                                                                       aaaaTGCCCGATCGAGGATCGGGCA‑TTTTGACTTTTACAGCTTATCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCttt                >  1:1331404/1‑150 (MQ=255)
                                                                                                                                          aTGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCa             <  1:102516/150‑1 (MQ=255)
                                                                                                                                           tGCCCGATCGAGGATCGGGCATTTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCa             >  2:2849002/1‑150 (MQ=255)
                                                                                                                                             cccGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAgcg          >  1:1605996/1‑150 (MQ=255)
                                                                                                                                              ccGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAgcgc         <  1:1754256/150‑1 (MQ=255)
                                                                                                                                              ccGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAgcgc         <  2:2434197/150‑1 (MQ=255)
                                                                                                                                              ccGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAgcgc         <  1:1416132/150‑1 (MQ=255)
                                                                                                                                               cGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGcc        >  1:2062030/1‑150 (MQ=255)
                                                                                                                                                  tCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAgg     <  1:1463196/150‑1 (MQ=255)
                                                                                                                                                  tCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAgg     <  1:1674684/150‑1 (MQ=255)
                                                                                                                                                   ccTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGt    >  2:2373848/1‑150 (MQ=255)
                                                                                                                                                   ccTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGt    >  1:781732/1‑150 (MQ=255)
                                                                                                                                                    cTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTa   <  2:1517073/150‑1 (MQ=255)
                                                                                                                                                     tCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAg  <  2:1677820/150‑1 (MQ=255)
                                                                                                                                                     tCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAg  <  2:540783/150‑1 (MQ=255)
                                                                                                                                                     tCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAg  <  2:541103/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGC‑AAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCA‑TTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTAG  >  NC_000913/1931947‑1932244

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: ‑