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

CAGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCC  >  NC_000913/1358613‑1359098
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
cAGACGGTGACAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATc                                                                                                                                                                                                                                                <  2:252034/251‑1 (MQ=255)
  gACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCca                                                                                                                                                                                                                                              <  1:426656/251‑1 (MQ=255)
        gCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCagcagc                                                                                                                                                                                                                                         <  2:401925/250‑1 (MQ=255)
              gCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACAGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGAcc                                                                                                                                                                                                                                  <  2:101664/251‑1 (MQ=255)
                          aTGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATa                                                                                                                                                                                                                       <  2:163601/250‑1 (MQ=255)
                                          aCCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGGTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGt                                                                                                                                                                                                       <  2:121975/250‑1 (MQ=255)
                                                    tAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCCGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGGCGAAGGTTACCACCCACACCCACGGCGGGACTTCCGGGGAGAGGGCGGGGGGATAGATTTTCGCCAACAAGACGGTAATCATCGGCAAAAAGAGACAATCCAGCCGCGATGACCAGCCGACCATAAATCCGACGGGCGGGTTAATCGAcctt                                                                                                                                                                                            >  2:62113/1‑248 (MQ=255)
                                                    aaaaCACCACCATGATGGAGATTTGCCCCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACTCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACttt                                                                                                                                                                                            <  2:173458/250‑1 (MQ=255)
                                                     aaaCACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACttt                                                                                                                                                                                            <  1:204867/250‑1 (MQ=255)
                                                         accaccATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGGCGAAGGTTACCACCCACACCCACGGCGGTACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGGTAATCAACGGCAAAAAGAGACAATCCCGCAGCGATGACCAGGCGACCATAAATCCGACGTGCGGGGTAATCGACTTTTgcgc                                                                                                                                                                                       >  1:109891/1‑251 (MQ=255)
                                                         accaccATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTgcgc                                                                                                                                                                                       <  2:275524/251‑1 (MQ=255)
                                                                                       aaCAGTGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCCCTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCCTAAATCCGACGTGCGGGTTAATCGCCTTTTGCGCGGAGGTATAAGCCGAACCGGCCTCCGGaa                                                                                                                                                          >  2:60718/1‑250 (MQ=255)
                                                                                        aCAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCCCACCCACCGCGGCACTTCCGGGGAGAGGGCGGAGAGATAGATTTTCGCCAAACAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCCGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGGAGGTATAGGCCGAACCGGCCTCCGGAAAc                                                                                                                                                       >  1:294047/1‑251 (MQ=255)
                                                                                        aCAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTGAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAAc                                                                                                                                                       <  1:62113/251‑1 (MQ=255)
                                                                                            ggTATTGAAGTTAGCGACCAGGTTGACGCGCTGCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGGGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGc                                                                                                                                                   <  2:272912/251‑1 (MQ=255)
                                                                                               aatGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTATAAAGGCGACGAAGGTTACCCCCCCCACCCACGGCGGCCCTTCCGGGGAGAGGGCGGAGAGGTAGATTTTCGCCCAAAAAAAGTTAATCATCGGGCAAAAGAGACAAATCAGAAGCGATGACCAGCCGACCATAAAAACGAAGCGCGCGGTAATCGACATTTGCGCGTAGGTATAAGCCGAAcaggggcccggaaaacggcga                                                                                                                                                  >  2:330741/3‑230 (MQ=255)
                                                                                                 tGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAAc                                                                                                                                               <  1:113982/250‑1 (MQ=255)
                                                                                                  gAAGTTAGCGACCAGGTTGACGCGCTTCAGATTCGCTGCGGTTAAAATGTCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGAGAGATTTTCGCCACCAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAAGCCGACGGGCGGGTTAATCGACTGTTGCGCTTAGGTATAGGCCTGACCGGCCACCGGAAACTTGCGAAcc                                                                                                                                              <  2:393630/250‑1 (MQ=255)
                                                                                                        aGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGc                                                                                                                                       <  2:392585/251‑1 (MQ=255)
                                                                                                             ccAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGGGAGATAGATTTTCGCCAACAAGACGTTAATCATAGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGCAACGGCCTCCGGAAACTGGCGAACCAGTTTGCCGt                                                                                                                                   >  1:389770/1‑249 (MQ=255)
                                                                                                              cAGGTTGACGCTCTTAAGATTCGCGGCGGTTAAAATGGCGACGAACGCTACCACCCAAACCCACGGCGGCACATCCGGGAAGCGGGCGGAGAGCTAGATTTTCGCCAACAAGACGTTAAACAATGAACAAAAACGACAATACACCAGAGATTACCAGCCGACCATAAATCCCACGAGCGGGGTAATCGAATTTTGCGCGAAGGTATAGGCCGAAGCGGCCTTCCGAAACTGAGGCACCAGTTTGCCgcaac                                                                                                                                >  1:185531/1‑247 (MQ=255)
                                                                                                                                                  ggCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCgg                                                                                                                                                            >  1:298704/1‑189 (MQ=255)
                                                                                                                                                  ggCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCgg                                                                                                                                                            <  2:298704/189‑1 (MQ=255)
                                                                                                                                                      acgaAGGTTCCCCCCCACCCCCACGGCGGCACTTCGGGGAAGAGGGCGGAGAGATAGATTTCCGCAAACAAGACGTTAACCATCGAAAAAAAGAGACAATCAAGCAGCGATGACCACCCGACCAAAAATCCGACGTGCGGGTGAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTTCCATAGCTGATGGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATa                                                                                       <  2:371973/251‑1 (MQ=255)
                                                                                                                                                                                cggcACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCCCACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTa                                                              >  2:243913/1‑250 (MQ=255)
                                                                                                                                                                                   cACTTCAGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCAAACCCGCCAGAGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGc                                                          >  2:274511/1‑251 (MQ=255)
                                                                                                                                                                                      ttCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCgg                                                       >  1:19220/1‑251 (MQ=255)
                                                                                                                                                                                        ccGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGCCGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTACCCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGAACACCAGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGcaggcg                                                     >  2:4121/1‑246 (MQ=255)
                                                                                                                                                                                         ggggAAGAGGGCGGAGAAATAGATTTTCGCCAACAAGACGTTAATCACCGGCAAAAAGAGACAATCCAGCAGCGAGGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCCCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGgaga                                                    <  2:389770/250‑1 (MQ=255)
                                                                                                                                                                                          gggAAGAGGGCGGAGAGATAGATTTGCGCCAACAAGACGTTAACCATCGACAAAAAGAGACAATCCAGCAGTGATGACCAGCCGACCAAAAATCCGCCGTGCGGGTAAACCGACTTTGGCGCTTAGGTATAGGCCAAACCGGCCCCGGAAATCGGGCGAACCAGCTTGCCGTAGCAGATAGGGATAAACAGCACACCCGCCAGCGCCAGCAAAGAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGAc                                                   <  2:67046/251‑1 (MQ=255)
                                                                                                                                                                                                               atTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGAAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAACTAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATAc                              >  2:131783/1‑251 (MQ=255)
                                                                                                                                                                                                                 ttacgcaaacaaGACGTTAATCATCGGCAAAAAGAGACAATCCAACAGCGATGCCAAGCGGACCATAAATCCGACGTGCGGGTTAAGCGGCTTTTGCGCGTAGGTATAGGCCGACCCGGCCCCCGGAAACGGGCGAACCAGTTTGCCGTAGCTGAGAGCGGTAAACAGCACGCCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAAGTCCAAAAGTATCAAATACg                             <  1:167256/243‑1 (MQ=255)
                                                                                                                                                                                                                   tgcgcaaaaaagcctttaaccctcGGCAAAAAGACAAAATCCAGCAGCGATGACCAGCCGACTAAAAATCCGACGTGCTGTTTAAGCGACTTTTGCGCGTAGGTATAGGCGGACCCGGCCTCCGGAACCTGGCTAACCAGTTGCTCGTAGCGGATAGTGGTAAAAAGCACAACCGCCAGCGTCAGCAAATAGGACGCCTGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACgg                            <  1:265152/228‑1 (MQ=255)
                                                                                                                                                                                                                       caacaaGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCg                      >  1:74794/1‑251 (MQ=255)
                                                                                                                                                                                                                            aGACGTTAATCATCTGCAAAAAGAGACAATCCAGCAGCGTTGACCAGCCGACCATAAATCCGCATTTCGGGTTAAGCGCCTTTAGCGCGTGGGTAGAGGCCGACCCGCCCTCCGGAAACGGGCGAACCAGTTTGCCGTAGTGTATAGCGGTAAACATCACACCCGCCAGCGCCAGCACATAGGACGCCGGAACGTGACAGTCGCTAATGCCGGAGACAAGACCAAAAGTATCAAATACGGTCATCGGCGTg                 <  1:437627/251‑1 (MQ=255)
                                                                                                                                                                                                                               cGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTg                                                                        >  2:77862/1‑193 (MQ=255)
                                                                                                                                                                                                                               cGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACGGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTg                                                                        <  1:77862/193‑1 (MQ=255)
                                                                                                                                                                                                                                           ggCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGAcc   >  1:130815/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGAccc  >  2:98226/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
CAGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCC  >  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: