Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,152,494 G→T A87S (GCG→TCG)  metC → cystathionine beta‑lyase/L‑cysteine desulfhydrase/alanine racemase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,152,4940GT97.8% 140.8 / NA 45A87S (GCG→TCG) metCcystathionine beta‑lyase/L‑cysteine desulfhydrase/alanine racemase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (29/15);  total (29/16)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

AGCTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGTT  >  NC_000913/3152249‑3152733
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
aGCTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCa                                                                                                                                                                                                                                             <  1:255943/251‑1 (MQ=255)
 gCTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGgtcggc                                                                                                                                                                                                                                              >  1:117456/1‑249 (MQ=255)
        aCTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGc                                                                                                                                                                                                                                     >  2:172287/1‑251 (MQ=255)
           cAACTGGTGAATGCAGGGCGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCATCTCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACGGGAAGGTGGCGCAGTCTGCGTGCTATTTCCCTGCGAGTCTGCTGCGGTTGCt                                                                                                                                                                                                                                    <  2:27953/249‑1 (MQ=255)
                 gTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCAt                                                                                                                                                                                                                            >  1:353307/1‑251 (MQ=255)
                           gACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTc                                                                                                                                                                                                                               >  1:257737/1‑238 (MQ=255)
                           gACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATAGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTc                                                                                                                                                                                                                               <  2:257737/238‑1 (MQ=255)
                                       aaTACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGc                                                                                                                                                                                                     <  2:284376/251‑1 (MQ=255)
                                               ctcGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCAtgt                                                                                                                                                                                             >  1:281367/1‑251 (MQ=255)
                                                 cGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCAtgt                                                                                                                                                                                             <  2:348004/249‑1 (MQ=255)
                                                    cgcgGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAGTTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTtg                                                                                                                                                                                         >  1:344709/1‑250 (MQ=255)
                                                     gcgGTCAATAGCGTGATTCAGCGCGATTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAACCCCACGCGACACGTACTCGCGCCAATGGAGAGTTGTTCTCTGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTtgat                                                                                                                                                                                       >  2:347814/1‑251 (MQ=255)
                                                     gcgGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTtgat                                                                                                                                                                                       >  1:64485/1‑251 (MQ=255)
                                                                cGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCg                                                                                                                                                                            >  1:211769/1‑251 (MQ=255)
                                                                   gATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCt                                                                                                                                                                         <  2:281367/251‑1 (MQ=255)
                                                                                     gCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGAtt                                                                                                                                                       >  2:366186/1‑251 (MQ=255)
                                                                                     gCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGAtt                                                                                                                                                       >  1:101650/1‑251 (MQ=255)
                                                                                            tttGACAGTGTAGAAGACAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCACCACCGCCTATGAACCGAGTCAGGATTTCTGt                                                                                                                                                  >  2:320391/1‑249 (MQ=255)
                                                                                                   gtgtAGAAGCCAAAAAACACGCGACACGTAATCGCGCAATTGGAGCGTTGTTCTATGGAAGGAGCGGAAAGTTAACCCATTTCTCATTATTACAAGCGCGGTGTGAACTGGAAGGTGGCGCAGGCTGGGTGCTATTTCCCTGCGGGTAGGCAGCGGTTGCTAATTCCTTTCTTGCTTTTATCGAACGGGGAGCTCATGTGTTGATGACCAACACCGACTATGGACCGAGTCAGGATTTCTGTAGCaaaagc                                                                                                                                         >  1:244780/1‑249 (MQ=255)
                                                                                                     gtAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCt                                                                                                                                       >  2:95892/1‑251 (MQ=255)
                                                                                                      tAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTc                                                                                                                                      >  1:120591/1‑251 (MQ=255)
                                                                                                       aGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTc                                                                                                                                      >  1:154141/1‑250 (MQ=255)
                                                                                                         aatacaaaaaatccgccgccaagTAATCGCGCCAGTGGAGAGTTGTTCTATGGACGGCGCGAAACGTAACCCCATTTCTCCTAACAACAAGGGTTGTGTGAACTGGAGGGTGGGCGATGCTGCGTTGGATTTCCTTGGGGGGGGGTTGCGTTTTTTAGTTCCATTCTTGCTTTTATCGAACAGGGGTATCATGGGTTGATGTCAAACGCTGTCGATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAg                                                                                                                                    <  1:212249/229‑1 (MQ=255)
                                                                                                          aatcAAAAAACCCCGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGACCGTTACCCCATTTCTCCTTACAACAAGCGATGTGGGAACTGGAGGGTGGCGCAGGCTGCGTGCTTTTCCCCGGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCAACCATGTTTTGATGACCACCACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGc                                                                                                                                   <  1:282998/248‑1 (MQ=255)
                                                                                                                                tAATCGCGCCAATGGAGAGTTGTTCTATGGATGGCGCGGAACGTTAACCCATTTCTCCTTACAACACGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCACCACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGCcatcat                                                                                                             >  2:164464/1‑251 (MQ=255)
                                                                                                                                    cgcgCCAATGGAGTGTTGTTCTATGGACGGCGCGAAACTTTACCCCTTTTCTCCTTCCACCAAGCGATTTGTGAATTGGAAGGTGGCGCAGGCTGCTTGTTTTTTCCCTGCGGGTCGGCAGCGGTTGTTAATTCCATTCTTGTTTTTATCGAACAGGCCTATCTGGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGt                                                                                                          <  1:356605/250‑1 (MQ=255)
                                                                                                                                       gcCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTg                                                                                                       >  2:362951/1‑250 (MQ=255)
                                                                                                                                                  gtttttCTATGGAGGGCGCGGACCTTTAACCAATTTCCCCTTACAACAAGCGAGGTGTGAACGGGAGGGGGGCGCAGGCTGCGTGTTTTTTCCCTGCGGGTCGGCAGGGGTTGCTAATCCCATTCTTGCTTTTATCGAACAGGGCAATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTg                                                                                           <  2:64485/247‑1 (MQ=255)
                                                                                                                                                          agggACGGCGGGGAACTTTACCCAATTTCTCCTTACAACAGGCGTGGTGTGAACGGGAAGGTGCCGCAGGCTGCGTGCTATTCCCCTGGGGGTCGGCAGCGGTGGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCTTGTGTTGATGACCACCACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGcc                                                                                      <  2:38510/246‑1 (MQ=255)
                                                                                                                                                          aTGGACGGCGCGGAACGTTACCCCATTTCTCCTTACAACAAGCGATGTGGGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGat                                                                                   <  1:172287/251‑1 (MQ=255)
                                                                                                                                                                     ggAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGc                                                                         >  1:356950/1‑250 (MQ=255)
                                                                                                                                                                     ggAACGTAAACCCATTTCTCCTTACAACAAGCGATGTGTGAACGGGAAGGTGGCGCAGGCGGCGTGCTTTTCCCCGGCGGGTCGGCAGCGGTTGCTAATCCCATTCTTGCTTTTAGCGAACAGGGCGATCATGTTTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCa                                                                        <  1:231785/251‑1 (MQ=255)
                                                                                                                                                                          gTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATTATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAGCGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGc                                                                   >  2:253069/1‑251 (MQ=255)
                                                                                                                                                                            tAACCCATTTCTCCTTACAACAAGCGATGTGTGACCTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGc                                                                   >  1:233922/1‑249 (MQ=255)
                                                                                                                                                                              aCCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGcc                                                               >  1:222652/1‑251 (MQ=255)
                                                                                                                                                                               cccATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGCGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGcc                                                               >  1:232237/1‑250 (MQ=255)
                                                                                                                                                                                ccATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCa                                                              >  1:16360/1‑250 (MQ=255)
                                                                                                                                                                                                  aGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTc                                           >  1:361802/1‑251 (MQ=255)
                                                                                                                                                                                                    cGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTg                                         >  1:245450/1‑251 (MQ=255)
                                                                                                                                                                                                              aCTGGAGGGTGGCGCAGGCTGCGTGTTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTTCTTTTATCGAACAGGGCGATCATTTGTTGATGACAAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAg                               <  1:82025/251‑1 (MQ=255)
                                                                                                                                                                                                                 ggAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCCACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCt                            >  2:129867/1‑251 (MQ=255)
                                                                                                                                                                                                                    acgttggcgcgttattctgtcTATTCCCCTGCGGTTCGACAGCTGTTGCAAATTCCATTCTTATTTTACCAGAACCGGGCTATCATGTTGTAATGACCAACTCCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGGTATCGTTAAGCATCTGCAGCCAAACACGATAATCGTGTTGCTAGAATCGCCAGGCTcc                          <  1:141114/231‑1 (MQ=255)
                                                                                                                                                                                                                        ggCGCAGGTTGCGGTCTATTTCCCTGCGGGTCGGCAGGGGTTGCTAATCCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTGGATGCCCAACTCCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCAc                     <  1:52755/251‑1 (MQ=255)
                                                                                                                                                                                                                                gCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGTTCCGCTGATTGGTGCCGATACCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGATGCGCCAGGCTCCATCACTATGGa               >  2:153203/1‑250 (MQ=255)
                                                                                                                                                                                                                                 cTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGa               >  1:320330/1‑249 (MQ=255)
                                                                                                                                                                                                                                            ttCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGtt  >  1:256355/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                
AGCTTGATACTCAACTGGTGAATGCAGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGTT  >  NC_000913/3152249‑3152733

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