Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,286,033 C→T A447T (GCG→ACG)  acs ← acetyl‑CoA synthetase (AMP‑forming)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,286,0330CT95.5% 182.3 / NA 68A447T (GCG→ACG) acsacetyl‑CoA synthetase (AMP‑forming)
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (33/31);  minor base G (1/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 6.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCAT  >  NC_000913/4285787‑4286279
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
tCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTAGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGt                                                                                                                                                                                                                                                    >  2:113845/1‑251 (MQ=255)
tCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGAGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGt                                                                                                                                                                                                                                                    >  1:242063/1‑251 (MQ=255)
    cAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGCAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGc                                                                                                                                                                                                                                                 >  1:172672/1‑250 (MQ=255)
     aGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTCCCTTCGTTATCGACCAGCGCCGGTTGccc                                                                                                                                                                                                                                               >  2:326842/1‑249 (MQ=255)
       aCGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGc                                                                                                                                                                                                                                             >  1:52622/1‑251 (MQ=255)
             aCCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCgaaga                                                                                                                                                                                                                                       <  1:159528/251‑1 (MQ=255)
                ggAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACg                                                                                                                                                                                                                                    <  1:48344/251‑1 (MQ=255)
                     cGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACgg                                                                                                                                                                                                                               >  2:265565/1‑250 (MQ=255)
                      gTTCAGCACGGCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCAGCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACgtg                                                                                                                                                                                                                              <  1:93607/251‑1 (MQ=255)
                        tCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACgtgt                                                                                                                                                                                                                             <  1:33670/250‑1 (MQ=255)
                         cAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGtt                                                                                                                                                                                                                            <  1:102215/250‑1 (MQ=255)
                         cAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTg                                                                                                                                                                                                                           >  1:393117/1‑251 (MQ=255)
                                  gtcCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCGACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCgg                                                                                                                                                                                                                  <  1:171226/251‑1 (MQ=255)
                                    ccACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCt                                                                                                                                                                                                                >  1:27067/1‑251 (MQ=255)
                                             cGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCg                                                                                                                                                                                                       <  2:331410/251‑1 (MQ=255)
                                                  aTCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTGCGGTGTCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAg                                                                                                                                                                                                   <  1:176942/250‑1 (MQ=255)
                                                    ccAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAgcgc                                                                                                                                                                                                <  1:42662/251‑1 (MQ=255)
                                                      aGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGcc                                                                                                                                                                                               >  1:104412/1‑250 (MQ=255)
                                                                    tCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGt                                                                                                                                                                                   <  1:290222/248‑1 (MQ=255)
                                                                     aaTCGCGAAGAGCACCGTCCCCACTGACATACATATTTATGACGGTGCAGAAGTAGGTCTGTTCAAAACGTTCGTCATCGCCAAACAGCGTACGCGCCTGACCAGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGGCGGTTGCACGCCGCAGAACGGAAGTGTTGCAGCACCGGCTTTCAGCGCGGTATCGCCCGGCCGCGGGGCGATc                                                                                                                                                                               >  2:40319/2‑251 (MQ=255)
                                                                        cgcgACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATg                                                                                                                                                                            >  1:139629/1‑251 (MQ=255)
                                                                           gACGCGCGCCGTCGCCGCTGAAATACATAGTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGAGCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTGCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCATCTTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCAAGaaa                                                                                                                                                                         <  1:373250/251‑1 (MQ=255)
                                                                                 cgcCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTCCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGGCGGTGATTCCCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAGccgccg                                                                                                                                                                   <  1:70668/251‑1 (MQ=255)
                                                                                   ccGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGt                                                                                                                                                                 >  1:67642/1‑251 (MQ=255)
                                                                                         ccgcTGAAATACATATTTTGGAAGGTGGAGAAGTGGGTCTTTTCAAAACGTTCGTGATCGCCAAACAGCGTACCCGCCTGCCCCGGCCAGGGGTCGGTGTTTCCCAGGCTCCCTTCGGGGGCCCCCTCCAGCGGGTTACCTTCGTTATCGCCCAGCGTCGGTTGCACGCCGAAGAACGGACGGTTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGt                                                                                                                                                           <  2:52622/251‑1 (MQ=255)
                                                                                          cgcTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGCTTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTc                                                                                                                                                          >  1:343499/1‑251 (MQ=255)
                                                                                              gAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAACCAGCGTACGCGCATGACCCGGCCCGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCCGCGGGTTACCTTCGTTATCGACCAGCGTAGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCCGCTCGGTAGGGACAGACAGCGGGGTGATCAAAAAACAGCCGGTTTCGTTCTGcc                                                                                                                                                      >  2:318286/1‑251 (MQ=255)
                                                                                                 atatatatTTTTGAAGGTGGAGAAGTATGTATGCTCAAACCTTTCGTGATCGCCAAACCGCGTTCGCCCCTGCCCCGGCCAGGAGTCGGTGATTACCAGGCTCCCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCCCCAGGCAGCGGGGGGTTCATGAAACCGCCGGTTTCGGTCTGccacc                                                                                                                                                   <  1:220648/247‑1 (MQ=255)
                                                                                                  tACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGccacc                                                                                                                                                   >  2:68532/1‑250 (MQ=255)
                                                                                                   aCATATTTTTGAAGGTGGAGAAGTAGGTCTGTCCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGGGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAg                                                                                                                                                 <  1:363648/251‑1 (MQ=255)
                                                                                                     atatTTTTGCAGGAGGAGAAGCAGGACTGTTAAAAACATTCGTGATCGCCAAACAGCGTACGCGCCTGACCAAGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGACCCATCCAGCGGGTTACCTTCGTTATAGCCCAGCGTCGGTTGCACGCAGACGAACGGAAGTGTTGCCGAACCGACTTTCAGCTCGGTCGCGCCAGGCAGCGGGGGGATCATGAAACCGCCGGTTTCGGTCTGCCACCAggg                                                                                                                                               >  2:402408/1‑250 (MQ=255)
                                                                                                         ttttGAAGGTGGAGAAGAAGGTCTGTTAAAAGCGTCCGGGATCGCCAAACCGCGTACCCGCCTGACCCGCCCATGAGTCGGTGATTCCCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTCCCTTCGTTATCGCCCAGCGTCGGTGGCACGCCGAAGACCGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGG‑‑TGCCACCAGGTATCGa                                                                                                                                          <  2:233568/250‑1 (MQ=255)
                                                                                                            tGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCAGGAAACCGCCGGTTTCGGTCTGCCACCACGTATCGacc                                                                                                                                        >  2:313220/1‑251 (MQ=255)
                                                                                                              aaGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGaccac                                                                                                                                      >  2:288497/1‑251 (MQ=255)
                                                                                                                 gTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCgg                                                                                                                                   >  1:266383/1‑251 (MQ=255)
                                                                                                                    gagaAGTAGGTCTGTTCAAACCGTCCGTGATCCCCAAACAGCTTACCCGCCGGCCCCGGCCAGGATTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTCCCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAAACCGGACGTGTTGCCGACCCGCCTTTCAGCTCGTTAGCCCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACa                                                                                                                                <  2:104412/251‑1 (MQ=255)
                                                                                                                            ggTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCGGACCCGGCCAGGAGTCGGTGATTCCCAGGCTACCTTCGGTGGCCCCCCCCAGCGGGTTACCTTCGTTAACGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGt                                                                                                                        <  1:247299/251‑1 (MQ=255)
                                                                                                                                 gTTCACAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCg                                                                                                                   >  2:204272/1‑251 (MQ=255)
                                                                                                                                  ttCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCg                                                                                                                   >  1:360055/1‑250 (MQ=255)
                                                                                                                                      aaaCGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCCCGCCGAAGAACGGACGGGTTGCCGCACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCTCCAGGTATCGACCACCGGACATTTCTCGTTGCCGAtctt                                                                                                              >  2:65602/1‑248 (MQ=255)
                                                                                                                                         cgttcggtaTCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGTTTACCAGTCTCCCTTCGGTGGCCCCCCCCAGCGGGTTACCTTCGTTATCGCCCAGCGTCGGTTGCACGCCGAAGAACGGACGGGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGGGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGAttttttt                                                                                                           <  1:391041/243‑1 (MQ=255)
                                                                                                                                                gATCGCCAAACAGCGTACGCGCCTGCCCCGGCCAGGGGGCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTCCCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGACCGGATGTGTTGCCGAACCGGCTTTCAGCTCGGGAGCGCCGGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGGATCGACCACCGGACATTTCTCGTTGCCGATTTTTTGCCAGTAc                                                                                                    <  2:263285/251‑1 (MQ=255)
                                                                                                                                                gATCGCCAAACAGCGTACGCGACTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGCTCATGACACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGCCATTTCTCGTTGCCGATTTTTTTCCAGTAc                                                                                                    >  2:184578/1‑251 (MQ=255)
                                                                                                                                                  tCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCAAGTAccc                                                                                                  >  2:297450/1‑250 (MQ=255)
                                                                                                                                                   cGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGGGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGGATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTCCCAc                                                                                                 >  2:192689/1‑251 (MQ=255)
                                                                                                                                                        aaCAGCGTCCGCGCCGGACCCGGCCAGGAGTCGGTGATTACCAGGCTCCCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCa                                                                                            <  1:159461/251‑1 (MQ=255)
                                                                                                                                                                   cgcCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCCGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCGGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCCGTACCACTCCCACCCTTCCGGGt                                                                                 >  1:330002/1‑251 (MQ=255)
                                                                                                                                                                    gcCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCACCTCCAGCGGGTTACCTTCGTTATCGACCAGAGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTCTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGCTTTTTTTCCAGTACCACTCCCACGCTTCCGGGtt                                                                                >  1:150875/1‑251 (MQ=255)
                                                                                                                                                                     ccTGACCCGCCCAGGTGTCGGTGTTTACCAGGCTGCCTTCTGTGGCCCCCTCCAGCGGGTTACCTCCTTTTGCGCCCTGCGTCGTTTCCACGCCGAATAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGGGATCAGGAAACCGCCGGTTTCGGTCTGCCACCAGGGATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTa                                                                               <  1:217987/251‑1 (MQ=255)
                                                                                                                                                                        gACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCCCTCCCCCGCTTCCGGGTTAAtt                                                                            >  1:8721/1‑251 (MQ=255)
                                                                                                                                                                                    gggtggtgtttacaaggtccccTTCGGTGGCCCCCACCAGCGGTTTTCCTTCGTTATCGCCCACCTTCGTTGGCCCGCCGAAGAACGGACGGTTTGCCGACCCGCTTTTCAGCTCGGTCCCGCTGGGCAGCGGGGTGATCATTAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTTCCCCTCCCACGCTTCCGGGTTAATTGGCTCGTCCACg                                                                <  2:287211/231‑1 (MQ=255)
                                                                                                                                                                                     gTCGGTGATTACCAGGCTCCCTCCGGTGGCCCCCTCCAGCGGGTTACCTCCGTTATCGCCCAGCGTCGGTTGCACGCGGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACgg                                                               <  1:73789/251‑1 (MQ=255)
                                                                                                                                                                                          tGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCCGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCg                                                          >  2:45402/1‑251 (MQ=255)
                                                                                                                                                                                              tACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCgaga                                                       <  1:288497/250‑1 (MQ=255)
                                                                                                                                                                                                      tACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGc                                              >  2:117785/1‑251 (MQ=255)
                                                                                                                                                                                                       aCCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGCACCGCGAATGCGCAGcc                                             >  2:102175/1‑250 (MQ=255)
                                                                                                                                                                                                           tCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTAGCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGTCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGcgac                                           <  1:129061/249‑1 (MQ=255)
                                                                                                                                                                                                                 gCCCCCTCCAGCGGGTTACCTTTGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGCCGTGTTGCCGACCCGGCTGTCAGCTCGGTAGCGCCAGCCAGCGGGGTTATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTTCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAAcggtc                                   <  1:240530/251‑1 (MQ=255)
                                                                                                                                                                                                                 ccccccTCCAGCGGGTTACCTTCGTTATCGACCAGCGGCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAAcggtc                                   <  1:281856/250‑1 (MQ=255)
                                                                                                                                                                                                                          aGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCGGGCAGCGGGGTGATCAGGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTc                          <  1:131275/251‑1 (MQ=255)
                                                                                                                                                                                                                          aGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTc                          <  1:68532/251‑1 (MQ=255)
                                                                                                                                                                                                                             gggTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGAt                       >  2:58050/1‑251 (MQ=255)
                                                                                                                                                                                                                               gTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTAAATTGGCTCg                                                                      >  1:248254/1‑202 (MQ=255)
                                                                                                                                                                                                                                     ttCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGCACGGTCGGTGCCTTCGCTCGCTTTAt               >  2:154641/1‑251 (MQ=255)
                                                                                                                                                                                                                                      tcgttcgcgcctaGCGGCGGTCGCACGCCGAAGACCGGACGTGTTGCCGACCCGTCTTTCAGCTCGGTAGCGCCAGGCAGCGGGTTGATCGTGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCTTTa                <  1:65602/237‑1 (MQ=255)
                                                                                                                                                                                                                                                ccAGCGTCGGTTGCGCGCGGAAGAACGGAGGTGTTGCCGAACCGGCTTTCAGCTCGGGAGCGCCAGGCAGGGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTTCCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGcc    <  2:278061/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  aGCGTCGGTTGCACGCCGAAGAACGGAGGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCAt  <  1:117785/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  aGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCAt  >  1:252229/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
TCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCAT  >  NC_000913/4285787‑4286279

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