Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,798,063 C→A A294S (GCC→TCC)  pheS ← phenylalanine‑‑tRNA ligase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,798,0630CA96.0% 148.7 / NA 50A294S (GCC→TCC) pheSphenylalanine‑‑tRNA ligase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (23/25);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GGAAGCTGCCGGCAACCGGTTCTACACCGTCAACTTCCAGGCCCGCCATAGTGATTTGATTTGCCAGCGCATCGCTATCAATCGCCGGGTTCACCCATTCGCGTAACCACAGTTCACTGAATTTCATAATCTATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGTAAAGCTGATGTTGGTATCAA  >  NC_000913/1797818‑1798302
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
ggacgCTGCCGCCCACCGGTTCTCCCCCGCCACCTCCCAGGCCCGCCATAGTGATTTGATTTGCCAGCGCATCGCTATCAATCGCCGGGTTCACCCATTCGCGTAACCACAGTTCACTGAATTTCATAATCTATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAAc                                                                                                                                                                                                                                            <  1:454827/247‑1 (MQ=255)
         cGGCAACCGGTCCTACACCGCCAACTTCCAGGCCCGCCATAGTGATTTGATTTGCCAGCGCATCGCTATCAATCGCCGGGTTCACCCATTCGCGTAACCACAGTTCACTGAATTTCATAATCTATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTaaa                                                                                                                                                                                                                                   <  1:499468/251‑1 (MQ=255)
            cAACCGGTTCTACACCGTCAACTTCCAGGCCCGCCATAGTGATTTGATTTGCCAGCGCATCGCTATCAATCGCCGGGTTCACCCATTCGCGTAACCACAGTTCACTGAATTTCATAATCTATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACtt                                                                                                                                                                                                                                >  1:292621/1‑251 (MQ=255)
                              caacctcagggCCCGCCATAGTGATTTGATTTGCCAGCGCATCGCTATCAATCGCCGGGTTCACCCATTCGCGTAACCACAGTTCACTGAATTTCATAATCTATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAAc                                                                                                                                                                                                                <  1:375210/240‑1 (MQ=255)
                               aaaccccgtccccGCCATAGTGATTTGATTTGCCAGCGCATCGCTATCAATCGCCGGGTTCACCCATTCGCGTAACCACAGTTCACTGAATTTCATAATCTATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGtt                                                                                                                                                                                                             <  2:232762/241‑1 (MQ=255)
                                    ccAGGCCCGCCATAGTGATTTGCTTTGCCAGCGCATCGCTATCACTCGCCGGGCTCACCCATTCGCGTAACCACAGTTCACTGAATTTCATAATCGATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGCAGCCGTAACGCAACATAGTCCGACGCTCCATCCCCATCCCGAAGGAGACACCAGAGTAGACTACCGAGTCGATGCCAACGTTACga                                                                                                                                                                                                         >  2:490862/1‑249 (MQ=255)
                                           cGCCATAGTGATTTGATTTGCCAGCGCATCGCTATCAATCGCCGGGTTCACCCATTCGCGTAACCACAGTTCACTGAATTTCATAATCTATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGtt                                                                                                                                                                                                 >  1:251475/1‑251 (MQ=255)
                                           cGCCATAGTGATTTGATTTGCCAGCGCATCGCTATCAATCGCCGGGTTCACCCATTCGCGTAACCACAGTTCACTGAATTTCATAATCTATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGtt                                                                                                                                                                                                 <  1:275953/251‑1 (MQ=255)
                                                                    gcATCGCTATCAATCGCCGGGTTCACCCATTCGCGTAACCACAGTTCACTGAATTTCATAATCTATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGc                                                                                                                                                                        >  1:25440/1‑251 (MQ=255)
                                                                     cactgCTATCAATCGCCGGGTTCACCCATTCGCGTAACCACAGTTCACTGAATTTCATAATCTATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCa                                                                                                                                                                       <  2:448740/247‑1 (MQ=255)
                                                                      gtCGCTATCAAACTCCGGGTTCAACCATTCGCGTAACCACAGTTCACTGAATTTCATAATCTATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCACGAAGGCGAAACCAGAGTAAACTTCAGGGCCGATGCCAACGTTACGCAACACGTTCGAATGCACCATCCCGCAGACCAGCAc                                                                                                                                                                      >  2:443957/2‑251 (MQ=255)
                                                                         gCTATCAATCGCCGGGTTCACCCATTCGCGTAACCACAGTTCACTGAATTTCATAATCTATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTc                                                                                                                                                                   >  1:65771/1‑251 (MQ=255)
                                                                            aTCAATCGCCGGGTTCACCCATTCGCGTAACCACAGTTCACTGAATTTCATAATCTATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAg                                                                                                                                                                <  1:48992/251‑1 (MQ=255)
                                                                               aaTCGCCGGGTTCACCCATTCGCGTAACCACAGTTCACTGAATTTCATAATCTATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCa                                                                                                                                                             <  1:499247/251‑1 (MQ=255)
                                                                                             cccATTCGCGTAACCACAGTTCACTGAATTTCATAATCTATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTAc                                                                                                                                               >  1:14618/1‑251 (MQ=255)
                                                                                             cccATTCGAGTAACCACAGTTCACTGCATTTCATAATCTATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGCAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCACCATCCCGAAGGAGAAACCCGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCCCCATCCCGCAGCCACGCACTTCCAGACATTTACCGATTTTAc                                                                                                                                               >  2:493880/1‑251 (MQ=255)
                                                                                                     cgTAACCACAGTTCACTGAATTTCATAATCTATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACCCATgccg                                                                                                                                       >  2:503510/1‑248 (MQ=255)
                                                                                                         accccaGTTCACGGAATTTCATAATCTATCCCTGCCTTAATTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGGGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGGAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACCCATGACGTCCa                                                                                                                                   <  2:69285/247‑1 (MQ=255)
                                                                                                                    cgtaaTTTCATAATCGATTCCTGCCTTATTTAAACGGTGTGAAGAAACGACGATCGTTTTCGAAGAATGAACGCCGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCAACCCCCGCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGAGGCCCACGTGACGCATCTCGTTCTGATGCACCATCCCGCAACCCAGCACGTACAGCCATTTTCCGTTTTTGCCCATCACCTCCCCTTCTGCCGaa                                                                                                                        <  2:75870/248‑1 (MQ=255)
                                                                                                                     tGAATTTCATAATCTATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCa                                                                                                                                                                       <  1:378225/203‑1 (MQ=255)
                                                                                                                     tGAATTTCATAATCTATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCa                                                                                                                                                                       >  2:378225/1‑203 (MQ=255)
                                                                                                                          ttCATAATCTATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGTTc                                                                                                                   <  2:25440/250‑1 (MQ=255)
                                                                                                                             aTAATCTATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGTTCGGTa                                                                                                               <  1:503510/251‑1 (MQ=255)
                                                                                                                                      tCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTGTTCGAAGAATGAACGCAGGTCGGGGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGTTCGGTAAACGGGaa                                                                                                       <  1:68321/250‑1 (MQ=255)
                                                                                                                                           ccTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTagg                                                                                                  >  1:63968/1‑250 (MQ=255)
                                                                                                                                              tATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTaggaag                                                                                               <  1:373304/250‑1 (MQ=255)
                                                                                                                                               aTTTAAACTGGTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGCCGCCGTAAAGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTCAACTGCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCATCCCAGCACTTCCAGCCATTTACCGTTTTTGCACATGACGTCCACGTCTGCAGAAGGTTCGGTAAACGGGACGTaggaagga                                                                                             <  2:94834/251‑1 (MQ=255)
                                                                                                                                                     aCTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAg                                                                                       <  2:285579/251‑1 (MQ=255)
                                                                                                                                                         gttGACGAAACGCAGATCGGTTTCCAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATTGTCATGCGTCCCATCCCCATCCCTAAGGAGACACCCGAGTAAACTTCCGGGTCTATGCCAACTTTACGCATCACTTTCTGATTCACCAGCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTTCCCATGACGTCCACTTCTGCAGAAGGTTCGGTAAACGGGTAGTAGGAAGGACGGAAGCGaa                                                                                   <  2:429581/250‑1 (MQ=255)
                                                                                                                                                             aGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACAATCCCGCAGCCCAGCACTTCCAGCACTTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGGTCGGTTAACGGGTAGAAGGAAGGACGGAAGGGGATCTg                                                                               >  2:276436/1‑251 (MQ=255)
                                                                                                                                                                           gTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTCGGAAGGACGGAAGCGAATCTGCAAATCTTACTCa                                                                  >  1:27717/1‑250 (MQ=255)
                                                                                                                                                                           gTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCaa                                                                 <  1:355245/251‑1 (MQ=255)
                                                                                                                                                                               tCGAAGAAGGAACGCAGGGCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAaaga                                                              <  1:504968/250‑1 (MQ=255)
                                                                                                                                                                                  aagaaTGCACGCAGGTCGGAGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGACGATGCCAACGTTACGCAACACGTTCAGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGATCGGTAAACGGGAAGGAGGAAGGACGGAAGCCAATATGCAAATCTTCCTCaaaacagtt                                                          >  2:92489/1‑245 (MQ=255)
                                                                                                                                                                                  aagaaTGAACGCAGGGCGGTTACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACCCGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTTCCCATGACGTACACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGtt                                                          <  1:263655/251‑1 (MQ=255)
                                                                                                                                                                                            gCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTAAGCAGGa                                                  >  2:439417/1‑249 (MQ=255)
                                                                                                                                                                                              aggtcggcgataccGTAAAGCACCATAGTCAGACGCTCCATTCCCATCCCGACGGCGACACCAGTGTAAACTACCGGGTCGATTCCAACGTAACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTTCGTCCATTTACCGTTGTTACCCATGACGCCCACTTCTGCAGAAGATTCGGTTACCGGGAAGTAGGAAGGACGGACGCGATACTGCAAATCTTCCTCAAAGAAGTTCCACAGGACGTc                                              >  2:60519/13‑251 (MQ=255)
                                                                                                                                                                                                gTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGACACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCCACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGCAGGTTCGGTAAACGGGAAGTCGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTAAGCAGGAAGTCgg                                            >  2:279631/1‑250 (MQ=255)
                                                                                                                                                                                                        cGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGTTCGGTAAACGGGAAGAAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTg                                    >  1:388643/1‑251 (MQ=255)
                                                                                                                                                                                                         gCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGc                                   >  1:476303/1‑251 (MQ=255)
                                                                                                                                                                                                            gTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCtt                                >  1:485203/1‑251 (MQ=255)
                                                                                                                                                                                                                   aaCATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGtt                         >  2:119223/1‑251 (MQ=255)
                                                                                                                                                                                                                   aaCATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAAATTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTACAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGACGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGtt                         >  1:279770/1‑251 (MQ=255)
                                                                                                                                                                                                                   aaCATAGTCAGACGCTCCATCCCCAGCCCGAAGGAGAAACCAGAGTAAACTCCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGtt                         <  2:251475/251‑1 (MQ=255)
                                                                                                                                                                                                                     cATAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTgg                       <  1:103965/251‑1 (MQ=255)
                                                                                                                                                                                                                      aTAGTCAGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGt                      <  1:465883/251‑1 (MQ=255)
                                                                                                                                                                                                                            aGACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGTAAAGCt                >  1:374052/1‑251 (MQ=255)
                                                                                                                                                                                                                             gACGCTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTCCGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACACATGACGTCCACTTCTGCAGCAGGTTCGGTAAACGGGAAATAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGTAAAGCTg               >  1:286865/1‑251 (MQ=255)
                                                                                                                                                                                                                                 cTCCATCCCCATCCCGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGTTCGGTAACCGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAACGAAGTTACGCAGGAAGTCGTGCAGCGGGCCTTTCCGGTTGGTAAACCTGATGt           >  2:4931/1‑251 (MQ=255)
                                                                                                                                                                                                                                          catcccGAAGGAGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGTAAAGCTGATGTTGGTATcaa  <  1:36080/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
GGAAGCTGCCGGCAACCGGTTCTACACCGTCAACTTCCAGGCCCGCCATAGTGATTTGATTTGCCAGCGCATCGCTATCAATCGCCGGGTTCACCCATTCGCGTAACCACAGTTCACTGAATTTCATAATCTATTCCTGCCTTATTTAAACTGTTTGAGGAAACGCAGATCGTTTTCGAAGAATGAACGCAGGTCGGTGACGCCGTAACGCAACATAGTCAGACGCTCCATCCCCATCCCGAAGGCGAAACCAGAGTAAACTTCCGGGTCGATGCCAACGTTACGCAACACGTTCGGATGCACCATCCCGCAGCCCAGCACTTCCAGCCATTTACCGTTTTTACCCATGACGTCCACTTCTGCAGAAGGTTCGGTAAACGGGAAGTAGGAAGGACGGAAGCGAATCTGCAAATCTTCCTCAAAGAAGTTACGCAGGAAGTCGTGCAGCGTGCCTTTCAGGTTGGTAAAGCTGATGTTGGTATCAA  >  NC_000913/1797818‑1798302

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