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.7% 196.0 / NA 72A447T (GCG→ACG) acsacetyl‑CoA synthetase (AMP‑forming)
Reads supporting (aligned to +/- strand):  ref base C (0/2);  major base T (46/20);  minor base G (2/1);  total (49/23)
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

CAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCG  >  NC_000913/4285791‑4286275
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
cAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGACAAACAGCGTGCGCGCCTGACCCGGCCAGGAGTCGGGGATTACCCGGCCACCTTCGGTGGCCCCACCCCGCGGGTTACCTTCGTTATCCAACAGCCTCGGTTGcc                                                                                                                                                                                                                                             >  2:461601/1‑250 (MQ=255)
 ggACGGTGACCGGCGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCAc                                                                                                                                                                                                                                            <  2:5816/250‑1 (MQ=255)
 aGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCTAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCAc                                                                                                                                                                                                                                            >  1:66961/1‑251 (MQ=255)
  gACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACg                                                                                                                                                                                                                                           >  1:442024/1‑251 (MQ=255)
      gTGAACGGTGACGTTCACCACGCCGGACACACTCCCGTTGATCCTGTAATAGCCACTTTAATTGCGACGCGCGCCGTCGCTGCGTAAATACAGATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAATCGTTCGTGATCGCCAAACAGCGGAGGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTGACTTCGGGGGCCCCCTCCAGCGGGTTACCCTCGTTATCGACCAGCGTAGGTTGCACGCCga                                                                                                                                                                                                                                       <  2:312602/251‑1 (MQ=255)
        gACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCga                                                                                                                                                                                                                                       >  2:345549/1‑249 (MQ=255)
              ggaCGTTCAGCCCGTCGTCCACACGCCCGGTTATCCAGTAATAGCCAGCTTCATCGCGGCGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGa                                                                                                                                                                                                                               <  2:542266/250‑1 (MQ=255)
                         aCGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGCAGAACGGACGTGTTTCCg                                                                                                                                                                                                                     >  2:14434/1‑250 (MQ=255)
                           gtcgtcCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGGGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGGGTTGcccacc                                                                                                                                                                                                                  >  2:46431/1‑247 (MQ=255)
                             cgtcCACCCGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCg                                                                                                                                                                                                                <  1:121618/251‑1 (MQ=255)
                             cgtcCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCg                                                                                                                                                                                                                <  1:345549/251‑1 (MQ=255)
                                    acGCCCGGTTATCAAGAACTAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGACATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGCGGCGATGATTACCAGGCTACCTTCGACGGCCCCCTCCCGCGGGTTACCTTCGTTATCGACCAGCGGCGGTTGCACGCCGAAGACCGGaacagtatccgaaccggatttac                                                                                                                                                                                                         >  2:268488/1‑229 (MQ=255)
                                       cccGGTTAGCCAGTAATAGCCAGCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGGTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCt                                                                                                                                                                                                      <  1:365760/251‑1 (MQ=255)
                                         cGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGGCGGTTGCACGCCGGAGAACGGACGTGTTGCAGAACCGGCGTTCAGCTCg                                                                                                                                                                                                    >  2:398685/1‑251 (MQ=255)
                                              aTCCAGTAATAGCCTACTCCAGCGCGACGCGCTCCGTCGCAGCTGAAATACAAATTTTTGAAGGTGGAGAAGTAGGTCGGGTCAAAACGGTCGTGATCGCCAAACAGCGGACGCGCCGGACCCGGCCAGGAGTCGGTGAGTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCTGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAg                                                                                                                                                                                                <  2:204355/250‑1 (MQ=255)
                                               tCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTATGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGAGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTCACGGCCAGAGTCGGTTGCAGGCCGAAGAACGGACGTGTTGACGAACAGGCTTTCCGCTCGGTaacg                                                                                                                                                                                              >  2:432910/1‑248 (MQ=255)
                                                        aGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGGGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGACAACCGGCCGTGGTCACGAACCGGCTTTCAGCtacgtaacaccagtcagc                                                                                                                                                                                     >  2:362196/1‑234 (MQ=255)
                                                               ttCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTg                                                                                                                                                                               >  1:224094/1‑250 (MQ=255)
                                                                       gACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGCCCAGCGTCGGTTGCACGCCGATGAACGGACGTGTTGCCGAACCGTCTTTCAGCTCGGTTGCGCCAGGTAGCGGGGTGGTCATGaa                                                                                                                                                                      >  2:148158/1‑250 (MQ=255)
                                                                             cgcCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGCACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAAccgccg                                                                                                                                                               >  1:38206/1‑251 (MQ=255)
                                                                             cgcCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGCACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGCGATCCAGAAAccgccg                                                                                                                                                               >  2:424837/1‑251 (MQ=255)
                                                                                    gccgcTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTCCCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACTGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAACGCCAGGCCGCGGGGTGATCACGAAACCGCCGGTTTCgg                                                                                                                                                        >  2:321468/1‑251 (MQ=255)
                                                                                        cTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTg                                                                                                                                                    >  1:374777/1‑251 (MQ=255)
                                                                                          gAAATACATATTTTGGAAGGTGTAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGCCCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTTTTGCCGAACCGGCTTTCAGCTCGTTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTg                                                                                                                                                    <  1:136925/249‑1 (MQ=255)
                                                                                            aaTACATATTTTTGACGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGACCCGGCTTTCCGCTCGGTAGCGCCAGGCCGCGGGGTGATCATGCAACCGCCGGTTTCGGTCTGccac                                                                                                                                                >  2:213227/1‑251 (MQ=255)
                                                                                            aaTACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGCCGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGccac                                                                                                                                                >  1:385072/1‑251 (MQ=255)
                                                                                                 atatTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGACTGCCACCAggg                                                                                                                                           >  1:12788/1‑250 (MQ=255)
                                                                                                    tttttGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCa                                                                                                                                                                                                         >  1:278056/1‑187 (MQ=255)
                                                                                                    tttttGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTCCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCa                                                                                                                                                                                                         <  2:278056/187‑1 (MQ=255)
                                                                                                       ttGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGAt                                                                                                                                                                             >  1:420134/1‑212 (MQ=255)
                                                                                                       ttGAAGGTGGAGAAGTAGGGCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCATCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGAt                                                                                                                                                                             <  2:420134/212‑1 (MQ=255)
                                                                                                          aaGGTGGAGAAGTAGGTCTGTTCAAAACGTTCTTGATCTCCAAACAGCGTCCGCGCCTGACCCGGCCAGGAGTCGGTGATTCCCAGGTTCCCTCCGGTGGCCCCCTCCAGCGGGTTACCTGCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGGTTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGaccac                                                                                                                                  <  1:369780/251‑1 (MQ=255)
                                                                                                           aGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTAGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATAATGAAACCGCCGGTTTCGGTCCGCCACCAGGTATCGaccacc                                                                                                                                 >  1:254379/1‑251 (MQ=255)
                                                                                                            ggTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCg                                                                                                                                >  1:387495/1‑251 (MQ=255)
                                                                                                                 agaAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACAt                                                                                                                           >  1:473461/1‑251 (MQ=255)
                                                                                                                  gaAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTCCGCGCTTGCCCCGGCCAGGAGTCGGTGATTACCAGGCTCCCTTCGGTGGCCCCCTCCAGCGGTTTACCGTCGTTATCGACCAGCGTCGTTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTAGCGACCACCGGACAtt                                                                                                                          <  1:319757/251‑1 (MQ=255)
                                                                                                                   aaGTAGTTCTGTCCAAAACGTTCGAGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTTATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGGCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTTATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACAttt                                                                                                                         <  1:23202/251‑1 (MQ=255)
                                                                                                                   aaGTAGGTCTGTTCAAACCGTTCGTGATCGCCAAACAGCGTACGCGCCGGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGGGGCCCCCTCCAGCGGGTTCCCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACAttt                                                                                                                         <  1:519154/251‑1 (MQ=255)
                                                                                                                     gTAGGTCTGTTAAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACACGGCCAGGAGTCGGGGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCATCGGGGGGGTCATGAGACCGCCGGTTTCGGTCTGTCACCAGGttagggacaacggacatttat                                                                                                                       >  2:153329/1‑231 (MQ=255)
                                                                                                                         gTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGtt                                                                                                                   >  1:400460/1‑251 (MQ=255)
                                                                                                                         gTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGAGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGtt                                                                                                                   >  1:211797/1‑251 (MQ=255)
                                                                                                                         gTCTGATCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCCGGCTACCTTCGGGGGCCCCATCCAACGGGTTCCCTTCGGTATCGGACAGCGACGGTTGAACGCCGAAGACCGGACGTGTTCCCGAACCGGCATTCAACTCGGTAGCGCCACGCAGCGGGGGAACGAAGAAACCACCGTGTTCGGTCTGCCCCCAGGtcacagacaacggccatttatcggt                                                                                                                   >  2:164863/1‑227 (MQ=255)
                                                                                                                           cTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGc                                                                                                                 >  1:115676/1‑251 (MQ=255)
                                                                                                                           cTGTTCAAAACGTTCGTGATCGCCAAACAGCGCACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACAAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGAACTTTCTCGTTGc                                                                                                                 >  2:349273/1‑251 (MQ=255)
                                                                                                                            tGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACAGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGcc                                                                                                                >  1:270141/1‑251 (MQ=255)
                                                                                                                                cAAAACGTCCGTGATCGCCAAACAGCGTCCGCGCCTGCCCCGGCCAGGAGTCGGTGATTCCCAGCCTCCCTTCGGTGGCCCCCTCCAGCGGGTTCCCTTCTTTATCGCCCAGCGTCGGTTGCACGCCGAAGAACGGACTTGTTGCCGAACCGGCTTTCAGCGCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGAtt                                                                                                            <  1:19013/251‑1 (MQ=255)
                                                                                                                                   aaCGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGAttttt                                                                                                         >  1:242301/1‑251 (MQ=255)
                                                                                                                                      gTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTc                                                                                                      >  1:215877/1‑251 (MQ=255)
                                                                                                                                       ttCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGACAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGCTCATGCAACCGCCGGTTTCGGTCTGCCACCGGGTATCGACCACCGTACATTTCCCGTTGCCGCTTTTTTTcc                                                                                                     >  2:20527/1‑251 (MQ=255)
                                                                                                                                         cGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCCGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCCCCAGGTATCGACCACCGGCCATTTCTCGTTGCCGATTTTTTTcccg                                                                                                   >  2:207511/1‑249 (MQ=255)
                                                                                                                                                   aaacggcgtccccgctggccccGGCCGGGAGGCGGTGATCACCAGGCGTCCTTCGGTGGCCCCCCCCGGCGGGTTACCTTCTTTATCGCCCAGCGCGGGTGGCACGCCGAAGAACGGACGGGTTGCGGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGGTATAATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTccc                                                                                         <  2:115676/232‑1 (MQ=255)
                                                                                                                                                   aaaCAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTccc                                                                                         >  1:130555/1‑251 (MQ=255)
                                                                                                                                                   aaaCAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTccc                                                                                         >  1:118529/1‑251 (MQ=255)
                                                                                                                                                         cGTACGCGCCTGCCCCGGCCAGGAGTCGGTGATTCCCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTACCGACCAGCGTCGGTGGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCtt                                                                                   <  1:86884/251‑1 (MQ=255)
                                                                                                                                                         cGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTATGCCACCAGGTATCGACCACCGGGCATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCtt                                                                                   >  1:355490/1‑251 (MQ=255)
                                                                                                                                                                    gccccGGCCAGGCGTCGGTGATTACCAGGTTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCTTTATCGCCCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCCCGGTAGCGCCAGGCAGCGGGTTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGGTCCACTCCCACGCTTCCGGGTTAAtt                                                                        <  1:461601/249‑1 (MQ=255)
                                                                                                                                                                                   cGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGGATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTAACACTCCCACGCTTCCGGGTTAATTGgcccgcccccagg                                                          >  2:422961/1‑245 (MQ=255)
                                                                                                                                                                                            ccAGGCTACCTTCGGTGGCCCCATCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTGGCCGAGTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTCATTGGCTCGACCCCGGGACCgggactg                                                >  2:269570/1‑245 (MQ=255)
                                                                                                                                                                                                   aCCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGACCCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAg                                           <  1:14434/249‑1 (MQ=255)
                                                                                                                                                                                                           gggCCCCCTCCAGCGGGTTCCCTGCGTTATCGTCCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTGGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGGGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAAcgg                                 <  1:349273/250‑1 (MQ=255)
                                                                                                                                                                                                            ggCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTTCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTAGGTAGCGCCAGGAAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGGATCGACCACCGGACCTTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCAc                                                                                       >  1:338617/1‑196 (MQ=255)
                                                                                                                                                                                                             gCCCCCTCCAGCGGGTTACCTTCGTTAGCGCCCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAAcggtc                               <  1:220033/251‑1 (MQ=255)
                                                                                                                                                                                                               ccccTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGTAGCGGGGGGATCATGAAACCGCCGGTTTCGGTCTGCCACCCGGTATCGACCACCGGCCATTTCTCGTGGCCGATTTTTTTCCAGAACCACTCCCAAGCTTCCGGGGTAATTGGCTCGCCCACGTGCCCGAGGCTGCACAGCGACGaagcggcgg                             >  2:413026/1‑244 (MQ=255)
                                                                                                                                                                                                                 ccTCCAGCGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGACCGGTCGGTg                           >  1:430808/1‑251 (MQ=255)
                                                                                                                                                                                                                        cGGGTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCg                     >  1:162748/1‑250 (MQ=255)
                                                                                                                                                                                                                         gggTTACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGGGATCATGAAACCGCCGGTTTCGGTATGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGGTTTTTTTTCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGGGACGCACGGTCGGGGCCTTCGa                    >  2:84401/1‑250 (MQ=255)
                                                                                                                                                                                                                             tACCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACgg                                                           >  1:137472/1‑207 (MQ=255)
                                                                                                                                                                                                                              aCCTTCGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGACCCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCGCGCCCACGGAACCGAGACTGCGCAGCGCCGAAAGGTCGGGGCCTTCGATAGCtt              >  2:121007/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gcGTTATGGACCAGCGTCGGTTGCCCGCGAAAGAAGGGAGGGGTGCCCGACCCGGCTTTCAGCGGGGTAGCGCCAGGCAGCGGGGTGATCAGGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCTTTAt           <  2:224094/249‑1 (MQ=255)
                                                                                                                                                                                                                                   cGTTATCGACCAGCGTCGGTTGCACGCCGAAGAACGGCCGTGTTGCCGACCCGGTTTTCACCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCTTTATCg         <  1:496603/251‑1 (MQ=255)
                                                                                                                                                                                                                                        tCGACCAGCGTCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCACCCACGGAACCGAGAATGCGCAGCGACGAACGGGCGGTGCCTTCGATCGCTTTATCGCCTTc    >  1:94634/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gACCAGCGTCGGTTGCACGCCGAAGCACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCCGCGACGAACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCg  >  2:374593/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
CAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCG  >  NC_000913/4285791‑4286275

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