Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 808,808 G→A G150G (GGC→GGT bioA ← adenosylmethionine‑8‑amino‑7‑oxononanoate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913808,8080GA98.0% 168.2 / ‑2.9 51coding (450/1290 nt)bioAadenosylmethionine‑8‑amino‑7‑oxononanoate aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (21/29);  minor base . (1/0);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

TTTTAACCATTCCGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCACCGCCGCGATTTCATGACGATGCGCCGCCATCAGGCGGGCAAAGCCCACCATATCGCGCTCATCCCATTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCGCCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCG  >  NC_000913/808571‑809049
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  
ttttAACCATTCCGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCACCGCCGCGATTTCATGACGATGCGCCGCCATCAGGCGGGCAAAGCCCACCATATCGCGCTCATCCCCTTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCCCAGACTGTGCATTTCGTTATCCGGATCGCACACCGACATCGCACCAAATGTATCGc                                                                                                                                                                                                                                      >  2:175502/1‑251 (MQ=255)
        aTTCCGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCACCGCCGCGATTTCATGACGATGCGCCGCCATCAGGCGGGCAAAGCCCACCATATCGCGCTCATCCCATTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATa                                                                                                                                                                                                                               >  1:86879/1‑250 (MQ=255)
          tCCGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCACCGCCGCGATTTCATGACGATGCGCCGCCATCAGGCGGGCAAAGCCCACCATATCGCGCTCATCCCATTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGGCTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAAcc                                                                                                                                                                                                                            >  1:209969/1‑251 (MQ=255)
          cccGGATGATACTCGCGCACCCCGCCTGCGCCCTGGACATTCGGCTCAAGGATCCCCGCCGCGATTTCATGACGATGCGCCGCCATCAGGCGGGCAAAGCCCACCATATCGCGCTCATCCCATTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAAc                                                                                                                                                                                                                             <  2:42213/249‑1 (MQ=255)
                              cccGCCTGCGCCCTGGACAATCGGCTCAATGATCACCGCCGCGATTTCATGACGATGCGCCGCCATCAGGCGGGCAAAGCCCACCATATCGCGCTCATCCCATTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAAc                                                                                                                                                                                                        >  1:138727/1‑251 (MQ=255)
                              cccGCCTGCGCCCTGGACAATCGCCTCAAGGACCACCGCCGCGATTTCATGACGATGCGCCGCCATCAGGCGGGCAAAGCCCACCATATCGCGCTCATCCCATTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAAc                                                                                                                                                                                                        <  2:3376/251‑1 (MQ=255)
                               ccGCCTGCGCCCTGGACAATCGGCTCAATGATCACCGCCGCGATTTCATGACGATGCGCCGCCATCAGGCGGGCAAAGCCCACCATATCGCGCTCATCCCATTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACg                                                                                                                                                                                                       >  2:164582/1‑251 (MQ=255)
                                cGCCGGCGCCCTGGACAATCGGCTAAATGATCTCCGCCGCGATTTCATGACGATGCGCCGCCATCAGGCGGGCAAAGCCCACCATATCGCGCTCATCCCATTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGc                                                                                                                                                                                                      <  1:364847/251‑1 (MQ=255)
                                  ccTGCGCCCTGGACAATCGGCTAAATGATCACCGCCGCGATTTCATGACGATGCGCCGCCATCAGGCGGGCAAAGCCCACCATATCGCGCTCATCCCATTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTg                                                                                                                                                                                                    <  2:61506/251‑1 (MQ=255)
                                       gcCCGGGACAAGCGGCACAATGATCACCGCCGCGATTTCATGACGATGCGCCGCCATCAGGCGGGCAAAGCCCACCATATCGCGCTCATCCCATTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGgcgcg                                                                                                                                                                                               <  1:46064/251‑1 (MQ=255)
                                             gACAATCGGCTCAATGATCACCGCCGCGATTTCATGACGATGCGCCGCCATCAGGCGGGCAAAGCCCACCATATCGCGCTCATCCCATTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGc                                                                                                                                                                                         >  1:382863/1‑251 (MQ=255)
                                               aaaacGGCTCAATGAACACCGCCGCGATTTCATGACGATGCGCAGCCATCAGGCGGGCAAAGCCCACCATATCGCGCTCATCCCATTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCt                                                                                                                                                                                       <  2:277288/247‑1 (MQ=255)
                                                   gggCTCAATGATCACCGCCGCGATTTCATGACGATGCGCCGCCATCAGGCGGGCAAAGCCCACCATATCGCGCTCATCCCATTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTgg                                                                                                                                                                                   <  2:86879/250‑1 (MQ=255)
                                                   cggccaaaggaccACCGCCGCGATTTCATGACGATGCGCCGCCATCAGGCGGGCAAAGCCCACCATATCGCGCTCATCCCATTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTgg                                                                                                                                                                                   <  2:334045/239‑1 (MQ=255)
                                                            gATCACCGCCGCGATTTCATGACGATGCGCCGCCATCAGGCGGGCAAAGCCCACCATATCGCGCTCATCCCATTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGGATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGt                                                                                                                                                                          >  2:354391/1‑251 (MQ=255)
                                                                                gccGATGCGCCGACATCAGGCGGGCAAAGCCCACCATATCGCGCTCATCCCATTCGCCATCCATGCGGCTTTGCGGGCCGGGAGCAAACAGTTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTCTCCGGATCGCACACCGACATCGCACAAACGGTATCGCCATGATAACCATTGCGGAAGGTCAGCAAACGCTGTCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATc                                                                                                                                                      >  2:23053/3‑251 (MQ=255)
                                                                                      gcgcCGCCATCAGGCGGGCAAAGCCCACCATATCGCGCTCATCCCATTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCAc                                                                                                                                                 <  1:184246/250‑1 (MQ=255)
                                                                                        gccgccATCAGGCGGGCAAAGCCCACCATATCGCGCTCATCCCATTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTc                                                                                                                                              <  1:362859/251‑1 (MQ=255)
                                                                                        gccgccATCAGGCGGGCAAAGCCCACCATATCGCGCTCATCCCATTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTc                                                                                                                                              <  1:149969/251‑1 (MQ=255)
                                                                                        cccgccATCAGGCGGGCAAAGCCCACCATATCGCGCGCATCCCATTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGAGAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTc                                                                                                                                              <  2:328233/250‑1 (MQ=255)
                                                                                                 aGGCGGGCAAAGCCCACCATATCGCGCTCATCCCATTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTAc                                                                                                                                     <  2:207252/251‑1 (MQ=255)
                                                                                                  ggCGGGCCAAGCCCACCATATCGCGCTCATCCCATTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTAGCCGGATCGCACACCGACATCGCACCAAATGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCACCGCCATTTTCATCGCCACTACCACCGCTACg                                                                                                                                    <  2:281008/251‑1 (MQ=255)
                                                                                                         aaaGCCCACCATATCGCGCTCATCCCATTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTAAGGAACCgg                                                                                                                             >  2:159031/1‑251 (MQ=255)
                                                                                                              ccaccaTATCGCGCTCATCCCATTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTcc                                                                                                                        >  2:379690/1‑251 (MQ=255)
                                                                                                                        gcgcTCATCCCATTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGaaaaa                                                                                                              <  1:186151/251‑1 (MQ=255)
                                                                                                                            tcttcctattagCCATGCATGCGGCTTTGCGGGGCGGGAGAAAACAGGTTTTATGGCAGGTAGCCTTTACACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCAATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGACCCGGAGTCCGCGAGAAAAACGc                                                                                                           >  2:265133/12‑250 (MQ=255)
                                                                                                                              aTCCCATTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACt                                                                                                        <  1:175502/251‑1 (MQ=255)
                                                                                                                                   aTTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGc                                                                                                   >  1:183584/1‑251 (MQ=255)
                                                                                                                                      cGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACGCCAGCgg                                                                                                 <  1:319672/250‑1 (MQ=255)
                                                                                                                                       gCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGtt                                                                                               <  2:203750/251‑1 (MQ=255)
                                                                                                                                            ccatGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTgcggc                                                                                          <  1:143490/251‑1 (MQ=255)
                                                                                                                                                      ttGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGcc                                                                                 <  1:204448/250‑1 (MQ=255)
                                                                                                                                                            ggcgggAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGCCTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCGCTCCAGCGGTTGCGGCGTCATCGCCACCAGtt                                                                          <  1:375072/251‑1 (MQ=255)
                                                                                                                                                               gggAGCCAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGc                                                                       >  1:356875/1‑251 (MQ=255)
                                                                                                                                                                     aaaCAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGcac                                                                  >  1:215514/1‑250 (MQ=255)
                                                                                                                                                                              ttCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTc                                                                                       >  1:419445/1‑220 (MQ=255)
                                                                                                                                                                              ttCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTc                                                                                       <  2:419445/220‑1 (MQ=255)
                                                                                                                                                                                  ggCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAAtggctg                                                     <  2:138727/250‑1 (MQ=255)
                                                                                                                                                                                  ggCAGGTAGCCGTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACTGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGTCGTCATCGCCACCAGTTTGCGGCACAGCTCAAtggct                                                      <  1:419688/249‑1 (MQ=255)
                                                                                                                                                                                        tAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATg                                              >  2:378609/1‑251 (MQ=255)
                                                                                                                                                                                          gCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACAGAAATCGCCCAACAGGTATCGCCATGATAAACATCGCGGAAGGTCAGACAACGCTGGCGCGCTTCGCCTTTGGCATCCCAGTACTGCCACGCCATTTTCATCGCCACTTCCACCGCTACTGAACCGGAGTCCGCGCGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGAACAGCTCACTGGCTGGCGCATggg                                            >  1:409626/1‑251 (MQ=255)
                                                                                                                                                                                                ccACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATa                                       <  1:379690/250‑1 (MQ=255)
                                                                                                                                                                                                 cacaGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATAcc                                     <  1:378609/251‑1 (MQ=255)
                                                                                                                                                                                                  aaaGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCg                                    <  2:95733/249‑1 (MQ=255)
                                                                                                                                                                                                    aGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACTCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGTGTGATACCGcc                                  <  1:13119/251‑1 (MQ=255)
                                                                                                                                                                                                          tgCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACAt                            >  2:354079/1‑251 (MQ=255)
                                                                                                                                                                                                          tgCATTGAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACAt                            >  1:135252/1‑251 (MQ=255)
                                                                                                                                                                                                                gAGTTATCCGGATCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATg                      >  1:296564/1‑251 (MQ=255)
                                                                                                                                                                                                                           aTCGCACACCGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCAt           >  2:392470/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cGACATCGCACCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCg  <  1:228234/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gACATCGCACCAACGGTATCGCCATGCTAACCATTGCGGAAGGTCAGAAAACGCTGGCGAGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACCTCACATGCGACAAGGCATCAATTTNCg  >  2:412890/1‑250 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  
TTTTAACCATTCCGGATGGTACATGCGCATCCCGCCTGCGCCCTGGACAATCGGCTCAATGATCACCGCCGCGATTTCATGACGATGCGCCGCCATCAGGCGGGCAAAGCCCACCATATCGCGCTCATCCCATTCGCCATCCATGCGGCTTTGCGGGGCGGGAGCAAACAGGTTTTCTGGCAGGTAGCCTTTCCACAGACTGTGCATTGAGTTATCCGGATCGCACACCGACATCGCGCCAAAGGTATCGCCATGATAACCATTGCGGAAGGTCAGAAAACGCTGGCGCGCTTCGCCTTTGGCTTGCCAGTACTGCAACGCCATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCG  >  NC_000913/808571‑809049

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