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.5% 132.9 / ‑6.6 40Y110C (TAT→TGT) puuPputrescine:H(+) symporter PuuP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (14/25);  minor base G (1/0);  total (15/25)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CAGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCATCATCA  >  NC_000913/1358613‑1359105
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
cAGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCTCCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATc                                                                                                                                                                                                                                                    <  2:380739/251‑1 (MQ=255)
cAGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATc                                                                                                                                                                                                                                                    <  1:50057/251‑1 (MQ=255)
    cGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCcagc                                                                                                                                                                                                                                                <  1:426336/251‑1 (MQ=255)
        gCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCagcagc                                                                                                                                                                                                                                             <  1:42553/250‑1 (MQ=255)
           aaCGCCTTCTGCTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCATCAGTACAAACAGGGTATTGAAGTTAGCGGCCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAAGCCAGCAGCGATGa                                                                                                                                                                                                                                        <  2:200084/251‑1 (MQ=255)
           aaCGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATg                                                                                                                                                                                                                                         <  1:429007/251‑1 (MQ=255)
            aCGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGCTAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGGTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATg                                                                                                                                                                                                                                         >  1:404830/1‑250 (MQ=255)
                  tctcCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGGAGGTTACCACCCACACCCCCGGCGGCACTTCCGGGGAGGGGGCCGAGAGATAGAATTTCGCCCACAAAGCGGTAATCATCGGTAAAAAAAGACAAACCAACAGCGATGAACAGcc                                                                                                                                                                                                                                  >  2:144449/1‑251 (MQ=255)
                        ttATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGATACCAACCACACCCACGGCGGCACTTACGGGCAGAGGGAGGAGAGATAGATTTTCGCCAAAAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCa                                                                                                                                                                                                                             >  1:151283/1‑250 (MQ=255)
                          aTGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTATCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCTCCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATTATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATaa                                                                                                                                                                                                                          <  1:313351/251‑1 (MQ=255)
                                          cccAGGAATATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATACGCTGCGGTTAAAAGGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGTGCGGAGAGATAGATTTTCGCCAACAAGACGTGAATCATCTGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGTCGTGAGGGtt                                                                                                                                                                                                          <  2:56729/250‑1 (MQ=255)
                                          aCCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCggg                                                                                                                                                                                                            <  1:56613/249‑1 (MQ=255)
                                              ggAAGATCAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATTGCGCCGAAGGTTACCACCCACCCCCACGGCGGCACTTCCGGGCAGAGGGCGGGGAGATAGCTTTTCCCCCACAAGACGTTAATCATCGGCAAAAAGAGACAATCCATCCGCGACTACCAGCCGACCCTAAATCCCACGTGCGGGGTAATc                                                                                                                                                                                                      >  2:58762/1‑251 (MQ=255)
                                                aaGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACATGGTATTTAAGTTAGCGTCCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATGGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGa                                                                                                                                                                                                    <  2:370615/251‑1 (MQ=255)
                                                           caccaTGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGt                                                                                                                                                                                         <  1:382208/251‑1 (MQ=255)
                                                                 gatgGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGTGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGta                                                                                                                                                                                    <  1:246426/250‑1 (MQ=255)
                                                                 gatgGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGtat                                                                                                                                                                                   <  2:369335/251‑1 (MQ=255)
                                                                      agaTTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACGTTTTCGCGTAGGTATAGGcc                                                                                                                                                                              <  2:341858/251‑1 (MQ=255)
                                                                       gaTTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCg                                                                                                                                                                             <  1:215580/251‑1 (MQ=255)
                                                                            gCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCg                                                                                                                                                                        <  2:96864/251‑1 (MQ=255)
                                                                               ccAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTTCCACCCACCCCCACCGCGGCACTTTCGGGAAGAGGGCGGCGAGATAGATTTTCGCCCACAAGACCTTAATCATCGGCAAAAAGAGACAATCCCGCCGCGATGACCACCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTTCGCGTAGGTATATGCCGAAccctcc                                                                                                                                                                     >  1:156622/1‑247 (MQ=255)
                                                                                cAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCt                                                                                                                                                                    <  2:294925/251‑1 (MQ=255)
                                                                                      aaaCAGGGTATTGAAGTTAGCGACCCGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCCTCCAATGGTACCACCCACCCCCACCGCGGTACTTCCAGGAGGAGGGCGGAGAGAGAGATTTTTGCCAACCAGACCTTGATAAGCCGCAGAAAGGGAGAATCCCGCAGCAAAGACCAGACGAACATAAAAGCGACGCGGAGGTTACTAGGCTGTTGCGCGTACGTATAGGACGAAcacggctaccga                                                                                                                                                               >  1:248831/1‑239 (MQ=255)
                                                                                                    aGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGGATAGGCCGAACCGGCCTCCGGAAAATGGCGAACCAGt                                                                                                                                                >  1:403354/1‑251 (MQ=255)
                                                                                                       tAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTTCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAAAAAAACGTTAATCATCGGCCAAAAGAGACAATCCAGCAGCGATGACCAACCGACCATAAATTCGACGCGCGGGGTAATCGACTTTTGCGCGTAGGTATAGGCCGTACCGGCCTCCGGCAAccggacaaaccatttt                                                                                                                                             >  2:224514/1‑236 (MQ=255)
                                                                                                        aGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTg                                                                                                                                             <  2:297597/250‑1 (MQ=255)
                                                                                                         gCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCAGCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGcc                                                                                                                                           <  1:58762/251‑1 (MQ=255)
                                                                                                                     ccGCTCTTCAGATTAGGGGCAGTTAAAATGGCGACGAAGGTTTCCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGACAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGCCCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGAGTAGGGATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGc                                                                                                                               <  2:248831/250‑1 (MQ=255)
                                                                                                                            ccAGATTCGCGGCGGTTAAAAAGGCGACGAAGGTTACCACCCGCACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGACAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAAc                                                                                                                        <  1:480094/250‑1 (MQ=255)
                                                                                                                                                     gacgaAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAAt                                                                                               >  1:208256/1‑251 (MQ=255)
                                                                                                                                                      acgaAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAAt                                                                                               >  1:297225/1‑250 (MQ=255)
                                                                                                                                                           ggTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACg                                                                                         >  1:488243/1‑251 (MQ=255)
                                                                                                                                                                ccccccacccccccGGCGGCCCTTCCGGAAAGAGGGCGGAGAAATAGATTTTCGCCAACAAGACGTTAATCACCGGCAAAAAGAGACAACCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAAAAGCCCACCCGCCAGCGCCAGCAAATAGGACGCCGGa                                                                                    <  1:396844/248‑1 (MQ=255)
                                                                                                                                                                     cacacccacGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGACAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTg                                                                               >  1:107604/1‑251 (MQ=255)
                                                                                                                                                                                cggcACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAAAAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGACGAACCGGCCTCCGGAAGCTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTAACCGTCGCTa                                                                     >  2:125492/1‑250 (MQ=255)
                                                                                                                                                                                            aaaGAGGGAGGAGAGAAAGATTTGCACAAAAAAGACGATAATAAACGGAAAAAAGAGACAATCCAGCAGCGATGACAAGCCGACCAAAAATCCGACGGGCGGGTAAAGCGACTTTGGGGCGAAGGTAAAGGCCGAACCGGCCGCCGGAAACTGGCGAACCAGTTGGCCGGAGCGGATAGGGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACaa                                                        <  2:156622/250‑1 (MQ=255)
                                                                                                                                                                                              agagggggaaaaaatagatTTCCGCCAAAAAGATTGTAATCACCGCAAAAAAGAGACAATCCAGCACCGAGGCCCAGCCGACCATAAACCCGACGTGCGGGTAAATCGATTTTTGCGGGTGGGTATAGGCCGAACCGGCCCCGGGAACCTGGCAACCCAGTTTGCCGTAGCTGATAGCGGAAAACAGCCCCCCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATa                                                      <  2:278347/238‑1 (MQ=255)
                                                                                                                                                                                                     ggAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCaaa                                                 >  1:467704/1‑249 (MQ=255)
                                                                                                                                                                                                                  tcccccAACAAGACGTAAACAATCGGCAAAAAGAGACAATCCAGCAGCGATGCCCAGCCGACCATAAATCCGACGTGCGGGTTAATCGTCTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACCCCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACgg                                   <  2:403354/246‑1 (MQ=255)
                                                                                                                                                                                                                                                  gagaCAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGATAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCatcatcc  >  2:389002/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
CAGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCATCATCA  >  NC_000913/1358613‑1359105

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