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,8590TC97.4% 117.4 / ‑6.7 39Y110C (TAT→TGT) puuPputrescine:H(+) symporter PuuP
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base C (17/20);  minor base A (0/1);  total (17/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 = 9.33e-01

CAGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCA‑TCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAA‑CTGGCGAACCAGTTTGC‑CGTAGCTG‑ATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCC  >  NC_000913/1358613‑1359098
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
cAGACGGTGACAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCA‑TCGGCAAAAAGAGACAATc                                                                                                                                                                                                                                                <  2:252034/251‑1 (MQ=255)
  gACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCA‑TCGGCAAAAAGAGACAATCca                                                                                                                                                                                                                                              <  1:426656/251‑1 (MQ=255)
        gCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCA‑TCGGCAAAAAGAGACAATCCagcagc                                                                                                                                                                                                                                         <  2:401925/250‑1 (MQ=255)
              gCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACAGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCA‑TCGGCAAAAAGAGACAATCCAGCAGCGATGAcc                                                                                                                                                                                                                                  <  2:101664/251‑1 (MQ=255)
                          aTGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCA‑TCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATa                                                                                                                                                                                                                       <  2:163601/250‑1 (MQ=255)
                                          aCCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGGTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCA‑TCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGt                                                                                                                                                                                                       <  2:121975/250‑1 (MQ=255)
                                                    tAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCCGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGGCGAAGGTTACCACCCACACCCACGGCGGGACTTCCGGGGAGAGGGCGGGGGGATAGATTTTCGCCAACAAGACGGTAATCA‑TCGGCAAAAAGAGACAATCCAGCCGCGATGACCAGCCGACCATAAATCCGACGGGCGGGTTAATCGAcctt                                                                                                                                                                                            >  2:62113/1‑248 (MQ=255)
                                                    aaaaCACCACCATGATGGAGATTTGCCCCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACTCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCA‑TCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACttt                                                                                                                                                                                            <  2:173458/250‑1 (MQ=255)
                                                     aaaCACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCA‑TCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACttt                                                                                                                                                                                            <  1:204867/250‑1 (MQ=255)
                                                         accaccATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGGCGAAGGTTACCACCCACACCCACGGCGGTACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGGTAATCA‑ACGGCAAAAAGAGACAATCCCGCAGCGATGACCAGGCGACCATAAATCCGACGTGCGGGGTAATCGACTTTTgcgc                                                                                                                                                                                       >  1:109891/1‑251 (MQ=255)
                                                         accaccATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCA‑TCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTgcgc                                                                                                                                                                                       <  2:275524/251‑1 (MQ=255)
                                                                                       aaCAGTGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCCCTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCA‑TCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCCTAAATCCGACGTGCGGGTTAATCGCCTTTTGCGCGGAGGTATAAGCCGAACCGGCCTCCGGaa                                                                                                                                                          >  2:60718/1‑250 (MQ=255)
                                                                                        aCAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCCCACCCACCGCGGCACTTCCGGGGAGAGGGCGGAGAGATAGATTTTCGCCAAACAGACGTTAATCA‑TCGGCAAAAAGAGACAATCCAGCCGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGGAGGTATAGGCCGAACCGGCCTCCGGAAA‑c                                                                                                                                                       >  1:294047/1‑251 (MQ=255)
                                                                                        aCAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCA‑TCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTGAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAA‑c                                                                                                                                                       <  1:62113/251‑1 (MQ=255)
                                                                                            ggTATTGAAGTTAGCGACCAGGTTGACGCGCTGCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCA‑TCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGGGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAA‑CTGGc                                                                                                                                                   <  2:272912/251‑1 (MQ=255)
                                                                                               aatGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTATAAAGGCGACGAAGGTTACCCCCCCCACCCACGGCGGCCCTTCCGGGGAGAGGGCGGAGAGGTAGATTTTCGCCCAAAAAAAGTTAATCA‑TCGGGCAAAAGAGACAAATCAGAAGCGATGACCAGCCGACCATAAAAACGAAGCGCGCGGTAATCGACATTTGCGCGTAGGTATAAGCCGAAcaggggcccggaaaacggcga                                                                                                                                                  >  2:330741/3‑230 (MQ=255)
                                                                                                 tGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCA‑TCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAA‑CTGGCGAAc                                                                                                                                               <  1:113982/250‑1 (MQ=255)
                                                                                                  gAAGTTAGCGACCAGGTTGACGCGCTTCAGATTCGCTGCGGTTAAAATGTCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGAGAGATTTTCGCCACCAAGACGTTAATCA‑TCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAAGCCGACGGGCGGGTTAATCGACTGTTGCGCTTAGGTATAGGCCTGACCGGCCACCGGAAA‑CTTGCGAAcc                                                                                                                                              <  2:393630/250‑1 (MQ=255)
                                                                                                        aGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCA‑TCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAA‑CTGGCGAACCAGTTTGc                                                                                                                                       <  2:392585/251‑1 (MQ=255)
                                                                                                             ccAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGGGAGATAGATTTTCGCCAACAAGACGTTAATCA‑TAGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGCAACGGCCTCCGGAAA‑CTGGCGAACCAGTTTGC‑CGt                                                                                                                                   >  1:389770/1‑249 (MQ=255)
                                                                                                              cAGGTTGACGCTCTTAAGATTCGCGGCGGTTAAAATGGCGACGAACGCTACCACCCAAACCCACGGCGGCACATCCGGGAAGCGGGCGGAGAGCTAGATTTTCGCCAACAAGACGTTAAACAATGAACAAAAA‑CGACAATACACCAGAGATTACCAGCCGACCATAAATCCCACGAGCGGGGTAATCGAATTTTGCGCGAAGGTATAGGCCGAAGCGGCCTTCCGAAA‑CTGAGGCACCAGTTTGC‑Cgcaac                                                                                                                                >  1:185531/1‑247 (MQ=255)
                                                                                                                                                  ggCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCA‑TCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCgg                                                                                                                                                            >  1:298704/1‑189 (MQ=255)
                                                                                                                                                  ggCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCA‑TCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCgg                                                                                                                                                            <  2:298704/189‑1 (MQ=255)
                                                                                                                                                      acgaAGGTTCCCCCCCACCCCCACGGCGGCACTTCGGGGAAGAGGGCGGAGAGATAGATTTCCGCAAACAAGACGTTAACCA‑TCGAAAAAAAGAGACAATCAAGCAGCGATGACCACCCGACCAAAAATCCGACGTGCGGGTGAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAA‑CTGGCGAACCAGTTTTC‑CATAGCTG‑ATGGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATa                                                                                       <  2:371973/251‑1 (MQ=255)
                                                                                                                                                                                cggcACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCA‑TCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAA‑CTGGCGAACCAGTTTGC‑CGTAGCTG‑ATAGCGGTAAACAGCCCACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTa                                                              >  2:243913/1‑250 (MQ=255)
                                                                                                                                                                                   cACTTCAGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCA‑TCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAA‑CTGGCGAACCAGTTTGC‑CGTAGCTG‑ATAGCGGTAAACAGCAAACCCGCCAGAGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGc                                                          >  2:274511/1‑251 (MQ=255)
                                                                                                                                                                                      ttCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCA‑TCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAA‑CTGGCGAACCAGTTTGC‑CGTAGCTG‑ATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCgg                                                       >  1:19220/1‑251 (MQ=255)
                                                                                                                                                                                        ccGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGCCGTTAATCA‑TCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTACCCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAA‑CTGGCGAACCAGTTTGC‑CGTAGCTG‑ATAGCGGTAAACAGAACACCAGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGcaggcg                                                     >  2:4121/1‑246 (MQ=255)
                                                                                                                                                                                         ggggAAGAGGGCGGAGAAATAGATTTTCGCCAACAAGACGTTAATCA‑CCGGCAAAAAGAGACAATCCAGCAGCGAGGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCCCCGGAAA‑CTGGCGAACCAGTTTGC‑CGTAGCTG‑ATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGgaga                                                    <  2:389770/250‑1 (MQ=255)
                                                                                                                                                                                          gggAAGAGGGCGGAGAGATAGATTTGCGCCAACAAGACGTTAACCA‑TCGACAAAAAGAGACAATCCAGCAGTGATGACCAGCCGACCAAAAATCCGCCGTGCGGGTAAACCGACTTTGGCGCTTAGGTATAGGCCAAACCGGCC‑CCGGAAATCGGGCGAACCAGCTTGC‑CGTAGCAG‑ATAGGGATAAACAGCACACCCGCCAGCGCCAGCAAAGAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGAc                                                   <  2:67046/251‑1 (MQ=255)
                                                                                                                                                                                                               atTTTCGCCAACAAGACGTTAATCA‑TCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGAAGGTATAGGCCGAACCGGCCTCCGGAAA‑CTGGCGAACCAGTTTGC‑CGTAGCTG‑ATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAACTAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATAc                              >  2:131783/1‑251 (MQ=255)
                                                                                                                                                                                                                 ttacgcaaacaaGACGTTAATCA‑TCGGCAAAAAGAGACAATCCAACAGCGATGCCAAGCGGACCATAAATCCGACGTGCGGGTTAAGCGGCTTTTGCGCGTAGGTATAGGCCGACCCGGCCCCCGGAAA‑CGGGCGAACCAGTTTGC‑CGTAGCTG‑AGAGCGGTAAACAGCACGCCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAAGTCCAAAAGTATCAAATACg                             <  1:167256/243‑1 (MQ=255)
                                                                                                                                                                                                                   tgcgcaaaaaagcctttaaccctcGGCAAAAAGACAAAATCCAGCAGCGATGACCAGCCGACTAAAAATCCGACGTGCTGTTTAAGCGACTTTTGCGCGTAGGTATAGGCGGACCCGGCCTCCGGAAC‑CTGGCTAACCAGTT‑GCTCGTAGCGG‑ATAGTGGTAAAAAGCACAACCGCCAGCGTCAGCAAATAGGACGCCTGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACgg                            <  1:265152/228‑1 (MQ=255)
                                                                                                                                                                                                                       caacaaGACGTTAATCA‑TCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAA‑CTGGCGAACCAGTTTGC‑CGTAGCTG‑ATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCg                      >  1:74794/1‑251 (MQ=255)
                                                                                                                                                                                                                            aGACGTTAATCA‑TCTGCAAAAAGAGACAATCCAGCAGCGTTGACCAGCCGACCATAAATCCGCATTTCGGGTTAAGCGCCTTTAGCGCGTGGGTAGAGGCCGACCCGCCCTCCGGAAA‑CGGGCGAACCAGTTTGC‑CGTAG‑TGTATAGCGGTAAACATCACACCCGCCAGCGCCAGCACATAGGACGCCGGAACGTGACAGTCGCTAATGCCGGAGACAAGACCAAAAGTATCAAATACGGTCATCGGCGTg                 <  1:437627/251‑1 (MQ=255)
                                                                                                                                                                                                                               cGTTAATCA‑TCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAA‑CTGGCGAACCAGTTTGC‑CGTAGCTG‑ATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTg                                                                        >  2:77862/1‑193 (MQ=255)
                                                                                                                                                                                                                               cGTTAATCA‑TCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAA‑CGGGCGAACCAGTTTGC‑CGTAGCTG‑ATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTg                                                                        <  1:77862/193‑1 (MQ=255)
                                                                                                                                                                                                                                           ggCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAA‑CTGGCGAACCAGTTTGC‑CGTAGCTG‑ATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGAcc   >  1:130815/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAA‑CTGGCGAACCAGTTTGC‑CGTAGCTG‑ATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGAccc  >  2:98226/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
CAGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCA‑TCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAA‑CTGGCGAACCAGTTTGC‑CGTAGCTG‑ATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCC  >  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: ‑