Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,286,001 (C)5→6 coding (1371/1959 nt) 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,0011.C87.2% 132.8 / 15.2 47coding (1371/1959 nt)acsacetyl‑CoA synthetase (AMP‑forming)
Reads supporting (aligned to +/- strand):  ref base . (2/4);  new base C (26/15);  total (28/19)
Fisher's exact test for biased strand distribution p-value = 2.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGG  >  NC_000913/4285753‑4286248
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                        
ggATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCTTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGcccccct                                                                                                                                                                                                                                                         <  1:236475/251‑3 (MQ=255)
 gATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGcccccctc                                                                                                                                                                                                                                                        <  1:316369/251‑4 (MQ=255)
    gcgcCACCCGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCGCGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGcccccctcca                                                                                                                                                                                                                                                      <  2:211915/250‑6 (MQ=255)
    gcgcCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGcccccctcc                                                                                                                                                                                                                                                       >  2:112701/1‑245 (MQ=255)
           caGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGtt                                                                                                                                                                                                                                              >  2:87530/1‑250 (MQ=255)
                   aCTCAATCTCTGCCCTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACCCGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGtt                                                                                                                                                                                                                                     <  1:215485/251‑1 (MQ=255)
                      cAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATc                                                                                                                                                                                                                                  >  1:178474/1‑251 (MQ=255)
                             tGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAgc                                                                                                                                                                                                                           >  1:17090/1‑251 (MQ=255)
                                           ggTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGcc                                                                                                                                                                                                             >  1:63026/1‑251 (MQ=255)
                                           ggTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGcc                                                                                                                                                                                                             >  1:285967/1‑251 (MQ=255)
                                               aCCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACTAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCgaag                                                                                                                                                                                                         <  2:201039/251‑1 (MQ=255)
                                                ccGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTACGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCgaaga                                                                                                                                                                                                        <  2:253473/251‑1 (MQ=255)
                                                   gggaCGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGAGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTGATCGACCAGCGCCGGTTGCACGCCGAAGAACg                                                                                                                                                                                                     <  2:256566/249‑1 (MQ=255)
                                                              cACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGcc                                                                                                                                                                                          <  1:517/251‑1 (MQ=255)
                                                                cgtcgtCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGAGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCGCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGa                                                                                                                                                                                        <  2:110356/251‑1 (MQ=255)
                                                                  tcgtcCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAAc                                                                                                                                                                                      >  1:92504/1‑251 (MQ=255)
                                                                       cacaCGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCt                                                                                                                                                                                 >  1:109627/1‑251 (MQ=255)
                                                                                        aGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGcc                                                                                                                                                                >  1:274116/1‑251 (MQ=255)
                                                                                           aaTAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAgg                                                                                                                                                             >  1:329247/1‑251 (MQ=255)
                                                                                           aaTAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAgg                                                                                                                                                             >  1:84150/1‑251 (MQ=255)
                                                                                                 cATCTGCATCGCGACGCGCGCCGGCGCCGCTGAAAGACATTTTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAGCGTTCGTGAGCGCCAGACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCGCCCCGCGAGCGGGTTACCTTCGTTATCGTCCAGCGCCGGTTGCACGCTGAAGAACGGACGGGTGGCCGAACCTGCTTTCAGCGCGGTAGCGCAAGGCAGCgg                                                                                                                                                       <  2:338534/251‑1 (MQ=255)
                                                                                                  aTCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCggg                                                                                                                                                      >  1:302847/1‑251 (MQ=255)
                                                                                                      tCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCggg                                                                                                                                                      <  1:175749/247‑1 (MQ=255)
                                                                                                            cgACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGa                                                                                                                                            <  2:149329/251‑1 (MQ=255)
                                                                                                                          gccgcTAAAATACATATTTTTGAAGGGGGAGAAGTAGCTCTGTTCAAAACGTTCGTGATCGCCAAACGGCGTACGCGCCTGCCCCGGCCAGGATTCGGTGATTACCAGGCTACCTTCGGTGCCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCCCGCCGAAGAACGGACGTTTTGCCGAACCGGCTTTCAGCGCGTTAGCGCCAGGCAGCGGGGTTATCATGAAACCGCCGGTTTCg                                                                                                                              <  2:169990/251‑1 (MQ=255)
                                                                                                                                  aaTACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGcac                                                                                                                      >  2:15931/1‑249 (MQ=255)
                                                                                                                                      cATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCACACAGCGTACGCGCATGACCCGGCCAGGAGTCGGAGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCCTGAAACCGCCGGTTTCGGTCTGccacc                                                                                                                    >  1:346704/1‑249 (MQ=255)
                                                                                                                                       atatTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAACCGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTGGCACGCCGAAGAACGGACGGGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCATGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAgg                                                                                                                 <  1:21351/251‑1 (MQ=255)
                                                                                                                                          tttttGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTAt                                                                                                              >  1:30186/1‑251 (MQ=255)
                                                                                                                                           ttttGAAGGTGGAGAAGTAGGTCTTTTAAAACCGTTCGTGATCGCCAAACAGCGTCCGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTCCCTTCTTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGACCCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTAt                                                                                                              <  2:317582/250‑1 (MQ=255)
                                                                                                                                                 aGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGaccac                                                                                                       >  1:53031/1‑251 (MQ=255)
                                                                                                                                                 aGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGaccac                                                                                                       >  1:178938/1‑251 (MQ=255)
                                                                                                                                                  ggTGGAGAGGTGGGTCTGTTAAAACCGTCCGTGAGCGCCAAACGGCGTACGCGCCTGCCCCGGCCAGGGGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGCCCAGCGCGGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGaccacc                                                                                                      <  2:64200/251‑1 (MQ=255)
                                                                                                                                                    tGGAGAAGAAGGTATGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCACTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCgg                                                                                                    >  1:312184/1‑251 (MQ=255)
                                                                                                                                                       agaAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGTGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCAAGAAACCGCCGGTTTCGGTCTGCCACCAGGAATCGACCACCGGaa                                                                                                  >  1:96472/1‑249 (MQ=255)
                                                                                                                                                             aGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTCCCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTtctc                                                                                           <  1:119754/251‑1 (MQ=255)
                                                                                                                                                                   gTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGCGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGcc                                                                                     >  1:374584/1‑251 (MQ=255)
                                                                                                                                                                                gTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTcc                                                                          >  2:328122/1‑249 (MQ=255)
                                                                                                                                                                                 tGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGt                                                                       >  1:333399/1‑251 (MQ=255)
                                                                                                                                                                                       caaaaCAGCGTGCGCGCCTGCCCCGGCCAGGGGTCGGTGATTACCAGGCTCCCTGCGGTGGCCCCCCCCCAGGGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGACCGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACt                                                                 <  2:109627/249‑1 (MQ=255)
                                                                                                                                                                                                      gccggccccGGCCGGGAGTCGGTGATTCCCAGGCTACCTGCGGGGGCCCCCCCCCAGCGGGTTACCTTCGTTATCGCCCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCCCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGt                                                  <  2:108013/245‑1 (MQ=255)
                                                                                                                                                                                                      gcCTGACCCGGCCCGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCgggg                                                  >  2:211334/1‑250 (MQ=255)
                                                                                                                                                                                                      gcCTGAACCGGACAGGAGTAGGTGATTACCAGGCTACCTTCGGTGGCACCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGt                                                  >  2:336072/1‑251 (MQ=255)
                                                                                                                                                                                                          gccccGGCCAGGAGTCGGTGATTCCCCGTCTCCCTTCGGTGGCCCCCCTCCAGGGGGTTCCCTTCGTTATCGCCCAGCGCCGGTTGCACGCCGAAGACCGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTa                                                <  2:96156/247‑1 (MQ=255)
                                                                                                                                                                                                          gACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCATCGCCGGTTGCACGCCGCAGCACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCCGCGGGGTGATCATGACACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGCCATTTCTCGTTGCCGATTTTTTTCCCGTACCACTCCCACGCTTCCGGGttcat                                              >  2:191344/1‑248 (MQ=255)
                                                                                                                                                                                                                gCCAGGGGGCGGTGATTACCAGGCTACCTGCGGTGGCCCCCCTCAAGGGGGTTACCTTCGTTAGCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGGAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTc                                        <  2:269429/251‑1 (MQ=255)
                                                                                                                                                                                                                 ccAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCCGTACCACTCCCACGCTTCCGGGTTAATTGGCTCg                                       >  1:183383/1‑251 (MQ=255)
                                                                                                                                                                                                                            tGATTACCAGGCTACCTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAAcc                            >  2:291768/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aCCAGGCTCCCTGCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTAGCGACCAGCGCCGGTGGCACGACGAAGAACGGGCGTGTTGCCGAACGGGTTTTCAGCTCGGTAGCGCCAGGCAGCGGGGGGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGaa                       <  2:63026/251‑1 (MQ=255)
                                                                                                                                                                                                                                           cTTCGGTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGCACCGAGAATGCGCAGcg              >  2:360651/1‑250 (MQ=255)
                                                                                                                                                                                                                                               ggTGGCCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGACCCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGGGAATGCGCAGcgacga          >  2:183592/1‑250 (MQ=255)
                                                                                                                                                                                                                                                gggccccccTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGAACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCGCGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAAcg       <  1:231013/247‑1 (MQ=255)
                                                                                                                                                                                                                                                   ccccccTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAAcggtc    <  1:229388/250‑1 (MQ=255)
                                                                                                                                                                                                                                                    cccccTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAAcggtcg   >  1:252035/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      ccccTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAAcggtcgg  >  1:234911/1‑251 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                        
GGATGCGCCACCAGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACGTTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGG  >  NC_000913/4285753‑4286248

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