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

CAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCC  >  NC_000913/1358623‑1359092
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          
cAACGCCTTCTCCTTTATGCCGTTCCTGAAATACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAATCAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACATGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGAt                                                                                                                                                                                                                              <  1:516197/251‑1 (MQ=255)
                aTGCAGTCCCTGAACTTCCAAGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTTAAGTTAGCGCCCAGGCTGACGCGCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACTTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATTACCAGCCGACCATa                                                                                                                                                                                                              <  2:119649/250‑1 (MQ=255)
                   cAGTCCCTGAACTACCAGGAAGATAAACAACACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGTCCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATc                                                                                                                                                                                                          <  1:450804/251‑1 (MQ=255)
                     gTCCCTGAACTACCAGGAAGGTAAACACCACCATAATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCCCTTCCGGGAAAAGGGGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGGTGACCAGCCGACCCAAAATCCgc                                                                                                                                                                                                       >  1:159872/1‑250 (MQ=255)
                            aaCTACCAGGAAGATAAACCCCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCg                                                                                                                                                                                                 <  1:495924/251‑1 (MQ=255)
                            aaCTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGGCCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCAGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCg                                                                                                                                                                                                 <  1:30537/251‑1 (MQ=255)
                             aCTACCAGGAAGATAAACACCCCCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAAAGGCGACGAAGGTTAAAACCCACACCCACGGCGGCACTTCCGGGAAGGGGGCGGAGAGATAAGTTTTCGCCAACAAGACGTTAAACATCGGCAAAAAGAGAAAATGCAGAAGTGATGAACAGCCGCACATAAATCCGACGTGCgg                                                                                                                                                                                                >  1:414089/1‑251 (MQ=255)
                                 ccAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTa                                                                                                                                                                                            <  1:145570/251‑1 (MQ=255)
                                   atgAAGATAAACACCACCATGATGGAGATTTGCACCTGTACAAACAGGGTATTGAAGTTAGCGACCAGGTGGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAATGGTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGAAGAGATAGATTTTCGCCAACAAGACGTTAAGCATCGGCAAAAAGAGACAAGCCAGCAGCGATGACCAGCCGACCATAAATCCGCCGTGCGGGttgat                                                                                                                                                                                          <  2:152950/249‑4 (MQ=255)
                                    ggAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGACGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGGGACGAAGGTTACCCCCCCCACCCACGGCGGGACTTCCGGGAAGAGGGGGGGGGGAAAGATTTTAGCCAACAAGACGTTAATACTCGGCAACAAGAGACAATCAAGCCGCCATGAAAAGCCGACCaacactccgaacgtcggggtaat                                                                                                                                                                                          >  2:84619/1‑228 (MQ=255)
                                      aaGATAAACACCACCATGCTGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCCCGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAAAGACAATCCAGCAGCGATGACCAGCCGACCCTAAATCCGACGTGCGGGTTAATCg                                                                                                                                                                                        >  1:327962/1‑250 (MQ=255)
                                           aaaCACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACtttt                                                                                                                                                                                  >  1:260022/1‑251 (MQ=255)
                                            aaCACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGAGGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTg                                                                                                                                                                                 <  1:185311/251‑1 (MQ=255)
                                                     atgatgGAGATTTGAACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGGTACCAACCCCACCCACGGCGGCACTTCCGGGGAGAGGGCGGAGAGAAAAATTTTCGCCAACAAAACGTTAATCATCGGGAAAAAGAGACAATCCAGCCGAGATGACCCGACGCCCAAAAATCCGACGGGCGGGGTAATCGACTTTTGAGCGTAggg                                                                                                                                                                        >  2:55115/1‑250 (MQ=255)
                                                          ggAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGAGTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGTTTAATCGACTTTTGCGCGTAGGTATAg                                                                                                                                                                    <  2:115711/250‑1 (MQ=255)
                                                          ggAGATGTGCACCATTACAAACAGGGTATTGAAGTTAGCGACCAGTTTGACGCTCTGAAGATTCGCGGCGGGTAAAATGGCGTCGAAGGTTTCCCCCCACACCCACGGCGGCGTTTCCGGGAATAGGGCGGAGAGATATATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACTATCCAGCAGTTATGACCTGCCGACCATAAATCCGACGTGCGGTTTAATCGGCTTGTGCGCGTAGGTATAg                                                                                                                                                                    <  2:64179/250‑1 (MQ=255)
                                                           gagaTTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGc                                                                                                                                                                  <  1:217805/251‑1 (MQ=255)
                                                            agaTTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATCGCGGCGAAGGGTACCACCCACACCCCCGGGGGGACTTCCGGGGAGAGGGCGGAGAGATAGATTTTCCCCCACAAGCAGTTAATAAACGGCCAAAAGAGACAAATCCACAGCGATGACCAGCCGCCCATAAATCCAACGGGCGGGTTAATCGACtttttagcggagggataagcc                                                                                                                                                                 >  2:33602/1‑234 (MQ=255)
                                                            agaTTTGCACCAGTAAAAACAGGGTATTGAAGTTAGCTACCTGGTTGACGCTCTTCAGATTCGCGGCGGTTAAGATGTAGACTAAGGTTGCCACCCACACCCACGGCGGCTCTTCCGGGAAGAGGGCGGAGAGATAGATTTGTGCCAACAAGACGTGAATCAACGGCAAAAAGAGTCAATCCAGCAGCGATGTCCAGCCGAGGATACGTCCGACGTTCGGGTGAATCGAGTTTTGCGCGTAGGAAAAGGcc                                                                                                                                                                 <  2:509799/251‑1 (MQ=255)
                                                                   cACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATTGCGACGGAGGTTACCACCCCCACCCACGGCGGCACTTCCGGGGAGAGGGCGGGGAGGTAGATTTTCTCCAACAAGACGTTAATCATCGGCAAAAAGAGAAAATCCAACAGCGAAGACCAGCCGAACATAAATCCGCCGTGCGGGTTAATCGAACTTTGCGCGTAGGcataagacgacccgg                                                                                                                                                          >  2:144709/1‑240 (MQ=255)
                                                                          ccAAACAGGGTATTGAAGTTAGCGACCAGGTTTACGCTCTTCAGATTCGCGGCGGTTAAAATTGCGACGTAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAACCCGACGTGCGGGTTAATCGACTATTGCTCTTAGGTATAGGCCGAACCGGCCTCCgg                                                                                                                                                   <  2:100130/250‑1 (MQ=255)
                                                                             aaCAGGGTATTGAAGTTAGAGACCAGGTTGACGCTCGTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCCCACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCAGCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGaaa                                                                                                                                                <  2:327962/251‑1 (MQ=255)
                                                                                       tGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAAAGGCGACGAAGGTTACCACCCACACCCCCGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAAcc                                                                                                                                      >  1:190230/1‑251 (MQ=255)
                                                                                        gAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGGTACCACCCACACCCACGGCGGGCCTTCCGGGAAGGGGGGGGAGAGATAGATTTTCGCCAACCAGACGTTAATCATCGGCAAAAAGAGACAATACAGCAGCGATGACCAGCCGACCATAAATCCGACGGGCGGGGTAATCGACTTTTGCGCGTAGGTATAGGCCCACCCGGACTCCGGAAAATGGCGAAccc                                                                                                                                     >  2:145411/1‑250 (MQ=255)
                                                                                                    cAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTGCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGc                                                                                                                         <  1:294897/251‑1 (MQ=255)
                                                                                                        ttGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGGTACCCCCCCCACCCACGGCGGCACTTCCGGGGAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCCTCGGCAAAAAGAGACAATCCAGCCGCGATGACCAACCGACCCTAAATCCGACGGGCGGGGTAATCGCCTTTTGGGCGTAGGTATAAGCCGAACCGGCCTCCGGAAACAGGCGAACCAGTTTGCCGTAGcgggt                                                                                                                     >  1:507191/1‑247 (MQ=255)
                                                                                                                   cAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCTGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACa                                                                                                          <  1:495309/251‑1 (MQ=255)
                                                                                                                            ggcggATAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGCGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGAGTAGGTATAGGCCGGACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAAACGCACACCcg                                                                                                 >  2:82988/1‑251 (MQ=255)
                                                                                                                             gcggTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCcgc                                                                                                <  1:55115/251‑1 (MQ=255)
                                                                                                                               ggTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGCAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGcc                                                                                                                                                                                                                     >  2:419900/1‑132 (MQ=255)
                                                                                                                               ggTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGcc                                                                                                                                                                                                                     <  1:419900/132‑1 (MQ=255)
                                                                                                                                   aaaaTGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCcgccagcgc                                                                                          >  1:93378/1‑251 (MQ=255)
                                                                                                                                    aaaTGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCcgccagcgcc                                                                                         >  2:175221/1‑251 (MQ=255)
                                                                                                                                             cgaAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGCGATAGATTTTCGCCACCAAGACGTTAATCATCGGCAAAAAGCGACAATCCAGCAGAGATGACCAGCCGACCATAAATACGACGTGCGGGTTAATCGACTTTTGAGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCAAGCCAAt                                                                                  >  1:484345/1‑249 (MQ=255)
                                                                                                                                                 ggTTCCCACCCACCCCCACGGCGCCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGGGCGGGTTAACCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACg                                                                            <  1:328348/251‑1 (MQ=255)
                                                                                                                                                        acccacacccacGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAACTAGGACGCCGGAAc                                                                     >  1:383588/1‑251 (MQ=255)
                                                                                                                                                                             tCCGGGAAGAGGGCGGAGAGATAGATTTGCGCCAACAAGCCGTTAATCATCGGCAAAAAGAGCCACTCCAGCAGCGATGACCAGCCGACCATAAATCCGACTTGCGGGTTAATCGCCTTTTGCGCGTAGGTATAGGCCGACCCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACACCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCgg                                                 <  1:367299/250‑1 (MQ=255)
                                                                                                                                                                                        ggCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCa                                      >  1:421934/1‑250 (MQ=255)
                                                                                                                                                                                          cGGAGAGATAGCTTTTCGCCAACCAGGCGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCAGTAGCTGATAGCGGTAACCAGCACACCCGCCAGCGCCAGCAAATCGGACGCCGGAACGTGACCGTCGCTAATGCAGGAGACAATAcccaa                                    >  2:311715/1‑247 (MQ=255)
                                                                                                                                                                                            gagagaTAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGt                                 >  1:354968/1‑251 (MQ=255)
                                                                                                                                                                                                   agatTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACACTCCAGCAGCGCTGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGCACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTCAACAGCACACCCGCCAGCGCCAGCAAATAGGACGACGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCaact                          >  1:201149/1‑249 (MQ=255)
                                                                                                                                                                                                         tCGCTAACAAGACGTAAATCATCGATAAAAAGAGACAATCCAGCAGCGTGCACCAGCCGACAATAAATCCGACGTGCGGGTGAACCGCCTTTTGCGCGGGGGTATAGGCCAAACCGGCCTCCGGAATCGGGCGAACAATTTGGCCGGAGCTGATAGCGGTAAACACCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAAGACCAAAATTATCAAATACGGt                     <  2:507191/250‑1 (MQ=255)
                                                                                                                                                                                                             caacaaGACGTTACTCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCg                >  2:75838/1‑251 (MQ=255)
                                                                                                                                                                                                                           ccagcGGCAAAAAGAGACAACCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACGGGCGACCCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGcc  <  2:93378/247‑1 (MQ=255)
                                                                                                                                                                                                                           aaaccGGCAAAAACAGACAATCCAGCAGCGCGGCCCAGCCGTCCAAAAATCCGACGGGCGGTTTAAGCGTCTTTTGCGCGTAGGTAGAGTCCGACCCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACCCCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGcc  <  2:507905/247‑1 (MQ=255)
                                                                                                                                                                                                                             atcGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGcc  >  1:442768/1‑249 (MQ=255)
                                                                                                                                                                                                                              tcGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGcc  >  1:238050/1‑248 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          
CAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCC  >  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: