Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,436,697 G→A E60E (GAG→GAA)  trkA → NAD‑binding component of Trk potassium transporters

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,436,6970GA97.8% 127.6 / ‑3.9 45E60E (GAG→GAA) trkANAD‑binding component of Trk potassium transporters
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (22/22);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-01

AAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTG‑ATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAGGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGATAATATTTACCGACTGATTGAGTATCCCGGCGC  >  NC_000913/3436453‑3436942
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
aaGAGGGCGACGTCTTCATTTACGCAAAGCTAATAAAAAAGTG‑ATGAGATAACGGGTCGCGACTGATGAAAATGATCAGTCGGGGTGCCGGCCAGGTTGGCGGCACACTTGCGGAAAACCGGGGTGGCGAGAACAACGATATTACAGTTGTCGATACCAACGGGGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGAGGGAAGCA‑GGt                                                                                                                                                                                                                                                 <  2:214353/251‑1 (MQ=255)
  ccggacaccGGCTTCTTTTACGCTAAGC‑AAACAAAAAGGGGATGAGAAAAGGGGTCGCGATTGTGGAAAATTACCATTCGGGGGGCCGGCCAGGTGGGCGGCCAACTGGCG‑AAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTg                                                                                                                                                                                                                                                <  2:97664/241‑1 (MQ=255)
          cGGCTTCTTTTACGCTAAGCTAATCAAAAAGTG‑ATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCG‑AAAACCTGGTTGGCGAGAACAAAGGTATTACTGTTGTCGGTACCAACGGGGAGCGTCTGCGGACCTTACAGGGTAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCCTCCACGCGTATTGCGGGAAGCA‑gggg                                                                                                                                                                                                                                                >  1:97664/1‑239 (MQ=255)
              ttcttTTACGCTAAGCTAATCAAAAAGTG‑ATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCg                                                                                                                                                                                                                                    >  1:154644/1‑250 (MQ=255)
                       gctaagctaaTCAAAAAGTG‑ATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGt                                                                                                                                                                                                                          <  1:256618/251‑1 (MQ=255)
                       gctaagctaaTCAAAAAGTG‑ATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGGATTTCGGGAAGCA‑GGGGCCGACGACGCCGATATGCTggg                                                                                                                                                                                                                          >  2:256618/1‑250 (MQ=255)
                       gctaagctaaTCAAAAAGAG‑ATGAGATAACGGGTCGAGACTGAGGAAAATTAGCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGCATAAATGTGACTTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTCTTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGt                                                                                                                                                                                                                          <  1:34425/251‑1 (MQ=255)
                               aaTCAAAAAGTG‑ATGAGATAACGGGTCGCGACTGATGAAACTTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGGTGCTGTa                                                                                                                                                                                                                   >  1:120772/1‑250 (MQ=255)
                                   aaaaaGTG‑ATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAg                                                                                                                                                                                                              <  1:85627/251‑1 (MQ=255)
                                      aaGTG‑ATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGGTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGAGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTc                                                                                                                                                                                                           >  1:318248/1‑251 (MQ=255)
                                             ggagaTAACGGGTCGCGACTGATGAAAATTATCATTCTGGGGGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGa                                                                                                                                                                                                     <  2:154644/250‑1 (MQ=255)
                                                                ggatgaAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTtgcctgc                                                                                                                                                                                  <  2:99998/250‑1 (MQ=255)
                                                                 gatgaAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCTCACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTtgcctgc                                                                                                                                                                                  <  1:357766/250‑1 (MQ=255)
                                                                                 cTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACt                                                                                                                                                                 >  1:372423/1‑251 (MQ=255)
                                                                                                   ggcggcACACTGGCGGAAAACCTGGTTGGCGAGCACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCCATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAat                                                                                                                                                 >  2:344787/1‑249 (MQ=255)
                                                                                                         acacTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGt                                                                                                                                                                                                                    <  2:314801/176‑1 (MQ=255)
                                                                                                         acacTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAatcgcatcgc                                                                                                                                         >  1:58268/1‑251 (MQ=255)
                                                                                                         acacTGGCGGAAAAACTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGt                                                                                                                                                                                                                    >  1:314801/1‑176 (MQ=255)
                                                                                                                                cGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACg                                                                                                                  <  2:74653/251‑1 (MQ=255)
                                                                                                                                    aacaacGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTgcg                                                                                                              >  1:381829/1‑251 (MQ=255)
                                                                                                                                               tACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGAt                                                                                                    <  2:97884/250‑1 (MQ=255)
                                                                                                                                                        cGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTAttt                                                                                           <  1:122056/250‑1 (MQ=255)
                                                                                                                                                                        cGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGcc                                                                          <  1:43545/251‑1 (MQ=255)
                                                                                                                                                                            tGCGGCCCTTACAGGATAAAATGGACCTGCGGGCCGTGCAGGGGCATGGCTCTCATCCACGCGTTTTGCGGGAATCA‑GGTGCCGACGACTCCGATATTCTGTTTGCTGTACCCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCGCACCAGACTACGTGTGCGACTCCGATAAGCTATTTCATCCAGATGCTGTGCCGAt                                                                       <  2:215407/250‑1 (MQ=255)
                                                                                                                                                                                  ccTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGatcat                                                                <  1:344787/251‑1 (MQ=255)
                                                                                                                                                                                      aCAGGATAAATTTGGCCTGCGGGGCGTGCAGGGGTATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCGGGTTGCTGTAACCAGTTCAGATGAATCCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAACCGCATCGCTCGCATCCGCTCACCAGACTACGTGCGCGATGCCGAGAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTg                                                             <  2:201106/250‑1 (MQ=255)
                                                                                                                                                                                         ggATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGCAGCA‑GGTGCCGCCGACGCCGATATGCTAGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGCTAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCg                                                         >  1:172191/1‑251 (MQ=255)
                                                                                                                                                                                         ggATAAATTTGACCTGCGGGCCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCg                                                         >  2:346892/1‑251 (MQ=255)
                                                                                                                                                                                                 ttGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGAGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCATACTCACTTTTCAACACCCCTAATCGCATCGATAGTATCCGCTAACCAGACTACGTGAGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATAGCACCag                                                   >  1:119689/1‑249 (MQ=255)
                                                                                                                                                                                                    aCCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCACCCCAGAGc                                                >  2:6839/1‑249 (MQ=255)
                                                                                                                                                                                                      cTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAGGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCACGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGtt                                            >  1:118852/1‑251 (MQ=255)
                                                                                                                                                                                                               gtgaaGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTGCCGGCGACGCCGATATGCTGGCTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATACTGTGCCGATTGATCATCTGCTCGCACCAGAGCAGATGGTTATCg                                    >  1:310704/5‑250 (MQ=255)
                                                                                                                                                                                                                 gCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCATGGT‑CCGACGACGCCGATATGCTGGTTGCTTTAACCAGGTCAGATGAAACCAATATGGTTGACTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGata                                 <  2:218160/251‑1 (MQ=255)
                                                                                                                                                                                                                   aGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCGGATGAAACCAATATGGTTGCCTTCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGTCTACGTGCGCGATGCCGGTAAGCTATTTCATTCAGATGCGGTACAGATTGATCATCTGATCGCACCAGCGCAGTAGGTTATCGgtaa                                >  2:138144/1‑250 (MQ=255)
                                                                                                                                                                                                                       gCATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGATAATAttt                           <  1:240549/251‑1 (MQ=255)
                                                                                                                                                                                                                        cATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCtg                                                                              >  2:347159/1‑199 (MQ=255)
                                                                                                                                                                                                                        cATGGCTCTCATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCtg                                                                              <  1:347159/199‑1 (MQ=255)
                                                                                                                                                                                                                              tctcATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGATAATATTTACCGACt                    >  1:222881/1‑251 (MQ=255)
                                                                                                                                                                                                                              tctcATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGCGCGCGATGCCGATAGGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCCCCAGGGCCGTTGGTTATCGTTAATATTTACCGACt                    >  2:163058/1‑251 (MQ=255)
                                                                                                                                                                                                                               ctcATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTAc                                                                                                                   <  1:139643/155‑1 (MQ=255)
                                                                                                                                                                                                                               ctcATCCACGCGTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTAc                                                                                                                   >  2:139643/1‑155 (MQ=255)
                                                                                                                                                                                                                                       cgcgTATTGCGGGAATCA‑TGTGCCGACGACGCCTATATGCTGGTTTCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCTCGTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGATTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGATAATATTTACCGACTGATTGAGTa           <  2:372476/251‑1 (MQ=255)
                                                                                                                                                                                                                                          gTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTAt                                                                                             <  2:189554/166‑1 (MQ=255)
                                                                                                                                                                                                                                          gTATTGCGGGAAGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTAt                                                                                             >  1:189554/1‑166 (MQ=255)
                                                                                                                                                                                                                                                cGGGAAGCA‑GGTGCAGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATAATCTGATCGCACCAGAGCAGTTGGTTATCGATAATATTTACCGACTGATTGAGTATCCCGgcgc  >  2:56737/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
AAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTG‑ATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAGGCA‑GGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGATAATATTTACCGACTGATTGAGTATCCCGGCGC  >  NC_000913/3436453‑3436942

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