Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,725,176 T→G E48A (GAG→GCG)  glyQ ← glycine‑‑tRNA ligase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,725,1760TG97.4% 109.3 / NA 39E48A (GAG→GCG) glyQglycine‑‑tRNA ligase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (19/19);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

AGCGTCGGGTTTTCCCAGTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGC  >  NC_000913/3724942‑3725407
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                       
aGCGTCGGGTTTTCCCAGTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTCCAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCAt                                                                                                                                                                                                                         >  2:176093/1‑251 (MQ=255)
             cccAGTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGt                                                                                                                                                                                                             <  1:375291/250‑1 (MQ=255)
              ccAGTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATTTTGTCAGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGAGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCCTCGGTTCTGGCCCCCGCGCGCGCAGGCAGGTCATTGGGTGTGCGGTTc                                                                                                                                                                                                           >  2:2461/1‑251 (MQ=255)
              ccAGTTATCGGCCACGAAACGGATGTCGTGAATAGTCGGGTCCATTCCCAGCTCTGTCAGAGAACCGAGGTACAGCTCCTGAATATAGTCCGGCGAGTGCTTAATGACCGCCTGGAACTGAGAGTAGTGCTGTAAACGGTTGGGGTTTTATCCGTAGCGACCATCGGTCGGGCGACGAGAAGTCTGAACTGAAGCAGCCGCCATCGGTTCTTGCCCCAGTGCACGCATACAGGGCATTGGGTGAGAGGTTc                                                                                                                                                                                                           <  2:19642/251‑1 (MQ=255)
               caggtATCTTCCAAGAGACGGAAGTCGTGAATCGTCGGGTCCAAGCCCCGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATCTTCCGGCGCTGGCTTAATGACCACCTGGAACTGATAGACGAGCAGTAAACGGTTGGGGATTTCGCCGTTGCGAACCTAGGTACGACGCCGAGAAAGCTGAACATAAGGAGCCGCAATAGGTTACGGCACCCGCGCGCGCAGACAGGTCATAGGGTGAGAGGttac                                                                                                                                                                                                          >  2:241812/5‑249 (MQ=255)
                        tCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGGCATTGGGTGAGAGGTTCCCGCGCCGAc                                                                                                                                                                                                 >  1:340167/1‑251 (MQ=255)
                          cACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACtt                                                                                                                                                                                               >  1:254950/1‑251 (MQ=255)
                           aCGAAACGGATGTCGTGAATAGTCGGGTCCCTGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGCCCATCGGTCGGGCGACGAGAAGGCTAAACATAAGCAGCCGCCATCGGTTCTGGCCCCCGCGCGCGCAGCCAGGTCATTGGGGGAGAGGTTCCCGCGCCGACTTc                                                                                                                                                                                              >  2:34482/1‑251 (MQ=255)
                               aaCGGATGTAGTGAGTAGTAGGGACCATGCCCAGCTCTTTCAGATAACCGAGCTACAGCTCCTGAATAATGCACAGCGATGGCTTCATGACCAACTGGAACTGATAGTAGTGCTGTAAGCGGTTGGGGTTTTCGCCGGAGCGACCATAGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCAGCGGTTCTTGCCCCAGCGCGAGCCAACAGGTCATTGTGTGGGAGGTTCCCGCGCCgcctgctac                                                                                                                                                                                           >  2:174801/1‑242 (MQ=255)
                                      gTCGTGAATAGTCGGGTCCATGCCCAGCGCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCCCCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAAt                                                                                                                                                                                    <  2:217621/250‑1 (MQ=255)
                                           gAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGa                                                                                                                                                                              <  1:176093/251‑1 (MQ=255)
                                             aTAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAAc                                                                                                                                                                            <  1:401056/251‑1 (MQ=255)
                                                  ggggTCCATGCCCAGCTCTTCCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGTCGATGGCTTAAGTACCACCTGGAACTGAGAGTTGGGCTGTAAACGGTTGGGGTTTTCGCCGTTGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCTGCCGCCATCGGTTCTGGCCCCAGCGCGAGCAGACAGGTCATGGGGGGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATgg                                                                                                                                                                       <  2:407572/250‑1 (MQ=255)
                                                  cGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATgg                                                                                                                                                                       <  1:317811/251‑1 (MQ=255)
                                                    ggTCCATGCCCAGCTCTTTCAGAGAACCGAGGTCCAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCCCCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTg                                                                                                                                                                     <  1:75657/251‑1 (MQ=255)
                                                       ccATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGc                                                                                                                                                                    >  1:326621/1‑249 (MQ=255)
                                                       ccATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCa                                                                                                                                                                   <  2:340167/250‑1 (MQ=255)
                                                                  tctTTCAGAGAACCGAGGAACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGc                                                                                                                                                       >  1:111554/1‑251 (MQ=255)
                                                                  gcgttCAGAGAACCGAGGTTCAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACGGATAGTAGTGCAGGAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGGGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTGCGGGCCCCAGCTCGCGCAGACAGGGCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGc                                                                                                                                                       <  2:200150/248‑1 (MQ=255)
                                                                   ctTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCTGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAg                                                                                                                                                        <  1:205206/249‑1 (MQ=255)
                                                                          ataaCCGAGGTGCAGCTCCTGAAGATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGGGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCTGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGATCCCAGTAAt                                                                                                                                               <  2:33026/249‑1 (MQ=255)
                                                                           gaACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATc                                                                                                                                              >  1:402603/1‑251 (MQ=255)
                                                                            aaCCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTCATcc                                                                                                                                             >  1:35687/1‑251 (MQ=255)
                                                                                               aaTATTGTCCGGCGATGGCTTCATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACAATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGTCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCAGCGCCGACTTCCATGACCAATGGTTGAACAATGGAGCAGCACTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGc                                                                                                                          >  2:89526/1‑251 (MQ=255)
                                                                                                    tGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTg                                                                                                                      <  1:151937/250‑1 (MQ=255)
                                                                                                          gCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCATACAGGTCATTGGGTGAGAGGTGCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGTCGAGCCCAGGAATGCTGTAAGGTCAGGATCAAGCCCTGGAAGGTc                                                                                                               <  1:249338/251‑1 (MQ=255)
                                                                                                                       gACCACCTGGAACTGATAGTAGTGCTGTAAACGGGTGGGGTTTTCGCCGTAGCGGCCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAAc                                                                                                   <  1:209916/250‑1 (MQ=255)
                                                                                                                        accaccTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTAAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACtt                                                                                                 >  1:440501/1‑251 (MQ=255)
                                                                                                                          caccTGGAACTGCTAGTAGTGCTGTAAACGGTTGGGGGTTTCGCCGCAGCGACCATCGGTCGCGCGACGAGAAGGCTGAACATAAGCAGCCGCCATAGGTTCTGGCGCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCAGACTTCCATGTCCAATGGTTGAACAATGGGGCAGCCATGGCGAGCCCGGTAATCCTGTAAGGCCAGGATCAAGCCCTGGAAGGTCCTCGTCTCAAACtttt                                                                                               >  2:91142/1‑251 (MQ=255)
                                                                                                                          caccTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCACAGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGGGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACtttt                                                                                               <  2:297220/251‑1 (MQ=255)
                                                                                                                                    tGATAGAAGTGCTGTAAACGGATGGGGGTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGATGAACATAAGCAGCCGCCATCGGTTCAGGCCCCAGCGTGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACACTGGTGAAGCCCTGGCGAGCCCAGTAATACTGTAAGGTAAGGATCAAGCCCTGGACGGTACTGGTATCAACCTTTTGCATattatt                                                                                     >  1:437653/1‑251 (MQ=255)
                                                                                                                                            gtGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTg                                                                             <  1:34482/251‑1 (MQ=255)
                                                                                                                                                    aaCGGTTGGGGGTTTCGCCGCAGCGACCATCGGTCGCGCGAAAAGAAGGCTGAACAACAGCCACCGACATCGGCTCTGGCCCCAGCGCGCGCAGCCAGGTCACTGGGTGAGAGGCTCCCGCGCCGACTTCCATGGCAAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGCCAGGATCAAGCCCTGGACGCGCCTGGTATCAAACTTCTGCATACTCTTTCGAGCTGGATAcc                                                                       >  2:354346/1‑248 (MQ=255)
                                                                                                                                                                     ccGTAGCGCTCATCTGTCGGGCGCAGAGAATGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGTTCATTGGGTGAGATGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGTTATCAAACTTTTGCATATTATTTCGTGCTGGATACTTGTGTATTTAAAGGAAggg                                                    <  1:111630/251‑1 (MQ=255)
                                                                                                                                                                      cGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAggg                                                    >  1:402933/1‑250 (MQ=255)
                                                                                                                                                                                  cgttaggtcgacgaGAAGGCTGAACATAAGCAGCCGCCATCGGTTCGTGCCCCAGTGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTTCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATCCCCg                                       <  2:299517/243‑1 (MQ=255)
                                                                                                                                                                                  cggtcggGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATAccc                                        <  1:391757/250‑1 (MQ=255)
                                                                                                                                                                                                                    cgccATCGGTTCTGGCCCCAGCGCTCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATAACCGCTGGATGGAAGATATACAGTACGAAACGGGaaa      >  2:359883/1‑250 (MQ=255)
                                                                                                                                                                                                                         tCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGAATTCCATGTCCAATGGTTGAACAATGGTGCAGCCCCGGCGAGCCCAGTAATCCTGTAAGGCCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGc  >  1:250603/1‑249 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                       
AGCGTCGGGTTTTCCCAGTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGC  >  NC_000913/3724942‑3725407

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