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.6% 208.7 / NA 69E48A (GAG→GCG) glyQglycine‑‑tRNA ligase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (28/40);  minor base A (1/0);  total (29/40)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CCAGGCACCCAGCGTCGGGTTTTCCCAGTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGCTTAACGCATG  >  NC_000913/3724932‑3725422
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       
ccAGGCACCCAGCGGCGGGTTTTCCCAGTTATCTTCCACGAAACGGATGGCGTGAATAGGCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCa                                                                                                                                                                                                                                                   <  2:279806/251‑1 (MQ=255)
  aGGCACCCAGCGTCGGGTTTTCCCAGTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTTTATACTGTTTGGGTTTTCGCCGTAGCT‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGa                                                                                                                                                                                                                                                 <  1:238038/251‑1 (MQ=255)
                    tttCCCAGTGATCTTCCACGAAACGGAGGTCGTGAATAGCCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAACCGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTgaga                                                                                                                                                                                                                                <  1:221220/250‑1 (MQ=255)
                        ccAGTTATCTTCCACGAAACGGCTGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGCGAGAGGTTc                                                                                                                                                                                                                           >  2:435995/1‑251 (MQ=255)
                        ccAGTTAACTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCGGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTTGGGCGACGAGAAGGCTGAACATAAGGGGCCGCCAGCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTc                                                                                                                                                                                                                           <  2:490619/251‑1 (MQ=255)
                           gTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCGTTCAGAGAACCGAGGTACAGCGCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCcg                                                                                                                                                                                                                        <  2:464540/251‑1 (MQ=255)
                                  tCCACGAAACGGATGTCGGGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATAGTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGGTTTCGCCGGAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCTGTTCTGTCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGcc                                                                                                                                                                                                                    <  1:283754/248‑1 (MQ=255)
                                   ccACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACt                                                                                                                                                                                                                >  1:426898/1‑251 (MQ=255)
                                   ccACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACt                                                                                                                                                                                                                <  1:94986/251‑1 (MQ=255)
                                    cACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACtt                                                                                                                                                                                                               <  1:269449/251‑1 (MQ=255)
                                          aCGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGAAGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATgg                                                                                                                                                                                                         >  1:83351/1‑250 (MQ=255)
                                                   gTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTGCCCGCGCCGACTTCCATGTCCAATGGtt                                                                                                                                                                                                <  2:215513/251‑1 (MQ=255)
                                                     gAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCATGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTAGCCGTAGCG‑ACCATCGGTCGGGAGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCACGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTg                                                                                                                                                                                               >  1:547570/1‑250 (MQ=255)
                                                          gTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACaag                                                                                                                                                                                         <  2:193746/251‑2 (MQ=255)
                                                           tctggTCCATACCAAGCTCTTTCAGAGGCCCGAAGTACAGCTCCTGAATAATGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTCGTGCTGTAAACGGTTGGGGATTTCGCCGTAGCG‑AACATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAACCTCCATAGGTTCTGGCCCCAGCGCGCGCAGCCAGGCAATTGGGTGAGAGGTTCCCGCGCCCACTTCCATGTCCAATGGTTGAACAAtt                                                                                                                                                                                        >  2:382749/4‑250 (MQ=255)
                                                               gTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGc                                                                                                                                                                                    >  1:64169/1‑251 (MQ=255)
                                                                 ccATGCCCAGCACTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGGAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTACCATGTCCAATGGTTGAACAATGGTGCAg                                                                                                                                                                                  <  2:138360/251‑1 (MQ=255)
                                                                  cATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGc                                                                                                                                                                                 <  1:114944/251‑1 (MQ=255)
                                                                   aTGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGcc                                                                                                                                                                                >  2:555558/1‑251 (MQ=255)
                                                                      cccAGCTCGTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTg                                                                                                                                                                             <  1:49212/251‑1 (MQ=255)
                                                                              gttCATAGAACCGAGGGACGTCTCCGTAATATTGTCCTGCTATGGCTTAATGACCACCTGGAACTGATAGTTGTGCTGTAATCGGTTGGGGGTTTCGCCGTAGCT‑ACCATCGGTCAGGCGACGAGAATGCTGAACATTAGCATCCGCTAGCGGTTCGGGCGCCAGCGCGCGCATACCGGTCATTGGGTGAGAGGTTCCCGCGCCGTCGTCCAAGGCCAATGTTCGAACAATGGTGTAGCCCGGTCGAGccc                                                                                                                                                                     <  2:509483/250‑1 (MQ=255)
                                                                               ttCAGAGAACCGAGGTACAGCTCCTGAAGATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGc                                                                                                                                                                       <  1:112499/248‑1 (MQ=255)
                                                                                tCAGAGAACCGAGGTTCAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTTATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCGGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGt                                                                                                                                                                                                                             <  1:326593/193‑1 (MQ=255)
                                                                                tCAGAGAACCGAGGTACAACTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAAAGGTTGTGGTTTTCGCCGTCGCG‑ACCATCGGTCGAGCGACTAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGATAGGt                                                                                                                                                                                                                             >  2:326593/1‑193 (MQ=255)
                                                                                   aagaACCGAGGTGCAGCTCCTGAATATTGTCCGGCGATGGCATAATGACCACCGGGAACTGATAGTATTACTGTAAACGGTTGTGGTTGTCCCCGTAGCG‑ACCAACGGTCGGGTGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCGGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCTCGCCGACTTCCATGGCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTaa                                                                                                                                                                <  2:506437/250‑1 (MQ=255)
                                                                                    agaACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAAt                                                                                                                                                               >  1:255652/1‑251 (MQ=255)
                                                                                           aGGTTCAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTaa                                                                                                                                                        <  1:69500/251‑1 (MQ=255)
                                                                                             gTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑gg                                                                                                                                                     <  1:190446/251‑1 (MQ=255)
                                                                                                cAGCTGCTGAATATTGTCCGGCGATGGCTTAATGACCCCCGGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTACGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCTCAGTAATCCTGTAA‑GGTCa                                                                                                                                                  <  2:101415/251‑1 (MQ=255)
                                                                                                          atagtGTCCGGCGATGGCTTAATGACCCCCTGGAACTGATAGTAGTGCTGGAAACGGGTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GTTCAGGATCAAGcc                                                                                                                                        <  1:450463/247‑1 (MQ=255)
                                                                                                            atTGGCCGGCGAGGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGccc                                                                                                                                       <  2:292263/250‑1 (MQ=255)
                                                                                                             ttGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATg                                                                                                                                                                                                   <  2:427696/190‑1 (MQ=255)
                                                                                                             ttGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATg                                                                                                                                                                                                   >  1:427696/1‑190 (MQ=255)
                                                                                                                   ggCGATGGCTTAATGACCACCTGGAACTGCTAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGt                                                                                                                               >  1:402524/1‑251 (MQ=255)
                                                                                                                    gCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGTc                                                                                                                              >  1:33640/1‑251 (MQ=255)
                                                                                                                         ggCTTAATGACCTCCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCT‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GTTCATGATCAAGCCCTGGAAGGTCCTgg                                                                                                                          <  2:469949/250‑1 (MQ=255)
                                                                                                                              aaTGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCaa                                                                                                                    <  2:100812/251‑1 (MQ=255)
                                                                                                                              aaTGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGGCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCGGGAAGGTCCTGGTATCa                                                                                                                     <  2:113472/250‑1 (MQ=255)
                                                                                                                                   ccaccTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACttt                                                                                                               <  1:555558/251‑1 (MQ=255)
                                                                                                                                      ccTGGAACTGATAGTAGTGATGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACAATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGCCTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCAAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAACCCTGGAAGGGCCTGGTTTCAAACTTTTGc                                                                                                            >  2:87970/1‑251 (MQ=255)
                                                                                                                                         ggAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCTGAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCat                                                                                                          <  2:426898/250‑1 (MQ=255)
                                                                                                                                         ggAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCata                                                                                                         <  2:64169/251‑1 (MQ=255)
                                                                                                                                         ggAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCat                                                                                                          >  1:189254/1‑250 (MQ=255)
                                                                                                                                             cTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCAtattat                                                                                                     >  1:409293/1‑251 (MQ=255)
                                                                                                                                              tGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATattatt                                                                                                    >  1:515883/1‑251 (MQ=255)
                                                                                                                                               gATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATattatt                                                                                                    <  1:477305/250‑1 (MQ=255)
                                                                                                                                               gATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGACAGGATCACGCACTGGAAGGTCCTGGTATCAAACTTTTGCATATTAttt                                                                                                   >  1:541662/1‑251 (MQ=255)
                                                                                                                                                           gcaaaCGGTCGGTGTTTCCGCCCTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACATTTCATTGGGTGAGAGGTTCCCGCGCCTACTGCTATGTCTAATGGTTGAACAATGGTGCAGCCCTGGTGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGGAGGTCCTGGTACCAAACTTTTTCATATTATTTCGGGCTTGATAc                                                                                       <  1:450837/249‑1 (MQ=255)
                                                                                                                                                           gTAAAGGGTGGGGGTTTTCGCC‑TAGGGGACCAGCGGTCGGGCGACGAAAAGGCGGAACATAGGCGGCGGCCACCGGTTCGGGCCCAGGGGCGGGCAGACAGGTCATTGGGTGAGAGGTTCCGGCGCCGACTTCCATGTCCAATGGTGGAACAATGGTGCAGCCCGGGCGAGCCCAGTAA‑CCCGTAAAGTTCAG‑ATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATAc                                                                                       <  1:171724/250‑1 (MQ=255)
                                                                                                                                                             aaaCGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACgg                                                                                     >  1:14853/1‑250 (MQ=255)
                                                                                                                                                                 ggTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTgg                                                                                 >  1:506848/1‑251 (MQ=255)
                                                                                                                                                                  gTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGa                                                                                <  1:374236/251‑1 (MQ=255)
                                                                                                                                                                        ggtgtCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCGGAACAAAAGCAGCCGCAATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTTA‑GGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTGGCATATTATTTCGTGCTGGATTCGTGTGGATTTaaa                                                                          <  2:111087/247‑1 (MQ=255)
                                                                                                                                                                              cccGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGTTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTTGAAGGTCCTGGTATCAAACTTTTGCATATTATTGCGTGCTGGATACGTGTGGATTTAaaggaagg                                                                    <  1:383461/250‑1 (MQ=255)
                                                                                                                                                                                 gTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGAt                                                                 >  2:127972/1‑251 (MQ=255)
                                                                                                                                                                                   aGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATc                                                                <  1:23839/250‑1 (MQ=255)
                                                                                                                                                                                        aCCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGCTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCAAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTCTTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGtat                                                           >  1:290448/1‑251 (MQ=255)
                                                                                                                                                                                             cggtcggGGGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATAGTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCg                                                      <  2:56210/251‑1 (MQ=255)
                                                                                                                                                                                                gcggggggcgaGGAGGGTGAACATAAGCAGCCGGCATCGGGTTTGGCCCCCGCGGGCGCAGACAGGTCATTGGGTGAGAGGGTCCCGCGCCGGCTTCCATGTCCAATGGTTGAACAATGGTGCAGGCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCCAGCCCTGGAAGGTCCTGGTATTAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCtg                                                   <  2:29145/243‑1 (MQ=255)
                                                                                                                                                                                                  gggCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCtgga                                                 <  2:33640/251‑1 (MQ=255)
                                                                                                                                                                                                          cgaAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCACCAGCGCGCGCAGCCAGGTCATTGCGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCACACTTTTGCATATTATTTCGTGCTGGCTACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCtggattga                                             >  2:508025/2‑244 (MQ=255)
                                                                                                                                                                                                                   aaCATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGCGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTAc                                >  2:30385/1‑251 (MQ=255)
                                                                                                                                                                                                                        aaGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAAc                           >  1:117105/1‑251 (MQ=255)
                                                                                                                                                                                                                          gcagcCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGCACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACgg                         >  1:363727/1‑251 (MQ=255)
                                                                                                                                                                                                                            agcCGGCATAGGTTCTGCCCACAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCGGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGGCCTGGTATCAAACTTATGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGAACAGTATCCCCGCtggatggagtagattcagtacgagagggg                        >  2:339947/1‑229 (MQ=255)
                                                                                                                                                                                                                             gccgccATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGa                       >  1:496334/1‑250 (MQ=255)
                                                                                                                                                                                                                             gccgcaATCGGTTCTGGCCCCAGGGCGCGCAGACAGGTCATTGGGGGAGAGGTGCCCGCGACGACTTCCATGTCCAATGGTTGACCAATGGTGCAGCCCTGGCGAGCCCAGTAATCATGTAA‑TGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACGTTAGAATATTAGTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGaa                      <  2:165330/251‑1 (MQ=255)
                                                                                                                                                                                                                                gccATCGGTCCTGGCCACAGCGCGCGCAGCCAGGTCATTGGGTGCGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGCGCAGCCATGGCGAGCCCAGTCATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATCCGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATAAAGTACGAAACGGGaaaa                    >  2:43497/1‑250 (MQ=255)
                                                                                                                                                                                                                                                 cAGCGCGCGAAGACCGGTCATTGGGGGAGAGGTTCCCGCGCCGACTTCCATGTCCCATGGCTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGGCAGGATCAAGCCCTGGAAGGTCCTGGTATCAACCTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAGGGAATCAGTATACCAGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGAGCTTAACGCATg  >  2:528432/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       
CCAGGCACCCAGCGTCGGGTTTTCCCAGTTATCTTCCACGAAACGGATGTCGTGAATAGTCGGGTCCATGCCCAGCTCTTTCAGAGAACCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCG‑ACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAA‑GGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGCTTAACGCATG  >  NC_000913/3724932‑3725422

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: ‑