Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,272,743 C→T C79C (TGC→TGT chaC → glutathione‑specific gamma‑glutamylcyclotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,272,7430CT97.3% 111.3 / NA 39C79C (TGC→TGTchaCglutathione‑specific gamma‑glutamylcyclotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (21/15);  minor base G (1/0);  total (22/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGCCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGCTTGATCTTGATGATGGACGCACAGTAAACGCCATTGTGTTTATTATGGACCCGCGACATCCAGA  >  NC_000913/1272509‑1272988
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                     
gATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGTCTGCGCCTGAccgcc                                                                                                                                                                                                                                        <  1:86214/250‑1 (MQ=255)
                        tGCCGATAGTAAACCGGCAATTGGTGCCATTGAAGAATAACTATTCTGGTCAGCAGAACAACGGGCGGCTTCGCTGGCAGCGAAGCTGGCTTGTCGACCTGATGAGGGACCGGCGTGGAACTTCGGCCATGGATCGTTGATGCAGAATCCGGCACTGGAGTTTACCGCATCGTGCACAGGTACAATGGTTGGAAGGCATCGCGCATTTTGGCGCCGCCTGACCGCCGGGCGCGGGACGGCGAACCAGCCgg                                                                                                                                                                                                               >  2:401676/1‑251 (MQ=255)
                         gCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCggg                                                                                                                                                                                                              >  1:38344/1‑251 (MQ=255)
                          ccGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGCCCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGa                                                                                                                                                                                                             <  2:124202/251‑1 (MQ=255)
                               tGTAAAAGGGCATTGGGTGCCATTGAAGCATCACTCTTATGGTCAGCAGAACAACGGGCGGATTAGCTGGCGGCGCCGCTGGCTTGTTGACCTGATGAGGGACCGGTGTGGATCTTCGGGTGTGTGGCGTTGATATGGAATCCGGCACTGGAGTTTACCGAATCGTGGACCGGTGCACTGGATGGATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGAGGGAGTGACCACCACACGGGACGGAt                                                                                                                                                                                                        >  2:209824/1‑251 (MQ=255)
                                  aaaaCGGCATTTGGTGCCATTGAAGAATCCCTCTTATGGTCATCAAAACAACGGACGGCTTCGCTGGCGGAAACGCTGGCTTGACGACCTGCTGAGGGACCGGAGTGGATCTTCGGCTTTGGATAGTTGATGTGGAATCATACACGGGAGTTTACTGAATCGTGCACCGGTACAAGGGTTGGATGGCCTCGCGAATCTTGTCTGCGCATGACCACCGATCGCGGCACAGCGCACCACCCGTGACGGAtcat                                                                                                                                                                                                     >  2:118693/1‑248 (MQ=255)
                                                    aTTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCgg                                                                                                                                                                                   >  1:237756/1‑251 (MQ=255)
                                                                       ggTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGGAAGAGGGCGGACGCACCACAGGCGGCGcc                                                                                                                                                                >  2:88360/1‑251 (MQ=255)
                                                                              agAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGCTCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACCCTGGTTGGATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGCAAGAGGGCGGACGCACCACCGGCGTAGCCTATCGAc                                                                                                                                                         >  1:392114/1‑251 (MQ=255)
                                                                                   aacGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCACAGGCGTCGCCTATCGACTGcc                                                                                                                                                     >  1:155402/1‑250 (MQ=255)
                                                                                   aacGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCACAGGCGTCGCCTATCGACTGCCa                                                                                                                                                    <  2:182900/251‑1 (MQ=255)
                                                                                    acGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCACAGGCGTCGCCTATCGACTGCCag                                                                                                                                                   <  1:73933/251‑1 (MQ=255)
                                                                                               cGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTg                                                                                                                                        >  1:221902/1‑251 (MQ=255)
                                                                                                                 cTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCtg                                                                                                                      >  1:368909/1‑251 (MQ=255)
                                                                                                                                   gtcctgtgtGGATCTTCGGTTATGGATCGTTGATGTTGAATCCGGCATTGGAGTTTTCCGAATCGTGCACCGGTTCACTGGTTGGATGTCAGCGCGAATCTTGCCTGCGCCTAACCGGCGGGCGGTGGACTTCGCTCCAGTCGGGACGGATGATTGCACTGAAAGAGGGCGGACGCACCACAGGCGTCGTCTACCGACTGCCAGAAGAGACGCTGGAGCATGAACTAACCCTGTTGTGGAAGCGTGAgatt                                                                                                    <  2:233199/246‑2 (MQ=255)
                                                                                                                                        gtgtGGATCTTCGGGTACGGATCGTTGATGGGGAATCCGGCAATGGAGTTTACCGAATCGTGCACCGGTACAGTGGTTGGATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCAGGCGCGGGATTGCGCACCAGCCGGGACGGATGCTTGCACTGAACGAGGGCGGACGCACCACAGGCGTAGCCTATCGACTGCCAGCAGAAACGCTGGCGCAGGAACCAACCCTGATGTTGAAAGGAGAGAGGATTAc                                                                                               >  2:321681/1‑251 (MQ=255)
                                                                                                                                                  taggCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCTTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGAGGCTTGCACTGAAAGAGGGCGGCCGCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGTGCAGGAACTAACCCTGTTGTGGAAGCGAAGGATGAATACCGGCTGtttg                                                                                     >  2:326614/3‑249 (MQ=255)
                                                                                                                                                  tCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCTCCGGTACACTGGTTGGATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTACACTGAAAGAGGGCGGACGCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGtt                                                                                       <  1:15058/249‑1 (MQ=255)
                                                                                                                                                  tCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAACCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGCCTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCACAGGCGTCGCTTAGCGACGGCCATAAGAGACGCTGGAGCAGGAACTAGCCCTGTTGTGCAAGCGAGATATGATTACCGGCTGtt                                                                                       <  2:215088/249‑1 (MQ=255)
                                                                                                                                                    ggCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCt                                                                                   >  1:402805/1‑251 (MQ=255)
                                                                                                                                                     gCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTg                                                                                  <  1:287509/251‑1 (MQ=255)
                                                                                                                                                                     gtgGAAGCCGGCACTGGAGTTTACCTAATCGTGCGCCGGGACACTGGTTGTATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTTTCAGc                                                                  <  2:373698/251‑1 (MQ=255)
                                                                                                                                                                          aTCCGGCACTGGAGTTTACCGAATCGTGCAGCGGTACACTGGTTGGATTGCATCGCGCATTTTGACTGCGCCTGACCGCCGGGCGCGGGACTGGGCACCAGCCGGGACGGATGCTTGCACTGAAAGATGGGGGACGCACCACAGGCGCCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGcttgat                                                             <  1:38191/251‑1 (MQ=255)
                                                                                                                                                                             cGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGcttgatctt                                                          <  1:382617/251‑1 (MQ=255)
                                                                                                                                                                                  cTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGCTTGATCTtgatg                                                      >  1:60089/1‑250 (MQ=255)
                                                                                                                                                                                   tGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGCTTGATCTTgatgat                                                    >  1:251900/1‑251 (MQ=255)
                                                                                                                                                                                    ggAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGCTTGATCTtgatga                                                     <  1:88360/249‑1 (MQ=255)
                                                                                                                                                                                           aCCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTTCCCAACTGGTGTCAACTTGATATTGATGATTGACGcac                                            >  2:430038/1‑251 (MQ=255)
                                                                                                                                                                                            ccGAATCGTGCATCGGTGCACTGGTTGGATGGCATCGCGCATTTTGTCTGCGCCTGTCCGGCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCACAGGCGTTGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGCTTGATCTTGATGATGGACGcaca                                           <  1:355458/251‑1 (MQ=255)
                                                                                                                                                                                            ccGAATCGGGCACCGGTACACTGGTTGGACTGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGCTTGATCTTGATGATGGACGcaca                                           <  2:368909/251‑1 (MQ=255)
                                                                                                                                                                                                             ccacTGGTTGGATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCACAGGCGGCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGCTTGATCTTGATGATGGACGCACAGTAAACGCCATTGTGtt                          <  1:53985/250‑1 (MQ=255)
                                                                                                                                                                                                                cTGGTTGGATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGCTTGATCTTGATGATGGACGCACAGTAAACGCCATTGTGTttat                       >  1:102796/1‑251 (MQ=255)
                                                                                                                                                                                                                     tGGATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAAACTGGTGTCAGCTTGATCTTGATGATGGACGCACAGTAAACGCCATTGTGTTTATTATgg                  >  1:101090/1‑251 (MQ=255)
                                                                                                                                                                                                                     tGGAGGGCATCGCGCAGTTTGTCTGCGGCTGACCGCCGGGCGCGTGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCAGAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTGGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGCTTGATCTTGAGGATGGACGCACAGTAAACGCCATTGTGTTTATTATg                   <  2:187820/250‑1 (MQ=255)
                                                                                                                                                                                                                      ggATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAgatgat                                                                                                  >  1:76883/1‑170 (MQ=255)
                                                                                                                                                                                                                      ggATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAgatgat                                                                                                  <  2:76883/170‑1 (MQ=255)
                                                                                                                                                                                                                      ggATGGCATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGCTTGATCTTGATGATGGACGCACAGTAAACACCATTGTGTTTATTATgg                  >  1:434947/1‑250 (MQ=255)
                                                                                                                                                                                                                            cATCGCGCATTTTGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGCTTGATCTTGATGATGGACGCACAGTAAACGCCATTGTGTTTATTATGGACCcgcg           >  1:426317/1‑251 (MQ=255)
                                                                                                                                                                                                                                      tttGTCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGCTTGATCTTGATGATGGACGCACAGTAAACGCCATTGTGTTTATTATGGACCCGCGACATCCAGa  >  1:7776/1‑250 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                     
GATAACGCGTGATTTCTTGATGAATGCCGATTGTAAAACGGCATTTGGTGCCATTGAAGAATCACTCTTATGGTCAGCAGAACAACGGGCGGCTTCGCTGGCGGCGACGCTGGCTTGTCGACCTGATGAGGGACCGGTGTGGATCTTCGGCTATGGATCGTTGATGTGGAATCCGGCACTGGAGTTTACCGAATCGTGCACCGGTACACTGGTTGGATGGCATCGCGCATTTTGCCTGCGCCTGACCGCCGGGCGCGGGACTGCGCACCAGCCGGGACGGATGCTTGCACTGAAAGAGGGCGGACGCACCACAGGCGTCGCCTATCGACTGCCAGAAGAGACGCTGGAGCAGGAACTAACCCTGTTGTGGAAGCGAGAGATGATTACCGGCTGTTATCTGCCAACCTGGTGTCAGCTTGATCTTGATGATGGACGCACAGTAAACGCCATTGTGTTTATTATGGACCCGCGACATCCAGA  >  NC_000913/1272509‑1272988

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: