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,1760TG98.5% 212.8 / NA 66E48A (GAG→GCG) glyQglycine‑‑tRNA ligase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (27/38);  total (27/39)
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

GGCACCCAGCGTCGGGTTTTCCCAGTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGCTTAACGCA  >  NC_000913/3724935‑3725420
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
ggCACCCAGCGTCGGGTTTTCCCAGTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGAc                                                                                                                                                                                                                                             <  1:322144/251‑1 (MQ=255)
   aCCCCGCGTCTGGCTTTCCCAGTTATCCTCCACGAAACGGATGTCGTGAATATTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGAGATTGCTTAATGACCACCTGTAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAgg                                                                                                                                                                                                                                          <  2:491299/251‑1 (MQ=255)
    cccAGCGTCGGGTTTTCCCATGTACCTTCCACGAAACGGATGTCGTGAATAGTCGGTTCCATGCCCAGCGCTTTCAGAGAACCTAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTGAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAAAATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGt                                                                                                                                                                                                                                         <  2:537468/251‑1 (MQ=255)
         cGTCGGGTTTTCCCAGTTATCTTCCACGCAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTAAGAGAACCGAGGTACAGATCCTGAATATTGTCCGGCGCTGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGCCCATCGGTCGGGCGACGAGAAGGCTGACCATAAGCAGCCGCCATCGGTTCTGGCCCCCGAGCGCGCCGACAGGTCATTg                                                                                                                                                                                                                                    >  2:688225/1‑251 (MQ=255)
                 tttCCCAGTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGagag                                                                                                                                                                                                                            <  1:389478/250‑1 (MQ=255)
                       aGTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTccc                                                                                                                                                                                                                      <  1:561068/251‑1 (MQ=255)
                                    gAAACGGTTGTCGGGAATAGTCGGGTCCATGACCAGCTCTTTCAGAGAAGAGAGGTAAAGCACCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTACTGTAAACGGTTGGGGTTATTGCCGAAGCGCCAATCGATCGGGCGATGCGAAGGCTGGACATTAGCAGCCGCAATAGGTTCTGGCCCCAGAGCGCGCAGCCAGGTCATTGTGCGATGGGTTGCGGCGCCGACTTccc                                                                                                                                                                                                         >  2:211753/1‑250 (MQ=255)
                                     aaaCGGATGTCGTGAATAGCCGGGTCCATACCCTGCGCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCTATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCAt                                                                                                                                                                                                        <  2:118361/251‑1 (MQ=255)
                                       ccGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCCGAGAACCGGGGTACAGCTCCTGAATATTGTCCGGAGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGAAAACGGTTTGGGTTTTCGCCGTAGCGACCATCGGTCGGGCACCGGGAAAGCTGAACATAAGCCGCCGCCCTCGGTTCTGGCCCCAGCGCGCGCAGACAGGCCATTGGGTGAGAGGTTCCCGCGACGACTTTCAttt                                                                                                                                                                                                      >  2:402599/2‑249 (MQ=255)
                                        cGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAAAATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGGCCCGACTTCCATGTc                                                                                                                                                                                                     >  2:588275/1‑251 (MQ=255)
                                               cGTGAATAGTCGGGTCCACGCCAAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATgg                                                                                                                                                                                               <  1:517179/250‑1 (MQ=255)
                                                      aggcGGGTCCATTCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTTAACATAAGCAGCCGCCATCGGTGCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACa                                                                                                                                                                                        <  1:659206/247‑1 (MQ=255)
                                                       gTCGGGTCCATGCCCTGCTCTTGCAGAGAACCGAGGTACAGCTCCTGAATATTTTCCGGCGATTGCTTAATTACCACCTTGAACTGATAGTAGTTCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGTCTTCCAGGTCCAATGGtttaccaat                                                                                                                                                                                      <  2:412721/251‑8 (MQ=255)
                                                           ggTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTg                                                                                                                                                                                  >  2:25407/1‑251 (MQ=255)
                                                              ccATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCa                                                                                                                                                                                <  1:688225/250‑1 (MQ=255)
                                                                  gCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCt                                                                                                                                                                           <  1:599470/251‑1 (MQ=255)
                                                                           tttCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGccc                                                                                                                                                                  >  1:12292/1‑251 (MQ=255)
                                                                              cAGAGAACCGAGGAACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGt                                                                                                                                                               >  1:215569/1‑251 (MQ=255)
                                                                                      cGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGTTCGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGAAGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGCCTTCCATGTCCAATGGTTAGACAATGGTGCAGCCCTGGCGAGCCCAGAAATCCTg                                                                                                                                                        >  2:525484/1‑250 (MQ=255)
                                                                                        aGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGAAAGGTCATTGGGTGAGAGGTTCCCGCGCCGCCTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGAAATCCTGTa                                                                                                                                                      >  2:91506/1‑250 (MQ=255)
                                                                                                     gAATATTGTCCGGCGATGGCTTAATGACCCCCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAg                                                                                                                                        <  2:12292/251‑1 (MQ=255)
                                                                                                         atTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATAGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGccc                                                                                                                                     <  1:698605/250‑1 (MQ=255)
                                                                                                           tGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTg                                                                                                                                   <  1:640687/250‑1 (MQ=255)
                                                                                                              ccGGCGATGGCTAAATGACCGCCTGGAACTGATAGTAGTGCTGTAAACGGTGGGGGTGTTCGCCGTAGCGACCACCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTGCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCGGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGaa                                                                                                                                <  2:149579/250‑1 (MQ=255)
                                                                                                                ggCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGt                                                                                                                             <  2:215569/251‑1 (MQ=255)
                                                                                                                   gATGGGTTAATGACCACCTGGAACTGATAGTAGTGCGGTAAACGGTTGGGGTTTTTGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCt                                                                                                                          <  2:652025/251‑1 (MQ=255)
                                                                                                                     tGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTACGGATCAAGCCCTGGCAGGTCCTgg                                                                                                                        >  2:146175/1‑251 (MQ=255)
                                                                                                                          tAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTAt                                                                                                                     <  1:25407/249‑1 (MQ=255)
                                                                                                                          tAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTGCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATc                                                                                                                    <  1:424876/250‑1 (MQ=255)
                                                                                                                                ccaccTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACttt                                                                                                             <  1:91013/251‑1 (MQ=255)
                                                                                                                                   ccTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTGTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGc                                                                                                          <  1:231852/251‑1 (MQ=255)
                                                                                                                                      ggAACTGATAGTAGTGCTGGAAACGGTTGGGGTTTTCGCCGCAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCata                                                                                                       <  2:249606/251‑1 (MQ=255)
                                                                                                                                        aaCTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCAtatt                                                                                                     >  1:422295/1‑251 (MQ=255)
                                                                                                                                            gATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCAtattat                                                                                                   <  1:146175/249‑1 (MQ=255)
                                                                                                                                             aTAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTc                                                                                                <  1:588275/251‑1 (MQ=255)
                                                                                                                                                   gtGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCGGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGTATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTg                                                                                          <  1:184337/251‑1 (MQ=255)
                                                                                                                                                       tGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGc                                                                                                                                                                        <  1:273663/169‑1 (MQ=255)
                                                                                                                                                       tGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGc                                                                                                                                                                        >  2:273663/1‑169 (MQ=255)
                                                                                                                                                          aaaCGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACgt                                                                                   <  1:501905/251‑1 (MQ=255)
                                                                                                                                                           aaCGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACgtg                                                                                  >  1:42288/1‑251 (MQ=255)
                                                                                                                                                                         tCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCACACTTTTGCATATTATTTCGTg                                                                                             >  2:366621/1‑226 (MQ=255)
                                                                                                                                                                         tCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTg                                                                                             <  1:366621/226‑1 (MQ=255)
                                                                                                                                                                         tCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAaagga                                                                     >  2:16909/1‑250 (MQ=255)
                                                                                                                                                                              ggaGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGGCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGAt                                                               >  2:515617/3‑251 (MQ=255)
                                                                                                                                                                                aGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATc                                                              <  1:91506/250‑1 (MQ=255)
                                                                                                                                                                                   gACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGt                                                           >  1:559425/1‑250 (MQ=255)
                                                                                                                                                                                   gACCATCGGTCGGGCGACGAGAAGGCTGAACAGAAGCAGCCGCCAGCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAAGGGTGCAGCCCTGGCGAGCCCAGTAATCCGGTAAGTTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAATGAAGGGATCAGta                                                          <  1:654507/251‑1 (MQ=255)
                                                                                                                                                                                      cATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATAc                                                       >  1:194894/1‑251 (MQ=255)
                                                                                                                                                                                      cATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTACCGCGCCGACTTCCAAGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGAATCAAACTTTTGCATATTCTTTCGTGCTGGATAAGTGTGGATTTAAAGGAAGGGATCAGTATAc                                                       >  2:418550/1‑251 (MQ=255)
                                                                                                                                                                                       aTCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGGGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATAcc                                                      <  1:323428/251‑1 (MQ=255)
                                                                                                                                                                                       aTCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATAcc                                                      >  2:250396/1‑251 (MQ=255)
                                                                                                                                                                                          ggtcggGCGACGAGAAGGCTGACCATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGc                                                   >  1:490396/1‑251 (MQ=255)
                                                                                                                                                                                          gatcggGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGAGCGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGCCAGGATCAAGCCCTGGAAGATCCTGGTCTCAAACTTTTGAATCTTATTTCGTGTTGGATACGTGTGGATTTAAAGGAAGGGATCAGTAAACCCgg                                                   >  2:481793/3‑250 (MQ=255)
                                                                                                                                                                                           gttcgcttacgatatggCTGAACATAAGCAGCCGCCATTGGTCCTGGCCCCAGCGTGCGTATACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCTGCCCTGTCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTTGAAGTTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGATGGGATCAGtctccccgtt                                                  <  2:20724/236‑10 (MQ=255)
                                                                                                                                                                                             cggGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCtgg                                                <  1:515617/251‑1 (MQ=255)
                                                                                                                                                                                                 ggacgggaatgcgtaaCATAAGCAGCCGCCATCGGTTCTGGCCCCTGCTCGCGCATACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGTATCAGTATACCCGCtggatgg                                            <  2:656088/237‑1 (MQ=255)
                                                                                                                                                                                                     gagaAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGAAAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGAGTTAAAGGAAGGGATCAGTATACCCGCtggatgga                                           >  2:254961/1‑248 (MQ=255)
                                                                                                                                                                                                     gagaAGGATGAACATAAGAAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGGTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGACGAGCCAAGTAATCCTGTAAGGTCAGGAGCAAGCCCTGGCAGGTCCTGGTATCCGACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACACGCtggatgga                                           >  1:18088/1‑248 (MQ=255)
                                                                                                                                                                                                        aaGGCTGAACATAAGCAGCCGCCATCTGATCTGGCCCAAGCGCGCGCAGACAGGTCATTGAGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGGGCAGCCCTGGCGAGACCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGCAGGTCCTGGTATCAAACTTTGGCATATTATTTCGTGCTGGATACGTGTGGATTTTATAGAAGGGAACAGGATACCCGCTGGATGGAAGatc                                      >  2:364285/1‑249 (MQ=255)
                                                                                                                                                                                                               aaCATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCATGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTAc                              >  1:171133/1‑251 (MQ=255)
                                                                                                                                                                                                                     agcagcCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACg                        <  1:16909/251‑1 (MQ=255)
                                                                                                                                                                                                                           cgccATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGaaaa                  <  1:431549/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAgg            <  1:254961/251‑1 (MQ=255)
                                                                                                                                                                                                                                           cccAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCTAATGGTTGAACAATGGTGCAGCCCTGGCGATCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATc                                                                                                                    <  2:342658/137‑1 (MQ=255)
                                                                                                                                                                                                                                           cccAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCTAATGGTTGAACAATGGTGCAGCCCTGGCGATCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATc                                                                                                                    >  1:342658/1‑137 (MQ=255)
                                                                                                                                                                                                                                           cccAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATTGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGCTTAACGCa  <  1:467919/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
GGCACCCAGCGTCGGGTTTTCCCAGTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGCTTAACGCA  >  NC_000913/3724935‑3725420

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