Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,286,370 G→C P334P (CCC→CCG 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,3700GC97.1% 235.1 / ‑6.2 68P334P (CCC→CCGacsacetyl‑CoA synthetase (AMP‑forming)
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base C (34/32);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGCGGATCGCCGTGGGTGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTC  >  NC_000913/4286124‑4286618
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                          
ggTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTGATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGATCGGTCGGTGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTa                                                                                                                                                                                                                                                        <  1:174501/251‑1 (MQ=255)
     gCCACCAGGAAACGACCAACGGCCATTTCTCGTTGACGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGATGATTGGCTCGCCCACGGAAAAGAGACTGCGCAGCGACGATAGGGCGGTGCCTTCAATCGCTTTATACCCTTCCGCCATCAGCGCGCGAATCGCCGTGGACGCGGTATAGAGAATATTGAACAGAAGCTTGACCAACAACTGCGCCATACGGGCAAGCGTCGGCCAGCTCGGTACGcc                                                                                                                                                                                                                                                    >  2:191208/1‑250 (MQ=255)
       caccaGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTGATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGATCGGTCGGTGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTc                                                                                                                                                                                                                                                 <  1:149880/251‑1 (MQ=255)
       caccaGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTGATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGATCGGTCGGTGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTc                                                                                                                                                                                                                                                 >  1:312465/1‑251 (MQ=255)
            ggTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTGATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGATCGGTCGGTGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATCGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCCTACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAAcc                                                                                                                                                                                                                                            >  2:367888/1‑250 (MQ=255)
               aTCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTGATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGATCGGTCGGTGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCa                                                                                                                                                                                                                                         <  1:440521/251‑1 (MQ=255)
                   cccaccGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTGATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGATCGGTCGGTGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGt                                                                                                                                                                                                                                     <  2:418915/250‑1 (MQ=255)
                    ccaccGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGGTGATTGTCTCGCCCACTGAACCGAGAATGCGCAGCGACGATCGGTCGGTGCCTTCAATCGCTTTATCCCCTTCCGCCAGCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCGTGTCCACCACCTGCGCCATTCGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTg                                                                                                                                                                                                                                    <  1:308134/251‑1 (MQ=255)
                      accGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTGATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGATCGGTCGGTGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGt                                                                                                                                                                                                                                  >  1:477124/1‑251 (MQ=255)
                       ccGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTGATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGATCGGTCGGTGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTc                                                                                                                                                                                                                                 <  1:120528/251‑1 (MQ=255)
                       ccGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTGATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGATCGGTCGGTGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTc                                                                                                                                                                                                                                 >  1:222736/1‑251 (MQ=255)
                            cATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTGATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGATCGGTAGGTGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGAACGCCTTCAAACATCAGCGTGGTCGCAcc                                                                                                                                                                                                                            >  2:57241/1‑251 (MQ=255)
                                   cGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTGAGTGGCTCTCCCTCGGAACCGAGAATGCGCATCGACGATCGGTCGGTGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCTCGCGAATCGCCGTTGGCGCGGTATAGAGAATATTGACCTGAGGCTTGTCCACCACCTGCGCCATACGAGCATGCTACGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCTCCGCAg                                                                                                                                                                                                                        <  1:432765/248‑1 (MQ=255)
                                    gTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTGATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGATCGGTCGGTGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTAGAGAGACTATTGACCTGATGCTTGTCCACCACCTGCGCCATAAGTGCAGGCGTCGGCCAGTTCGGTACGATTTAAAACATCAGCGTGGTCGGACCGCAggac                                                                                                                                                                                                                     >  1:146509/1‑248 (MQ=255)
                                      tGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTGATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGATCGGTCGGTGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGc                                                                                                                                                                                                                  >  1:173497/1‑251 (MQ=255)
                                       gCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTGGTTGGCGCGCCCACGGAACCGAGAATGCGCAGCGACGATCGGTCGGTGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGc                                                                                                                                                                                                                  <  2:25121/250‑1 (MQ=255)
                                          gATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTGATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGATCGGTCGGTGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGc                                                                                                                                                                                                              >  1:330535/1‑251 (MQ=255)
                                          gATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTGATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGATCGGTCGGTGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGc                                                                                                                                                                                                              >  1:202/1‑251 (MQ=255)
                                                  cccAGTACCACTCCCCCGCGTCCTGGTTGATTGGCCCGCCCACGGAACCGAGAATGCGCAGCGACGATCGGTCGGTGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTAGAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGTTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGc                                                                                                                                                                                                      <  2:486784/250‑1 (MQ=255)
                                                                 aCGCTTCCGGGTTGATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGATCGGTCGGTGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCg                                                                                                                                                                                       <  1:149469/251‑1 (MQ=255)
                                                                         gggatgataggcacacccACGGCCCCGAGAATGCGCAGCGACGATCGGTCGGTGCCTTCCATCGCTTTATCCCCTTCCGCCATCAGCGCGCGCATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCCGGCGTCGGCCAGCTCGGTACGCCTTCAACCAGCAGCGTGGGCGCACCGCAGGCAAGCGCGCCGTACAGCACGTAACCGTGTCCGGCCAccc                                                                                                                                                                                >  2:181600/16‑250 (MQ=255)
                                                                             tgaTTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGATCGGTCGGTGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGcc                                                                                                                                                                            <  1:209174/248‑1 (MQ=255)
                                                                                                   gagaATGCGCAGCGACGATCGGTCGGTGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGGGCACCAGt                                                                                                                                                       >  1:257515/1‑248 (MQ=255)
                                                                                                         gcgcAGCGACGATCGGTCGGTGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAAGCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGTGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCa                                                                                                                                               <  2:182825/250‑1 (MQ=255)
                                                                                                         gcgcAGCGACGATCGGTCGGGGCCGTCAATCGCTTTATCCCCGTCCGCCATCAGCGCGCGAATCGCCGTGTGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGGCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCAc                                                                                                                                              <  2:40825/251‑1 (MQ=255)
                                                                                                          cgcAGCGACGATCGGTCGGTGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCCGCGTGGTCGCACCGCCGGCCAGCGGGCCGTACAGCGAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGGGCAcgagaggaggagggg                                                                                                                                             >  2:399698/1‑237 (MQ=255)
                                                                                                          cgcAGCGACGATCGGTCGGTGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGGCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACGGTGTCCGGTCACCCAGCCCACATCGGCGGTGCa                                                                                                                                                            <  1:399698/236‑1 (MQ=255)
                                                                                                               cgacgaTCGGTCGGTGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGgat                                                                                                                                        <  1:450786/251‑1 (MQ=255)
                                                                                                                 acgaTCGGTCGGTGCCTGCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGgatga                                                                                                                                      <  2:127876/251‑1 (MQ=255)
                                                                                                                  cgatcggtcggtGCCTTCCATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAACCATCCGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCACGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCCCCCGTAGATATCCCCCggcttct                                                                                                                                     >  2:182088/5‑246 (MQ=255)
                                                                                                                  cgatcggtcggtGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGgatgat                                                                                                                                     <  1:335907/247‑1 (MQ=255)
                                                                                                                  cgatcggtcggtGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGGTCACCAAGAGATATCACCCGgatgat                                                                                                                                     >  2:355114/5‑251 (MQ=255)
                                                                                                                   gatcggtcggtGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGGAGGTATCACCCGgatgat                                                                                                                                     >  1:112260/4‑250 (MQ=255)
                                                                                                                      cggtcggtGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATc                                                                                                                                 <  1:341882/251‑1 (MQ=255)
                                                                                                                        gtcggtGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCaa                                                                                                                               <  1:419396/251‑1 (MQ=255)
                                                                                                                          cggtGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGCATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGAAGATATCACCCGGCTGATAATCaaa                                                                                                                              >  1:169833/1‑250 (MQ=255)
                                                                                                                          cggtGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAg                                                                                                                             <  1:271329/251‑1 (MQ=255)
                                                                                                                          cggtGCCTTCAATCGCTTTATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAg                                                                                                                             <  1:49275/251‑1 (MQ=255)
                                                                                                                                           ttATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGTTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATTCGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGTTCGCACCTCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGTTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTc                                                                                                             <  1:384415/249‑1 (MQ=255)
                                                                                                                                           ttATCCCCTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTc                                                                                                             <  1:388795/249‑1 (MQ=255)
                                                                                                                                             atcccctcccGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCCCCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCAgcg                                                                                                         <  2:73726/243‑1 (MQ=255)
                                                                                                                                              tcccctcccGCCATCAGCTCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATTCGGGCAGGCGTCGGCCAGGTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCCGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCAgcg                                                                                                         <  1:19367/243‑1 (MQ=255)
                                                                                                                                               ccccTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCAGcgcc                                                                                                       <  1:7027/249‑1 (MQ=255)
                                                                                                                                                 ccTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCAGcgccgc                                                                                                     >  1:82603/1‑251 (MQ=255)
                                                                                                                                                  cTTCCGCCATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATAAGGGCAGGCATCGGCCCGTTCGGTACGCCTTCAAACATCAGCGTGGACGCACAGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGCCCGGTCACCCAGCCCACATCGGAGGTGCACCCGTAGATATCACCCGGCTGATACTCAAAAACCTATTTAAAGGTCAGCGCcgcg                                                                                                    >  1:461310/1‑251 (MQ=255)
                                                                                                                                                     cctccATCAGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGAGATCACCCGGATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTac                                                                                                 <  2:406602/248‑1 (MQ=255)
                                                                                                                                                             aGCGCGCGAATCGCCGTGGGCGCGGTATAGCGAATATTGACCTGATGCTTGTCCACCACCTGCGCAATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCCGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATTTCACCCGGATGATAATCCAAGACATATTTAAAGGt                                                                                                              >  1:225185/1‑230 (MQ=255)
                                                                                                                                                             aGCGCGCGAATCGCCGTGGGCGCGGTATAGAGAATATTGTCCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGTGCCGTACAGCAAGTAACGGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGt                                                                                                              <  2:225185/230‑1 (MQ=255)
                                                                                                                                                               cgcgcgAATCGCCGGGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTGCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAAc                                                                                       <  2:247540/251‑1 (MQ=255)
                                                                                                                                                                  gcgaaTCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCCGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAAccgc                                                                                    >  1:183318/5‑251 (MQ=255)
                                                                                                                                                                   cgaaTCGCCGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCCCCCCCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAAccgcc                                                                                   <  1:299675/248‑1 (MQ=255)
                                                                                                                                                                          ccGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTCCGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCCGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGgtagt                                                                             >  1:474811/1‑250 (MQ=255)
                                                                                                                                                                          ccGTGGGCGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGgtagt                                                                             >  1:188198/1‑250 (MQ=255)
                                                                                                                                                                            gTGGGCGCGGTATAGAGAATCTTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCCGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCCGCCCACATCGGCGGTGCCCCAGTAGATATCCCCCGGATGATAATCCAAGACCATTTTAAACGTCAGCGCCGCGAACACACGATAACCACCTgtagac                                                                            >  2:476010/1‑247 (MQ=255)
                                                                                                                                                                                     gTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACCTATTTAAAGGGCAGCGCCGCGTACACCAGATAAACGCCGGTAGAATGCAGCCCAcc                                                                 >  2:142338/1‑251 (MQ=255)
                                                                                                                                                                                              atatTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTCGATATCACCCGGATGATAATCACAGACCTATTTTAAGGTCAGCCCCGCCTACACCAAATTACCGCCCGTAGTATGCCGCACCCCTTTTGGCtt                                                        >  2:169343/1‑251 (MQ=255)
                                                                                                                                                                                                atTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTAc                                                      >  1:364560/1‑251 (MQ=255)
                                                                                                                                                                                                      cTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGt                                                  >  1:392095/1‑249 (MQ=255)
                                                                                                                                                                                                      cTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGGCGCACCGCAGGCCAGCGGGCCGGACAGCAAGTAACTGTGTCCGGTGACCCAGCCCACAGCGGCGGTGCACCAGTAGAGATCACCCg                                                                                                                                           <  2:44014/161‑1 (MQ=255)
                                                                                                                                                                                                      cTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGAGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCg                                                                                                                                           >  1:44014/1‑161 (MQ=255)
                                                                                                                                                                                                            cTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTCGGAACGCCTTCCAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACCTCGGCGGTGCACCAGTAGCTATCACCCGGCTGATAATCAACGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGAATGCAGCACACCTTTTGGCTTACCGGTAGCACCgg                                          >  1:306065/1‑251 (MQ=255)
                                                                                                                                                                                                               gccccccccccgcgcCATAGGGGCGGGCGGCGGCCAGTGGGGTACGCCTACAAACATCAGCGTGGTCGCCCCGCGGGCCAGCGGGCCGTACAGAAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATACCACCCGGATGATAATAAAAGACATATTAAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAg                                        <  1:355114/239‑1 (MQ=255)
                                                                                                                                                                                                                               aTACGGGCAGGCGTCGGCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGGGTAGAGAATAAAc                         >  1:27908/1‑249 (MQ=255)
                                                                                                                                                                                                                                    ggCAGGCGTCGGCCCGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCCGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAAGTCAGCGCCGCGTACACCAGATAAACGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGGGTAGAGAATAAACAGCgg                    >  2:142949/1‑249 (MQ=255)
                                                                                                                                                                                                                                        ggCGTCGGCCAGTCGGGTCCGCCTTCAAACATCAGCGTGGCCGCCCCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCCCCCGGCCCACACCGGGGGTCCACCAGTAGATATCCCCCGTATGATAATAAAAAACATATTTAAAGGTCAGCGCCGTGTCCACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTc              <  1:57241/251‑1 (MQ=255)
                                                                                                                                                                                                                                             cGGCCAGTGCGGTCCGCCTCAAAACATCGGCGGGGCCGCACCGCAGGCCAGCGGGCCGTCCAGCAAGTAACTGTGTCCGGTCCCCCAGCCCACATCGGCGGTCCACAATTAAATAGCACCCGGATGATAATCAAAGAAATATTTAAAGGTCAGCGCCGCGTACCCCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCg          <  1:367406/250‑1 (MQ=255)
                                                                                                                                                                                                                                              ggCCAGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCCGTAGCTATCACCCGGATGATAATCAAAGACATATTTAAAGGTCAGCGACGCGGACACCAGATAACCGCCGGTAGTATGCAGCACCCCTTTTGGCTTACCGGTAGAAACGGAGGTGTAGAGAATAAACATCGGATCTTCGGCGtt        >  2:386337/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  aGTTCGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCCCCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCCGCCCACCTTTTGGCATACCGGTAGAACCGGCGGTGTAGGGAATAAACAGCGGATCTTCGGCGATCAtt    >  2:319579/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     ggggTACGCCTTCAAACACCAGCGGGGCCGCACCGCGGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGTTCACCCAGCCCACATCGGGGTGGCCCCAGTAGATATCCCCCGGAGGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGGTCAtctc  <  1:355054/249‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                          
GGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCTTTATCGCCTTCCGCCATCAGCGCGCGGATCGCCGTGGGTGCGGTATAGAGAATATTGACCTGATGCTTGTCCACCACCTGCGCCATACGGGCAGGCGTCGGCCAGTTGGGTACGCCTTCAAACATCAGCGTGGTCGCACCGCAGGCCAGCGGGCCGTACAGCAAGTAACTGTGTCCGGTCACCCAGCCCACATCGGCGGTGCACCAGTAGATATCACCCGGATGATAATCAAAGACATATTTAAAGGTCAGCGCCGCGTACACCAGATAACCGCCGGTAGTATGCAGCACACCTTTTGGCTTACCGGTAGAACCGGAGGTGTAGAGAATAAACAGCGGATCTTCGGCGTTCATCTC  >  NC_000913/4286124‑4286618

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