Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,358,859 T→C Y110C (TAT→TGT)  puuP ← putrescine:H(+) symporter PuuP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,358,8590TC93.6% 124.3 / ‑5.9 47Y110C (TAT→TGT) puuPputrescine:H(+) symporter PuuP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (21/23);  minor base A (3/0);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

CAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑TAATCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCC  >  NC_000913/1358623‑1359092
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          
cAACGCCTTCTCCTTTATGCCGTTCCTGAAATACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAATCAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACATGACGTTAATCATCGGCAAAAAGAGA‑CAATCCAGCAGCGAt                                                                                                                                                                                                                              <  1:516197/251‑1 (MQ=255)
                aTGCAGTCCCTGAACTTCCAAGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTTAAGTTAGCGCCCAGGCTGACGCGCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACTTTAATCATCGGCAAAAAGAGA‑CAATCCAGCAGCGATT‑ACCAGCCGACCATa                                                                                                                                                                                                              <  2:119649/250‑1 (MQ=255)
                   cAGTCCCTGAACTACCAGGAAGATAAACAACACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGTCCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑CAATCCAGCAGCGATG‑ACCAGCCGACCATAAATc                                                                                                                                                                                                          <  1:450804/251‑1 (MQ=255)
                     gTCCCTGAACTACCAGGAAGGTAAACACCACCATAATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCCCTTCCGGGAAAAGGG‑GGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑CAATCCAGCAGCGGTG‑ACCAGCCGACCCAAAATCCgc                                                                                                                                                                                                       >  1:159872/1‑250 (MQ=255)
                            aaCTACCAGGAAGATAAACCCCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑CAATCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGACGTGCg                                                                                                                                                                                                 <  1:495924/251‑1 (MQ=255)
                            aaCTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGGCCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCAGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑CAATCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGACGTGCg                                                                                                                                                                                                 <  1:30537/251‑1 (MQ=255)
                             aCTACCAGGAAGATAAACACCCCCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAAAGGCGACGAAGGTTAAAACCCACACCCACGGCGGCACTTCCGGGAAGGGGGCGGAGAGATAAGTTTTCGCCAACAAGACGTTAAACATCGGCAAAAAGAGA‑AAATGCAGAAGTGATG‑AACAGCCGCACATAAATCCGACGTGCgg                                                                                                                                                                                                >  1:414089/1‑251 (MQ=255)
                                 ccAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑CAATCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGACGTGCGGGTTa                                                                                                                                                                                            <  1:145570/251‑1 (MQ=255)
                                   atgAAGATAAACACCACCATGATGGAGATTTGCACCTGTACAAACAGGGTATTGAAGTTAGCGACCAGGTGGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAATGGTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGAAGAGATAGATTTTCGCCAACAAGACGTTAAGCATCGGCAAAAAGAGA‑CAAGCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGCCGTGCGGGttgat                                                                                                                                                                                          <  2:152950/249‑4 (MQ=255)
                                    ggAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGACGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGGGACGAAGGTTACCCCCCCCACCCACGGCGGGACTTCCGGGAAGAGGGGGGGGGGAAAGATTTTAGCCAACAAGACGTTAATACTCGGCAACAAGAGA‑CAATCAAGCCGCCATG‑AAAAGCCGACCaacactccgaacgtcggggtaat                                                                                                                                                                                          >  2:84619/1‑228 (MQ=255)
                                      aaGATAAACACCACCATGCTGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCCCGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAAAGA‑CAATCCAGCAGCGATG‑ACCAGCCGACCCTAAATCCGACGTGCGGGTTAATCg                                                                                                                                                                                        >  1:327962/1‑250 (MQ=255)
                                           aaaCACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑CAATCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACtttt                                                                                                                                                                                  >  1:260022/1‑251 (MQ=255)
                                            aaCACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑CAATCCAGCAGCGAGG‑ACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTg                                                                                                                                                                                 <  1:185311/251‑1 (MQ=255)
                                                     atgatgGAGATTTGAACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGGTACCAACCCCACCCACGGCGGCACTTCCGGGGAGAGGGCGGAGAGAAAAATTTTCGCCAACAAAACGTTAATCATCGGGAAAAAGAGA‑CAATCCAGCCGAGATG‑ACCCGACGCCCAAAAATCCGACGGGCGGGGTAATCGACTTTTGAGCGTAggg                                                                                                                                                                        >  2:55115/1‑250 (MQ=255)
                                                          ggAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGAGTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑CAATCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGACGTGCGGTTTAATCGACTTTTGCGCGTAGGTATAg                                                                                                                                                                    <  2:115711/250‑1 (MQ=255)
                                                          ggAGATGTGCACCATTACAAACAGGGTATTGAAGTTAGCGACCAGTTTGACGCTCTGAAGATTCGCGGCGGGTAAAATGGCGTCGAAGGTTTCCCCCCACACCCACGGCGGCGTTTCCGGGAATAGGGCGGAGAGATATATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑CTATCCAGCAGTTATG‑ACCTGCCGACCATAAATCCGACGTGCGGTTTAATCGGCTTGTGCGCGTAGGTATAg                                                                                                                                                                    <  2:64179/250‑1 (MQ=255)
                                                           gagaTTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑CAATCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGc                                                                                                                                                                  <  1:217805/251‑1 (MQ=255)
                                                            agaTTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATCGCGGCGAAGGGTACCACCCACACCCCCGGGGGGACTTCCGGGGAGAGGGCGGAGAGATAGATTTTCCCCCACAAGCAGTTAATAAACGGCCAAAAGAGACAAATCCA‑CAGCGATG‑ACCAGCCGCCCATAAATCCAACGGGCGGGTTAATCGACtttttagcggagggataagcc                                                                                                                                                                 >  2:33602/1‑234 (MQ=255)
                                                            agaTTTGCACCAGTAAAAACAGGGTATTGAAGTTAGCTACCTGGTTGACGCTCTTCAGATTCGCGGCGGTTAAGATGTAGACTAAGGTTGCCACCCACACCCACGGCGGCTCTTCCGGGAAGAGGGCGGAGAGATAGATTTGTGCCAACAAGACGTGAATCAACGGCAAAAAGAGT‑CAATCCAGCAGCGATG‑TCCAGCCGAGGATACGTCCGACGTTCGGGTGAATCGAGTTTTGCGCGTAGGAAAAGGcc                                                                                                                                                                 <  2:509799/251‑1 (MQ=255)
                                                                   cACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATTGCGACGGAGGTTACCACCCCCACCCACGGCGGCACTTCCGGGGAGAGGGCGGGGAGGTAGATTTTCTCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑AAATCCAACAGCGAAG‑ACCAGCCGAACATAAATCCGCCGTGCGGGTTAATCGAACTTTGCGCGTAGGcataagacgacccgg                                                                                                                                                          >  2:144709/1‑240 (MQ=255)
                                                                          ccAAACAGGGTATTGAAGTTAGCGACCAGGTTTACGCTCTTCAGATTCGCGGCGGTTAAAATTGCGACGTAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑CAATCCAGCAGCGATG‑ACCAGCCGACCATAAACCCGACGTGCGGGTTAATCGACTATTGCTCTTAGGTATAGGCCGAACCGGCCTCCgg                                                                                                                                                   <  2:100130/250‑1 (MQ=255)
                                                                             aaCAGGGTATTGAAGTTAGAGACCAGGTTGACGCTCGTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCCCACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCAGCGGCAAAAAGAGA‑CAATCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGaaa                                                                                                                                                <  2:327962/251‑1 (MQ=255)
                                                                                       tGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAAAGGCGACGAAGGTTACCACCCACACCCCCGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑CAATCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAAcc                                                                                                                                      >  1:190230/1‑251 (MQ=255)
                                                                                        gAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGGTACCACCCACACCCACGGCGGGCCTTCCGGGAAGGGGGGGGAGAGATAGATTTTCGCCAACCAGACGTTAATCATCGGCAAAAAGAGA‑CAATACAGCAGCGATG‑ACCAGCCGACCATAAATCCGACGGGCGGGGTAATCGACTTTTGCGCGTAGGTATAGGCCCACCCGGACTCCGGAAAATGGCGAAccc                                                                                                                                     >  2:145411/1‑250 (MQ=255)
                                                                                                    cAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTGCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑CAATCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGc                                                                                                                         <  1:294897/251‑1 (MQ=255)
                                                                                                        ttGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGGTACCCCCCCCACCCACGGCGGCACTTCCGGGGAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCCTCGGCAAAAAGAGA‑CAATCCAGCCGCGATG‑ACCAACCGACCCTAAATCCGACGGGCGGGGTAATCGCCTTTTGGGCGTAGGTATAAGCCGAACCGGCCTCCGGAAACAGGCGAACCAGTTTGCCGTAGcgggt                                                                                                                     >  1:507191/1‑247 (MQ=255)
                                                                                                                   cAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCTGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑CAATCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACa                                                                                                          <  1:495309/251‑1 (MQ=255)
                                                                                                                            ggcggATAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGCGA‑CAATCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGAGTAGGTATAGGCCGGACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAAACGCACACCcg                                                                                                 >  2:82988/1‑251 (MQ=255)
                                                                                                                             gcggTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑CAATCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCcgc                                                                                                <  1:55115/251‑1 (MQ=255)
                                                                                                                               ggTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGCAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑CAATCCAGCAGCGATG‑ACCAGcc                                                                                                                                                                                                                     >  2:419900/1‑132 (MQ=255)
                                                                                                                               ggTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑CAATCCAGCAGCGATG‑ACCAGcc                                                                                                                                                                                                                     <  1:419900/132‑1 (MQ=255)
                                                                                                                                   aaaaTGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑CAATCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCcgccagcgc                                                                                          >  1:93378/1‑251 (MQ=255)
                                                                                                                                    aaaTGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑CAATCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCcgccagcgcc                                                                                         >  2:175221/1‑251 (MQ=255)
                                                                                                                                             cgaAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGCGATAGATTTTCGCCACCAAGACGTTAATCATCGGCAAAAAGCGA‑CAATCCAGCAGAGATG‑ACCAGCCGACCATAAATACGACGTGCGGGTTAATCGACTTTTGAGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCAAGCCAAt                                                                                  >  1:484345/1‑249 (MQ=255)
                                                                                                                                                 ggTTCCCACCCACCCCCACGGCGCCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑CAATCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGACGGGCGGGTTAACCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACg                                                                            <  1:328348/251‑1 (MQ=255)
                                                                                                                                                        acccacacccacGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑CAATCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAACTAGGACGCCGGAAc                                                                     >  1:383588/1‑251 (MQ=255)
                                                                                                                                                                             tCCGGGAAGAGGGCGGAGAGATAGATTTGCGCCAACAAGCCGTTAATCATCGGCAAAAAGAGC‑CACTCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGACTTGCGGGTTAATCGCCTTTTGCGCGTAGGTATAGGCCGACCCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACACCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCgg                                                 <  1:367299/250‑1 (MQ=255)
                                                                                                                                                                                        ggCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑CAATCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCa                                      >  1:421934/1‑250 (MQ=255)
                                                                                                                                                                                          cGGAGAGATAGCTTTTCGCCAACCAGGCGTTAATCATCGGCAAAAAGAGA‑CAATCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCAGTAGCTGATAGCGGTAACCAGCACACCCGCCAGCGCCAGCAAATCGGACGCCGGAACGTGACCGTCGCTAATGCAGGAGACAATAcccaa                                    >  2:311715/1‑247 (MQ=255)
                                                                                                                                                                                            gagagaTAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑CAATCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGt                                 >  1:354968/1‑251 (MQ=255)
                                                                                                                                                                                                   agatTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑CACTCCAGCAGCGCTG‑ACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGCACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTCAACAGCACACCCGCCAGCGCCAGCAAATAGGACGACGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCaact                          >  1:201149/1‑249 (MQ=255)
                                                                                                                                                                                                         tCGCTAACAAGACGTAAATCATCGATAAAAAGAGA‑CAATCCAGCAGCG‑TGCACCAGCCGACAATAAATCCGACGTGCGGGTGAACCGCCTTTTGCGCGGGGGTATAGGCCAAACCGGCCTCCGGAATCGGGCGAACAATTTGGCCGGAGCTGATAGCGGTAAACACCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAAGACCAAAATTATCAAATACGGt                     <  2:507191/250‑1 (MQ=255)
                                                                                                                                                                                                             caacaaGACGTTACTCATCGGCAAAAAGAGA‑CAATCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCg                >  2:75838/1‑251 (MQ=255)
                                                                                                                                                                                                                           ccagcGGCAAAAAGAGA‑CAACCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACGGGCGACCCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGcc  <  2:93378/247‑1 (MQ=255)
                                                                                                                                                                                                                           aaaccGGCAAAAACAGA‑CAATCCAGCAGCGCGG‑CCCAGCCGTCCAAAAATCCGACGGGCGGTTTAAGCGTCTTTTGCGCGTAGGTAGAGTCCGACCCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACCCCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGcc  <  2:507905/247‑1 (MQ=255)
                                                                                                                                                                                                                             atcGGCAAAAAGAGA‑CAATCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGcc  >  1:442768/1‑249 (MQ=255)
                                                                                                                                                                                                                              tcGGCAAAAAGAGA‑CAATCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGcc  >  1:238050/1‑248 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          
CAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGA‑TAATCCAGCAGCGATG‑ACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCC  >  NC_000913/1358623‑1359092

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