Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,564,385 C→T S102N (AGC→AAC)  yjiL ← putative ATPase, activator of (R)‑hydroxyglutaryl‑CoA dehdratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,564,3850CT98.5% 226.8 / ‑4.7 65coding (305/768 nt)yjiLputative ATPase, activator of (R)‑hydroxyglutaryl‑CoA dehdratase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (33/31);  minor base . (0/1);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

TGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGCTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTGCTGTTGCCGGGCGAAAGGGGGTTGGAACGAGGAAACGGCGCGTAATCACGCCG  >  NC_000913/4564141‑4564627
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
tGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCt                                                                                                                                                                                                                                              >  2:375444/1‑251 (MQ=255)
       aTCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGaccacc                                                                                                                                                                                                                                       <  2:101158/251‑1 (MQ=255)
         cGCTTCGGGCGCTGCGCCCGCTTGGCGCAGGTTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGTATCTCTTTGTTGCCCTGACCACCg                                                                                                                                                                                                                                      <  2:173815/250‑1 (MQ=255)
          gCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGAt                                                                                                                                                                                                                                    <  1:375444/251‑1 (MQ=255)
          gCTTATGGCGCGACGCCCGCGGAGCGCAGGCTGATCGCTTCGGATTCAGCAACCACTGTGCATAGACTCGTGATGGCGTGCGGTGTGACATTTTCGTTAATGCTGTCGAGTTGCGCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATAAAGCTGATTCACTTTTTTGTCCTGACCACCGAt                                                                                                                                                                                                                                    <  2:421114/251‑1 (MQ=255)
          ccTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGTTCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATTCTGTCGAGTTGCCCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTTCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTTTTGTCCTGACCACCGAt                                                                                                                                                                                                                                    <  2:384859/250‑1 (MQ=255)
                ggCGCGAGGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCg                                                                                                                                                                                                                                <  2:344078/249‑1 (MQ=255)
                          ccGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCTTGCGCTAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACcgcgcgc                                                                                                                                                                                                                    <  1:56290/251‑1 (MQ=255)
                                aGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCt                                                                                                                                                                                                              >  1:351857/1‑251 (MQ=255)
                                     aGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGCAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGgcgc                                                                                                                                                                                                         >  1:149557/1‑251 (MQ=255)
                                      ggCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCa                                                                                                                                                                                                        <  1:38831/251‑1 (MQ=255)
                                      ggCTGATCGCATCTGATTCAGCACACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGCGTTGCTCGACGCTGGTGCCAGGCGTGCGCGAGATCACCTCCAGGAAACGACCGGTGCCCGCCGCGCATTTGTCATTCATCAGGACATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCTATTACCTCGCGCGTTGCTGgcgc                                                                                                                                                                                                         >  2:7972/1‑250 (MQ=255)
                                            acGCTTCTTATTCAGCAAAAACTGTGCACTTATGCGTGAGGGCGTGCGGTGTGACATTTTCGGGAATGCGGTCTGGTTGCTCGTCGGTGGTCCCAAGCGTGAGCTTGATCACCTCCAGTAAACGCCTGGTGCCCGCCTCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGACTCTCTTTGTTGTCCTGACCACCGATGTCGATTCCCGCGCTCGTGGCTTGCGCAAGaaa                                                                                                                                                                                                   <  2:394998/249‑1 (MQ=255)
                                             cGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAAcc                                                                                                                                                                                                 <  2:108975/251‑1 (MQ=255)
                                                 tCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTgc                                                                                                                                                                                             >  2:356396/1‑251 (MQ=255)
                                                   tgcttCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTgcgc                                                                                                                                                                                           >  1:84130/4‑251 (MQ=255)
                                                              acacTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGa                                                                                                                                                                                <  2:27253/251‑1 (MQ=255)
                                                                   gtgCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCACGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAgg                                                                                                                                                                            >  1:290584/1‑250 (MQ=255)
                                                                       ccaTACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCGCTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGAtt                                                                                                                                                                       <  1:435190/250‑1 (MQ=255)
                                                                        caTACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGAttt                                                                                                                                                                      <  1:107652/251‑1 (MQ=255)
                                                                           aCTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGccc                                                                                                                                                                                    >  2:21088/1‑234 (MQ=255)
                                                                           aCTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGccc                                                                                                                                                                                    <  1:21088/234‑1 (MQ=255)
                                                                                 gATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGACACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCt                                                                                                                                                             >  1:155339/1‑251 (MQ=255)
                                                                                 gATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCGCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACATGTTACCGTCATCATCAAGCTTAATCACTTTGTTGTTCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTAcc                                                                                                                                                              <  1:191270/250‑1 (MQ=255)
                                                                                  aTGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGAGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGGGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGGGATTTCCGTTACCTg                                                                                                                                                            >  2:2648/1‑251 (MQ=255)
                                                                                      cGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGtttt                                                                                                                                                        >  1:235260/1‑250 (MQ=255)
                                                                                      cGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCAAGGAAACGCCCGGTGCACGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGCGATTTCCGTTACCTGttt                                                                                                                                                         >  1:420265/1‑250 (MQ=255)
                                                                                            gcgTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCTCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCa                                                                                                                                                  <  1:418385/251‑1 (MQ=255)
                                                                                                 aCATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTAACGCGCGCGTTGCTCGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATc                                                                                                                                             >  2:147740/1‑251 (MQ=255)
                                                                                                        cGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGAGTACCGCGCGCGGTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGGGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGt                                                                                                                                      <  1:211942/251‑1 (MQ=255)
                                                                                                               gCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTgccgcc                                                                                                                               <  1:39452/251‑1 (MQ=255)
                                                                                                                 tGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCg                                                                                                                             >  1:152842/1‑251 (MQ=255)
                                                                                                                  gTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCg                                                                                                                             <  2:235260/250‑1 (MQ=255)
                                                                                                                          gCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGCAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCAGTTACCTGTTGATCGGCAAAATCCACCAGTCGCCGCCCGTAGCAGGTg                                                                                                                    >  2:163546/1‑251 (MQ=255)
                                                                                                                           cTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCAGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCggcg                                                                                                                    >  1:466082/1‑248 (MQ=255)
                                                                                                                               aCGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGt                                                                                                               <  1:75793/251‑1 (MQ=255)
                                                                                                                                cgcgggTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACCCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTAGCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTc                                                                                                              <  2:38170/247‑1 (MQ=255)
                                                                                                                                cGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTCCCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTc                                                                                                              >  1:380252/1‑251 (MQ=255)
                                                                                                                                 gCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCa                                                                                                             <  1:241635/251‑1 (MQ=255)
                                                                                                                                  cTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGTAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGCTGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGACCGTc                                                                                                              >  1:211240/1‑249 (MQ=255)
                                                                                                                                              gTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGt                                                                                                >  1:158443/1‑251 (MQ=255)
                                                                                                                                                             cccAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTcgcg                                                                                 <  2:276783/250‑1 (MQ=255)
                                                                                                                                                              ccAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGAGTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTCCCGTTACCTGTTTATCGGCAAAATCCTCCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTcgcg                                                                                 <  1:40256/250‑1 (MQ=255)
                                                                                                                                                                 ggAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCgcgc                                                                                >  1:78052/1‑248 (MQ=255)
                                                                                                                                                                 ggAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCag                                                                              <  2:78052/250‑1 (MQ=255)
                                                                                                                                                                            gTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGTCAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGc                                                                  >  2:375844/1‑251 (MQ=255)
                                                                                                                                                                             tGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACATGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCt                                                                 >  1:69275/1‑251 (MQ=255)
                                                                                                                                                                                 cgccgcGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTc                                                               >  2:255244/1‑249 (MQ=255)
                                                                                                                                                                                  gccgcTCATTTGTCATTCAGCAGGAAATCGCACAGGTTACCGTCATCAGCAATCTGAGTCACTTTGTGTTCATGACCACCGATGTCGATTACCGCGCGCGGTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGAGTTCCGTTTCCTGTTTATCGTCAAAATCCACCAGTTGCCGCCCGTAGCCGTTGAGCGTCAGAAACGGCGTTGTCTCGAACTCTTCGCGCAGAGTTTCCCAGGCTTCAGt                                                            <  2:399782/251‑1 (MQ=255)
                                                                                                                                                                                       gcATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGGCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTg                                                       <  2:351857/251‑1 (MQ=255)
                                                                                                                                                                                            tGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGGTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGCGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTGCTGt                                                   <  1:465024/250‑1 (MQ=255)
                                                                                                                                                                                             gTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTGCTGTTg                                                 <  1:406871/251‑1 (MQ=255)
                                                                                                                                                                                                attcatCAGGAAATCGCACAGGTTACCGCCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGATGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTGCTGTTGCCg                                              >  1:17816/1‑251 (MQ=255)
                                                                                                                                                                                                        ggAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGCGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTCATTGCTGTTGCCGGGCGaa                                        >  2:139032/1‑249 (MQ=255)
                                                                                                                                                                                                         gAAATCGCACAGGTTACCGTAATCATCAAGCTGACTCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTGCTGTTGCCGGGCGAAAgg                                     >  1:143748/1‑251 (MQ=255)
                                                                                                                                                                                                          aactCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTACTTGCTGTTGCCGGGCGAAAggg                                    >  2:254609/4‑251 (MQ=255)
                                                                                                                                                                                                              cGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGACAGAAACGTCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTGCTGTTGCCGGGCGAAAGGGGGtt                                >  1:368205/1‑251 (MQ=255)
                                                                                                                                                                                                                  caGGTTACCGTCATCATCAAGCTGACTCACTTTGTTGTCCTGACCACCGCTGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTGCTGTTGCCGGGCGAAAGGGGGTTGGa                             >  2:315599/1‑250 (MQ=255)
                                                                                                                                                                                                                   aGGTTACCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCGCCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTGCCCAGGCTTCAGTAATTTCGGTTGCGCGGCGAAAGGGGGTTGGAAc                           <  2:290584/251‑1 (MQ=255)
                                                                                                                                                                                                                        aCCGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTGCTGTTGCCGGGCGAAAGGGGGTTGGAACGAgg                       >  2:212849/1‑250 (MQ=255)
                                                                                                                                                                                                                         ccGTCATCATCAAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTGTCCCAGGCTTCAGTAATTGCTGTTGCCGGGCGAAAGGGGGTTGGAACGAGGaa                     <  2:155339/251‑1 (MQ=255)
                                                                                                                                                                                                                           gTCATCATCAAGCTGAATCGCTTTGTTGTCCTGTCCACCGATGGCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTTCCTGTTTATCGGCAAAATCCGCCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGTCGTTGTCTCGAACCCTTCGCGCAGAGTTTCCCATGCTTCAGTAATTGCTGTTGCCGGGCGAAAGGGGGTTGGAACGAGGAAAc                   <  1:455168/251‑1 (MQ=255)
                                                                                                                                                                                                                                   caAGCTGAATCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTGCTGTTGCCGGGCGAAAGGGGGTTGGAACGAGGAAACGGCGCGt            >  1:127983/1‑250 (MQ=255)
                                                                                                                                                                                                                                            tCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGCGTTTCCCAGGCTTCAGTAATTGCTGTTGCCGGGCGAAAGGGGGTTGGAACGAGGAAAAGGAGCGTAATCACGCCg  >  2:160631/1‑251 (MQ=255)
                                                                                                                                                                                                                                            tCACTTTGTTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGCGCGCCCAGCCCGTGACAGGCGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTACTTGCTGTTGCCGGGCGCAAGGGGGTTGGACCGAGGAAACGGCGCGTAATCccgccg  >  2:127526/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
TGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCAAGCTGAATCACTTTGCTGTCCTGACCACCGATGTCGATTACCGCGCGCGTTGCTGGCGCAAGAAACCGTGCGCCCAGCCCGTGACAGGAGATTTCCGTTACCTGTTTATCGGCAAAATCCACCAGTTGCCGCCCGTAGCCGGTGAGCGTCAGAAACGGCGTTGTCTCTAACCCTTCGCGCAGAGTTTCCCAGGCTTCAGTAATTGCTGTTGCCGGGCGAAAGGGGGTTGGAACGAGGAAACGGCGCGTAATCACGCCG  >  NC_000913/4564141‑4564627

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: