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,8590TC95.7% 136.8 / NA 46Y110C (TAT→TGT) puuPputrescine:H(+) symporter PuuP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (19/25);  minor base A (1/1);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

CAGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGG‑AGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGG‑CGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCC  >  NC_000913/1358613‑1359098
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
cAGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGG‑AGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATc                                                                                                                                                                                                                                                <  1:303353/251‑1 (MQ=255)
        gCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGG‑AGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCg                                                                                                                                                                                                                                        <  1:382247/251‑1 (MQ=255)
                            gCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAAAGGCGACGAAGGGTAACACCCACACCCACGGCGGCCCCTCCGGGGAGAGGGAGG‑AGAGATAAATTTTCGCCAACAAGACGTTAATCCTTGGCAAAAGAAGAAAAATCAGCAGCGTTTACCATCCGACCATaaa                                                                                                                                                                                                                     >  1:249360/1‑250 (MQ=255)
                              aGTCCCTGAAATACCAGGAAGATAAACCCCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGGTGACGCTCTTCAGATTCGCGGCGGTTAAAAGGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGG‑AGAGATAGATTTGCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATcc                                                                                                                                                                                                                  <  2:311674/251‑1 (MQ=255)
                                  ccTGAACTACCAGGAAGATAAACACCCCCATGATGGAGATTTGCGCCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGG‑AGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACg                                                                                                                                                                                                             <  2:620093/251‑1 (MQ=255)
                                     gAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGG‑AGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGTg                                                                                                                                                                                                           <  2:487830/250‑1 (MQ=255)
                                                                       gaTTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGG‑AGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCg                                                                                                                                                                       <  1:48546/251‑1 (MQ=255)
                                                                                cAGTACAAACAGGGTAGTGAAGTTAGCGACCAGGTTGACGCGCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACCCCCACGGCGGCACTTCCGGGAAGAGGGCGG‑AGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCt                                                                                                                                                              <  2:90374/251‑1 (MQ=255)
                                                                                cAGTACAAACAGGGGATTGAAGTTAGCTACTAGTTTGACGCTCTTCCGATTCGCGGCGGTTAAAATGGCGACGAAGGGTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGG‑AGAGATAGATTGTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGTCCAGCCGACCATAAATCC‑GACGTGCGGGTTCATCG‑ACGTGTGCGCGTAGGTATAGGCCGAACCGGcgt                                                                                                                                                              <  1:487205/251‑3 (MQ=255)
                                                                                  gTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGG‑AGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTcc                                                                                                                                                            <  1:282792/251‑1 (MQ=255)
                                                                                             gTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATGCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGG‑AGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGG‑Cg                                                                                                                                                <  1:613021/251‑1 (MQ=255)
                                                                                                   aaGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGG‑AGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTcc                                                                                                                                                            <  1:379269/234‑1 (MQ=255)
                                                                                                       gaGCGACCAGGTGTACGCTTTTCAGATTCGCGGAGGTTAAAAGGGCGACGAAGGTTACCACCCACACCCCCGGCGGCACTTCCGGGAAGATGGCGG‑AGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGTGCGGTTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGG‑CGAGCCAGttt                                                                                                                                       <  2:408068/249‑1 (MQ=255)
                                                                                                           gACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCCCCCACACCCCCGGCGGCCCCTCCGGGGAGAGGGCGG‑GGGGATAGATTTTCGCCAACAAGACGTTAATCATCGGCCAAAAGAGACAATCCCGCAGCGATGACCAGCCGACCATAAATCC‑GACGTGCGGGGTAATCG‑ACTTTTGCGCGGAGGTATAGGCCGAACCGGCCTCCGGAAAATGG‑CGGACCAGATTGCCGt                                                                                                                                  >  1:276449/1‑251 (MQ=255)
                                                                                                               aGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGGTACCACCCACACCCACGGCGGCACTTCCGGGGAGAGGGCGG‑AGAGATAGATTTTCGCCAAAAAGACGTTAATTATCGGCAAAAAGAGACAATCCAGCAGCGGTGACCAGCCGACCATAAATCC‑GACGTGCg                                                                                                                                                                                                         >  1:408986/1‑178 (MQ=255)
                                                                                                               aGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAAGGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGG‑AGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGTGCg                                                                                                                                                                                                         <  2:408986/178‑1 (MQ=255)
                                                                                                                    gACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGG‑AGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGG‑CGAACCAGTTTGCCGTAGCTGATAg                                                                                                                         <  1:475040/251‑1 (MQ=255)
                                                                                                                            tCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGG‑AGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGG‑CGAACCAGTTTGCCGTAGCTGATAGCGGTa                                                                                                                    <  1:325353/248‑1 (MQ=255)
                                                                                                                                aTTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGG‑AGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGG‑CGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGca                                                                                                             <  1:320255/251‑1 (MQ=255)
                                                                                                                                   cgcgGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGG‑AGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGG‑CGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGcaca                                                                                                           <  1:453535/250‑1 (MQ=255)
                                                                                                                                          gTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGG‑AGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGAAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGG‑CGAACCAGTTTGCCGTAGCGGATAGCGGTAAACAGCACACCcgccag                                                                                                   <  1:153661/251‑1 (MQ=255)
                                                                                                                                                         aaggatacccccccaccccacGGCGGAACT‑CCGGGAAGGGGGGGGTAGAAATATATTCGCGCAAAAAAGACGTTAACCCCCGGAAAAAAGAGACAATCCGGCAGCGAAGACGAGCCG‑CCAAAAATCCCGGCGGGCGGGTTAACCG‑ACTTTTGCGCGTAGTTAGAGTCATACCCGGCCTCCGGAAACTGG‑CGACCCAGTTTGGAGTAGCAGATTCGGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGa                                                                                    <  2:188452/235‑1 (MQ=255)
                                                                                                                                                            gTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGG‑AGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGCATAGGCCGAACCGGCCTCCGGAAACTGG‑CGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGc                                                                                 >  1:81725/1‑251 (MQ=255)
                                                                                                                                                             ttACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGG‑AGAGATAGATTTTCCCCAACAAGCCGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GAAGTGCGGGGTAATCG‑ACTTTTGCGCGTAGGCATAGGCCGAACCGGCCTCCGGCAACTGG‑CGAACCAGTTTGCCGAAGCTCATAGCGGTACACGGCACACCAGACCGCGCCGGCCAATAGGACGcc                                                                                >  2:181142/1‑251 (MQ=255)
                                                                                                                                                              tcccaccaacACCCACGGCGGCACTTCCGGGACAAGGGCGG‑AGAGATAGATTTTCCCCAACAAGCCGTTAATCATCGGCAAAAGGAGACAATCCAGCAGCGATGACCAGCCGCCCATAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGATTAGGCCGAACCGGCCTCCGGAAACGGC‑CGACCCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACCCCCGCCAGCGCcaacaagtcggacccgg                                                                               >  2:178211/3‑240 (MQ=255)
                                                                                                                                                                 cacccacacccaCGGCGGCACTTCCGGGAAGAGGGCGG‑AGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGG‑CGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGa                                                                             >  1:192937/1‑250 (MQ=255)
                                                                                                                                                                        ccccccggcggCCCTTCGGGGAAGAGGGCGG‑AGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAACCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGGGCGGGTTAATCG‑ACTTTGGCGCGTAGGTATAGGCCGAACCGGCCCCCGGAAACGGG‑CGAACCATTTTGCCGTAGGTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGAcc                                                                     <  1:619698/246‑1 (MQ=255)
                                                                                                                                                                                     cTTCCGGGAAGAGGGCGG‑AGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGCTGACCAGCCGACCATAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGG‑CGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCg                                                        >  1:8335/1‑251 (MQ=255)
                                                                                                                                                                                          gggAAGAGGGCGG‑AGAGATAGATTTTCGCCAACAAGCCGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAACTCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGACTAGGCCGAACCGGCCTCCGGAAACTAG‑CGAACCAGTTTGCCGTAGCTGATAGCGGAAACCAGCACACCCGCCAGACCCAGCCATTAGGACGCCGGACCGTGACCGTCACTAATGCCGGGGAc                                                   >  2:680944/1‑251 (MQ=255)
                                                                                                                                                                                          gggAAGAGGGCGG‑AGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGG‑CGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGAc                                                   >  1:283571/1‑251 (MQ=255)
                                                                                                                                                                                           ggAAGAGGGCGG‑AGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGG‑CGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGgagcc                                                   >  2:239883/1‑248 (MQ=255)
                                                                                                                                                                                           ggAAGAGGGCGG‑AGAGATAGATTTCCGCCAACAAGACTTAAATCACCGGAAAAAAGAGACAACCCAGCAGCGATGACCAGCCGACCAAAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGG‑CGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACCCCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACa                                                  <  1:239883/251‑1 (MQ=255)
                                                                                                                                                                                                agGGCGG‑AGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGCCCATAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGG‑CGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCcgccag                                                                                                   <  1:585497/197‑1 (MQ=255)
                                                                                                                                                                                                agGGCGG‑AGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGG‑CGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCcgccag                                                                                                   >  2:585497/1‑197 (MQ=255)
                                                                                                                                                                                                        agagaTAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGG‑CGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAACTAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAGAAGt                                       >  1:277091/1‑250 (MQ=255)
                                                                                                                                                                                                         gagataatttgtCGCCAAAAAGACGTTAACAATCGGCAAAAAGAGACAATCCAGCGGCGATGACCAGCCGCCCATAAATCC‑GACGTGCGGGTTAAGCG‑ACTTTTGCGCGTAGTTATAGGCCGAACCGGCCTCCGGAAACTGG‑CGAACCAGTTTGCGGTGGCTGATAGCGGTAAACAGCACCCCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTAt                                     <  1:185659/240‑1 (MQ=255)
                                                                                                                                                                                                              agatTTTCGCCAAAAAGACGTTAATCATCGGAAAAAAGAACAAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GTCGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGG‑CGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAAt                                <  1:220868/251‑1 (MQ=255)
                                                                                                                                                                                                               gattgccGCCAACAAGACTTTAATCTTCGGAAAAAAGAGACAATCCAGCAGCGATGACCAGCCGTCCATAATCCC‑GACGTGCGGGTTAACCG‑ACTTTTGCGCGTAGGTATAGGCCAAACCGGCCTCCGGAA‑CTGGGCTAACCAGTTTTCCGTAGCGGAGAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATa                               <  1:618414/245‑1 (MQ=255)
                                                                                                                                                                                                                atTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGAACAGCCGACCATAAATCC‑GACATGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACGGG‑CGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGACCGTGACCGTAGCTAATGCCGGAGACAATACCACACGTATCACATAc                              >  1:647272/1‑251 (MQ=255)
                                                                                                                                                                                                                  tttCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGG‑CGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACgg                            >  1:377755/1‑251 (MQ=255)
                                                                                                                                                                                                                       ccAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGG‑CGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCGAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAAAACCAAAAGCATCAAATACGGTCATc                       >  2:537164/1‑251 (MQ=255)
                                                                                                                                                                                                                            aaGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGG‑CGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGt                  >  1:297997/1‑251 (MQ=255)
                                                                                                                                                                                                                                cgtcaatccccgacAAAAAAAGACAATCCAGCAGCGAAGCCCTGCCGACC‑TAAACCCCGACGTG‑GGGGTAAACGCACTTTTGGGCGTGGGTAAAGGCCGAACGGGCCTCGGGAAACGGG‑GGAACCAGTTGGCCGGAGGGGAAAGCGGTAACCAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTgaga              <  2:61101/238‑1 (MQ=255)
                                                                                                                                                                                                                                 gTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGG‑CGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTgag               >  1:83642/1‑249 (MQ=255)
                                                                                                                                                                                                                                       catcGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGG‑CGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGAATCAAATACGGTCATCGGCGTGAGATAGGcc        >  1:107826/1‑250 (MQ=255)
                                                                                                                                                                                                                                             cAAAAAGCGACAATCCAGCAGCGATGACCAGCCGACCATAAATAC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGG‑CGAACCAGTTTGACGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAATGCGACCGTCGCTAATGCCGGAGACAATACCAACAGTATCAAAAACGGTCATCGGCGTGCGAGAGGcacgacac  >  2:2327/1‑243 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
CAGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGG‑AGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCC‑GACGTGCGGGTTAATCG‑ACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGG‑CGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCC  >  NC_000913/1358613‑1359098

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