Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,694,102 C→A C481F (TGC→TTC)  purL ← phosphoribosylformyl‑glycineamide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,694,1020CA98.8% 310.1 / ‑3.3 83C481F (TGC→TTC) purLphosphoribosylformyl‑glycineamide synthetase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (42/40);  total (43/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-01

TCTCGCGCAGTTCAAATTTACCGCCGCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGCCGGGCCGCCGAGAACGACCAGCTTCGCACCG  >  NC_000913/2693956‑2694248
                                                                                                                                                  |                                                                                                                                                  
tctcGCGCAGTTCAAATTTACCGCCGCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCg                                                                                                                                                >  2:220889/1‑151 (MQ=255)
tctcGCGCAGTTCAAATTTACCGCCGCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCg                                                                                                                                                >  1:937922/1‑151 (MQ=255)
tctcGCGCAGTTCAAATTTACCGCCGCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCg                                                                                                                                                >  1:891768/1‑151 (MQ=255)
 ctcGCGCAGTTCAAATTTACCGCCGCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCg                                                                                                                                                >  2:934652/1‑150 (MQ=255)
   cgcgCATTTAAATTTTCCCGCCGCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACg                                                                                                                                             <  2:130222/151‑1 (MQ=255)
    gcgcATTACAAATTTCCCGCCGCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACg                                                                                                                                             <  2:620412/150‑1 (MQ=255)
    gcgcAGTTCAAATTTACCGCCGCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGc                                                                                                                                            >  1:419058/1‑151 (MQ=255)
        aGTCCAAATTTACCGCCGCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGGCCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCt                                                                                                                                           <  1:229691/148‑1 (MQ=255)
          ttaaaaTTTCCCGCCGCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACg                                                                                                                                             <  2:709610/141‑1 (MQ=255)
          ttaaaaTTTCCCGCCGCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTc                                                                                                                                          <  2:342598/144‑1 (MQ=255)
          ttCAATTTTCCCGCCGCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCt                                                                                                                                           <  1:956180/146‑1 (MQ=255)
             aaaTTTACCGCCGCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCAtctc                                                                                                                                    <  1:496175/150‑1 (MQ=255)
              attttCCCGCCGCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCg                                                                                                                                  <  1:155845/149‑1 (MQ=255)
              attttACCGCCGCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCg                                                                                                                                  <  1:943628/149‑1 (MQ=255)
                 ttACCGCCGCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGt                                                                                                                               >  1:1022499/1‑151 (MQ=255)
                   aCCGCCGCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTg                                                                                                                             >  2:734060/1‑151 (MQ=255)
                    ccgccgCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTg                                                                                                                             <  1:244571/150‑1 (MQ=255)
                    ccgccgCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGt                                                                                                                            >  2:964145/1‑151 (MQ=255)
                     cgccggGCCCGCCGTCGCTCACCAGTTCGGGCATGGCGTTAGAAAGACCGCCAGCGACAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCTCTTCCTGGAAGCGACGCTCCATCGCCggg                                                                                                                                <  2:364072/146‑1 (MQ=255)
                        cgcgCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTcgcg                                                                                                                        >  2:685069/1‑151 (MQ=255)
                        cgcgCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTcgcg                                                                                                                        <  2:480117/151‑1 (MQ=255)
                        cgcgCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTcgcg                                                                                                                        >  2:784893/1‑151 (MQ=255)
                         gcgcCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCgcgc                                                                                                                       >  1:303924/1‑151 (MQ=255)
                          cacccgccGTCGCTCACCAGTTCCGGCATGGCGGTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCAtct                                                                                                                                     <  2:972212/134‑1 (MQ=255)
                             ccgccgTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGt                                                                                                                    >  1:1099464/1‑150 (MQ=255)
                                    cGCTCACCAGTTCCGGCATGGCGTTAGAAAGCCCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAg                                                                                                             <  2:798943/150‑1 (MQ=255)
                                    cGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAg                                                                                                             <  2:1248714/150‑1 (MQ=255)
                                      cTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCa                                                                                                           <  1:819548/150‑1 (MQ=255)
                                              ttCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAg                                                                                                   <  1:1145771/150‑1 (MQ=255)
                                                ccGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTc                                                                                                <  1:31829/151‑1 (MQ=255)
                                                   gCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCgg                                                                                              >  1:1237730/1‑150 (MQ=255)
                                                    cATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGc                                                                                             >  2:977958/1‑150 (MQ=255)
                                                    cATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGc                                                                                             >  2:219840/1‑150 (MQ=255)
                                                       ggCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATc                                                                                          >  2:915931/1‑150 (MQ=255)
                                                              gAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGAc                                                                                  >  1:149419/1‑151 (MQ=255)
                                                                aaGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGAcc                                                                                 >  1:1072192/1‑150 (MQ=255)
                                                                aaGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCa                                                                                <  2:518003/151‑1 (MQ=255)
                                                                 aGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGAcc                                                                                 >  2:854144/1‑149 (MQ=255)
                                                                   aCCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGAc                                                                             >  2:646327/1‑151 (MQ=255)
                                                                       ccAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGcc                                                                          >  1:1050952/1‑150 (MQ=255)
                                                                       ccAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCa                                                                         <  2:228973/151‑1 (MQ=255)
                                                                        cAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCATGCTGCCAGCAACGGTCGATCGCTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCAt                                                                        <  2:1006810/151‑1 (MQ=255)
                                                                        cAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCAt                                                                        >  2:338445/1‑151 (MQ=255)
                                                                            gcCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCCGACGCCATa                                                                       >  2:1154974/1‑148 (MQ=255)
                                                                            gcCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATagaa                                                                    <  1:401822/151‑1 (MQ=255)
                                                                               aaCGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGAc                                                                 >  2:903456/1‑151 (MQ=255)
                                                                               aaCGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGAc                                                                 <  1:499272/151‑1 (MQ=255)
                                                                               aaCGTCGTGGATAAACAGGATTGAGTTGGCATCACCAAGCTGCCCGCATCGGTCGATCACTTCCTGGAAGCGACGCTCCAGCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGGCCAGACGCCATagaaga                                                                  <  2:31899/150‑1 (MQ=255)
                                                                                  gtcgtGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCt                                                              >  1:483905/1‑151 (MQ=255)
                                                                                  gtcgtGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCt                                                              <  2:584248/151‑1 (MQ=255)
                                                                                    cgtGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGc                                                            >  1:721382/1‑151 (MQ=255)
                                                                                      tGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGGTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCGGcac                                                          >  2:1175546/1‑151 (MQ=255)
                                                                                        gATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGcacc                                                         >  1:30883/1‑150 (MQ=255)
                                                                                         aTAAACATGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGcaccac                                                       <  2:929096/151‑1 (MQ=255)
                                                                                         aTAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGcaccac                                                       >  2:441477/1‑151 (MQ=255)
                                                                                         aTAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGcaccac                                                       >  2:1109648/1‑151 (MQ=255)
                                                                                            aaCAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGc                                                    <  2:346792/151‑1 (MQ=255)
                                                                                             aCAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGcc                                                   >  2:303955/1‑151 (MQ=255)
                                                                                                 gATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCa                                                  >  2:973181/1‑148 (MQ=255)
                                                                                                      ggTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGAt                                          >  2:975837/1‑151 (MQ=255)
                                                                                                         tGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCACCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGt                                        <  2:1108114/150‑1 (MQ=255)
                                                                                                            cATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCAt                                    >  2:206846/1‑151 (MQ=255)
                                                                                                             aTCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCAtt                                   >  1:748878/1‑151 (MQ=255)
                                                                                                                 ccAAGCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGCCg                               <  1:605035/151‑1 (MQ=255)
                                                                                                                   aaGCTGCCAGCAACGTTCGATCACTTCCTGCAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAGGTCGAGTTCGGCATCAGACTGACCAGACGCCATAGAAGCCGCTGCACCACCGCCAAGACCGATGTTCATTGCCggg                             <  1:1130302/151‑1 (MQ=255)
                                                                                                                     gCTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGCCGGgcc                           <  1:881909/151‑1 (MQ=255)
                                                                                                                     gCTGCCAGCAACGGTCGATCACTCCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGCCGGgcc                           <  1:306404/151‑1 (MQ=255)
                                                                                                                      cTGCCAGCAACGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGTTGCACCACCGCCAAGACCGATGTTCATTGCCGGgcc                           <  2:675020/150‑1 (MQ=255)
                                                                                                                              aaCGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGCCGGGCCGCCGAGAAc                  >  1:454803/1‑151 (MQ=255)
                                                                                                                               aCGGTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGCCGGGCCGCCGAGAACg                 >  2:838863/1‑151 (MQ=255)
                                                                                                                                  gTCGATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGCCGGGCCGCCGAGAACGAcc              >  1:987763/1‑151 (MQ=255)
                                                                                                                                     gATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGCCGGGCCGCCGAGAACGACCa             <  2:1152941/149‑1 (MQ=255)
                                                                                                                                     gATCACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGCCGGGCCGCCGAGAACGACCAGc           >  1:386038/1‑151 (MQ=255)
                                                                                                                                     gATCACTTCCTGGAAGCGACGCTCCATCTCCGGGGTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGCCGGGCCGCCGAGAACGACCAGc           >  2:248240/1‑151 (MQ=255)
                                                                                                                                        cACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGCCGCTGCACCACCGCCAAGACCGATGTTCATTGCCGGGCCGCCGAGAACGACCAGCtt         <  2:798571/150‑1 (MQ=255)
                                                                                                                                        cACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGCCGGGCCGCCGAGAACGACCAGCtt         <  2:1188297/150‑1 (MQ=255)
                                                                                                                                        cACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGCCGGGCCGCCGAGAACGACCAGCTTc        <  2:521854/151‑1 (MQ=255)
                                                                                                                                        cACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGCCGGGCCGCCGAGAACGACCAGCTTc        <  2:1099411/151‑1 (MQ=255)
                                                                                                                                        cACTTCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGCCGGGCCGCCGAGAACGACCAGCTTc        <  1:1109701/151‑1 (MQ=255)
                                                                                                                                           ttCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCCAGACCGAAGTTCATTGCCGGGCCGCCGAGAACGACCAGCTTCGc      >  1:627484/1‑150 (MQ=255)
                                                                                                                                           ttCCTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGCCGGGCCGCCGAGAACGACCAGCTTCGc      >  2:287418/1‑150 (MQ=255)
                                                                                                                                             ccTGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGCCGGGCCGCCGAGAACGACCAGCTTCGCa     <  2:478639/149‑1 (MQ=255)
                                                                                                                                               tGGAAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGCCGGGCCGCCGAGAACGACCAGCTTCGCACCg  >  2:712724/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TCTCGCGCAGTTCAAATTTACCGCCGCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGCCGGGCCGCCGAGAACGACCAGCTTCGCACCG  >  NC_000913/2693956‑2694248

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