Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,815,908 C→A 5.8% A178S (GCG→TCG)  gshA ← glutamate‑‑cysteine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,815,9080CA 5.8% 265.2 / 2.9 79A178S (GCG→TCG) gshAglutamate‑‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base C (44/30);  new base A (3/2);  total (47/32)
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

GTAATACATACCGCACTCGGTTTTCTCAAACGGCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAGCGATTTTTCAGCCCTTCACGATACAGCGTTT  >  NC_000913/2815759‑2816057
                                                                                                                                                     |                                                                                                                                                     
gTAATACATACCGCACTCGGTTTTCTCAAACGGCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGc                                                                                                                                                       <  2:233048/150‑1 (MQ=255)
 tAATACATACCGCACTCGGTTTTCTCAAACGGCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCa                                                                                                                                                      >  2:1339979/1‑150 (MQ=255)
    tacatacCGCACTCGGTTTTCTCAAACGGCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGaa                                                                                                                                                   <  1:539495/150‑1 (MQ=255)
      catacCGCACTCGGTTTTCTCAAACGGCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAAt                                                                                                                                                 <  2:1237678/150‑1 (MQ=255)
      catacCGCACTCGGTTTTCTCAAACGGCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAAt                                                                                                                                                 >  2:1198273/1‑150 (MQ=255)
       atacCGCACTCGGTTTTCTCAAACGGCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAAGGATAGTAATTGCGGATAACGCGGAATTAGCCCGCAGAAAtt                                                                                                                                                >  1:1114974/1‑150 (MQ=255)
          ccGCACTCGGTTTTCTCAAACGGCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTtc                                                                                                                                             >  1:813336/1‑150 (MQ=255)
          ccGCACTCGGTTTTCTCAAACGGCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCAAACCGAAACGATAGTAAATGCGGATAACGCGGAAAGAGCCCGAAGAAAtattg                                                                                                                                             >  1:1605589/1‑149 (MQ=255)
             cACTCGGTTTTCTCAAACGGCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTctct                                                                                                                                          >  2:1429325/1‑150 (MQ=255)
                 cGGTTTTCTCAAACGGCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGAAGAAATTTTCTCTTTgg                                                                                                                                      >  1:1404012/1‑150 (MQ=255)
                 cGGTTTTCTCAAACGGCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGAAGAAATTTTCTCTTTgg                                                                                                                                      >  1:1409804/1‑150 (MQ=255)
                  ggTTTTCTCAAACGGCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGc                                                                                                                                     >  2:1460940/1‑150 (MQ=255)
                      ttCTCAAACGGCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCa                                                                                                                                 <  2:1330109/150‑1 (MQ=255)
                      ttCTCAAACGGCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCa                                                                                                                                 <  1:1533176/150‑1 (MQ=255)
                        ctcAAACGGCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGCACCGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAgc                                                                                                                               >  2:915547/1‑150 (MQ=255)
                        ctcAAACGGCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAgc                                                                                                                               >  2:878809/1‑150 (MQ=255)
                          cAAACGGCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAgcgc                                                                                                                             >  1:1514332/1‑150 (MQ=255)
                          cAAACGGCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAgcgc                                                                                                                             >  2:849358/1‑150 (MQ=255)
                          cAAACGGCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAgcgc                                                                                                                             >  2:870642/1‑150 (MQ=255)
                             aCGGCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCg                                                                                                                          <  1:1633906/150‑1 (MQ=255)
                               ggCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAATCAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCgag                                                                                                                        <  2:813336/150‑1 (MQ=255)
                                                      cAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCa                                                                                                 <  1:759700/150‑1 (MQ=255)
                                                      cAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGAGCCCGAGATATCACCGCACTTCGCTTGCCa                                                                                                 <  1:878809/150‑1 (MQ=255)
                                                       aGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAg                                                                                                >  1:334754/1‑150 (MQ=255)
                                                        ggAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGa                                                                                               >  2:128092/1‑150 (MQ=255)
                                                          aaaGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACACAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAAt                                                                                             >  1:995741/1‑150 (MQ=255)
                                                             gaagaaCAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCgag                                                                                                                        <  2:618135/120‑1 (MQ=255)
                                                             gaagaaCAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCgag                                                                                                                        >  1:618135/1‑120 (MQ=255)
                                                             gaagaaCAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGcc                                                                                          >  2:283944/1‑150 (MQ=255)
                                                             gaagaaCAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGcc                                                                                          >  1:1650416/1‑150 (MQ=255)
                                                             gaagaaCAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACAGCACTTCGCTTGCCAGAATGcc                                                                                          >  1:1192165/1‑150 (MQ=255)
                                                               agaaCAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCAt                                                                                        >  1:878132/1‑150 (MQ=255)
                                                                gaaCAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCAtt                                                                                       >  1:355141/1‑150 (MQ=255)
                                                                 aaCAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTg                                                                                      <  1:1460940/150‑1 (MQ=255)
                                                                  aCAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAATTAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGATATATCACCGCACTTCGCTTGCCAGAATGCCATTgg                                                                                     <  2:1230179/150‑1 (MQ=255)
                                                                  aCAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTgg                                                                                     <  1:1184792/150‑1 (MQ=255)
                                                                   cAAATCGCCGGAGATTCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGc                                                                                    <  1:1387725/150‑1 (MQ=255)
                                                                   cAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGc                                                                                    >  1:8223/1‑150 (MQ=255)
                                                                     aaTCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCaa                                                                                  <  1:515188/150‑1 (MQ=255)
                                                                             gagaTGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAAt                                                                          <  1:253492/150‑1 (MQ=255)
                                                                               gaTGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTg                                                                        >  1:1150676/1‑150 (MQ=255)
                                                                                aTGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGt                                                                       <  1:1183275/150‑1 (MQ=255)
                                                                                aTGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGt                                                                       <  2:196499/150‑1 (MQ=255)
                                                                                  gCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAg                                                                     <  2:1530018/150‑1 (MQ=255)
                                                                                   cACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGt                                                                    >  2:1582168/1‑150 (MQ=255)
                                                                                   cACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGt                                                                    >  2:1320194/1‑150 (MQ=255)
                                                                                     ccAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGc                                                                  >  1:495438/1‑150 (MQ=255)
                                                                                        aaCAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACg                                                               >  1:1243605/1‑150 (MQ=255)
                                                                                               aaGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGaaa                                                        >  1:89380/1‑150 (MQ=255)
                                                                                                 ggAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATg                                                      >  1:1374922/1‑150 (MQ=255)
                                                                                                  gAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATgg                                                     >  2:349721/1‑150 (MQ=255)
                                                                                                   aaTGACCCAACCGAAACGATAGTAATTGCGGATAACTCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAg                                                                     <  2:918803/133‑1 (MQ=255)
                                                                                                   aaTGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAg                                                                     >  1:918803/1‑133 (MQ=255)
                                                                                                       aCCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGATATGGTTTGc                                                >  2:910515/1‑150 (MQ=255)
                                                                                                        cccAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCCCGCCGGAAATGGTTTGCa                                               >  2:941862/1‑150 (MQ=255)
                                                                                                          aaaCCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATc                                             <  2:995741/149‑1 (MQ=255)
                                                                                                               gAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAgcgc                                        >  1:204179/1‑150 (MQ=255)
                                                                                                                 aaCGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAgcgcgc                                      <  1:870642/150‑1 (MQ=255)
                                                                                                                  aCGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCCCGCCGGAAATGGTTTGCATCAGCGCGcc                                     >  1:1503563/1‑150 (MQ=255)
                                                                                                                     aTAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTa                                  <  1:1614688/150‑1 (MQ=255)
                                                                                                                       aGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAGc                                <  1:803517/150‑1 (MQ=255)
                                                                                                                        gTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAGCg                               >  2:64310/1‑150 (MQ=255)
                                                                                                                                ggATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAGCGATTTTTCa                       >  1:1444223/1‑150 (MQ=255)
                                                                                                                                     aCGCGGAAAAAGCCCGCAGAAATTTTCTCTTTGGCATAAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATGGTAGAGCACGCCGGAAATGGATTGCATCAGCGAGCCGGAGCGATATTTAAGaccg                  >  2:764302/1‑149 (MQ=255)
                                                                                                                                      cgcgGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAGCGATTTTTCAGCCCtt                 >  2:1319227/1‑150 (MQ=255)
                                                                                                                                      cgcgGAAATAGCCAGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCTATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAGCGATTTTTCAGCCCtt                 <  2:752211/150‑1 (MQ=255)
                                                                                                                                       gcgGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAGCGATTTTTCAGCCCTTc                >  1:99007/1‑150 (MQ=255)
                                                                                                                                       gcgGAAATAGCCCGAAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAGCGATTTTTCAGCCCTTc                <  2:1404012/150‑1 (MQ=255)
                                                                                                                                       gcgGAAATAGCCCGAAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAGCGATTTTTCAGCCCTTc                <  2:1409804/150‑1 (MQ=255)
                                                                                                                                         ggAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAGCGATTTTTCAGCCCTTCAc              <  1:283944/150‑1 (MQ=255)
                                                                                                                                         ggAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAGCGATTTTTCAGCCCTTCAc              <  2:1374922/150‑1 (MQ=255)
                                                                                                                                          gAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTCCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAGCGATTTTTCAGCCCTTCACg             >  1:553475/1‑150 (MQ=255)
                                                                                                                                             aTAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCTAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAGCGATTTTTCAGCCCTTCACGATa          <  2:699323/150‑1 (MQ=255)
                                                                                                                                              tAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAGCGATTTTTCAGCCCTTCACGATAc         <  1:849358/150‑1 (MQ=255)
                                                                                                                                                gCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAGCGATTTTTCAGCCCTTCACGATACAg       <  2:449774/150‑1 (MQ=255)
                                                                                                                                                gCCCGCAGAAATTTTATCTTTGGCATCAGAGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAGCGATTTTTCAGCCCTTCACGATAccg       >  2:1192527/1‑148 (MQ=255)
                                                                                                                                                 cccGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAGCGATTTTTCAGCCCTTCACGATACAGc      >  2:730575/1‑150 (MQ=255)
                                                                                                                                                  ccGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTg                                                 >  2:233421/1‑106 (MQ=255)
                                                                                                                                                  ccGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTg                                                 <  1:233421/106‑1 (MQ=255)
                                                                                                                                                  ccGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAGCGATATTTCAGCCCTTCACGATACAGCg     >  2:1292706/1‑150 (MQ=255)
                                                                                                                                                     cAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCTCTCTTCGCTTTCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAGCGTTTTTTCAGCCCTTCACGATACAGCGttt  <  2:825097/150‑1 (MQ=255)
                                                                                                                                                     cAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAGCGATTTTTCAGCCCTTCACGATACAGCGttt  <  1:1582168/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTAATACATACCGCACTCGGTTTTCTCAAACGGCAGCGACGTTGGTTTTCCTTGCAGGAAAGAAGAACAAATCGCCGGAGATGCACCAAACAGATAAGGAATGACCCAACCGAAACGATAGTAATTGCGGATAACGCGGAAATAGCCCGCAGAAATTTTCTCTTTGGCATCAGCGCCCGAGATATCACCGCACTTCGCTTGCCAGAATGCCATTGGCAAAGAGAAATTGTAGTGCACGCCGGAAATGGTTTGCATCAGCGCGCCGTAGCGATTTTTCAGCCCTTCACGATACAGCGTTT  >  NC_000913/2815759‑2816057

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