Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,358,859 T→C Y110C (TAT→TGT)  puuP ← putrescine:H(+) symporter PuuP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,358,8590TC93.2% 144.4 / NA 44Y110C (TAT→TGT) puuPputrescine:H(+) symporter PuuP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (17/24);  minor base A (2/1);  total (19/25)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGAC‑TTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACC  >  NC_000913/1358621‑1359097
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                               
tccAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGATACAATCCagcagc                                                                                                                                                                                                                                      <  1:17816/249‑1 (MQ=255)
gcaaaCGCCGTCTCCTTTATGCAGTCCCGGAACTACCAGGAAGATAAACACCACCATGATGGAGATTGGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGGTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTGACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGAGACAAACCagcagc                                                                                                                                                                                                                                      <  2:283857/247‑1 (MQ=255)
gCCAACGCCTTTTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCACGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCCCCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATGTTCGCC‑AACAAGACGTTAATCAGCGGCAAAAAGAGACAATCCAGCAGCg                                                                                                                                                                                                                                     <  2:117094/251‑1 (MQ=255)
             ccTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGcc                                                                                                                                                                                                                           <  2:511638/248‑1 (MQ=255)
             ccTTTATGCAGGCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGAGTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGAc                                                                                                                                                                                                                        <  2:46510/251‑1 (MQ=255)
                                   ccAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGTTTACCACCCACACCCACGGCGTCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCTACGTGCGGGTTa                                                                                                                                                                                                  <  1:294902/251‑1 (MQ=255)
                                        aaGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGTCGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGa                                                                                                                                                                                             <  1:514790/251‑1 (MQ=255)
                                         aGATAACCACCTCCATGATGGAGATTTGCACCCGTACAAACAGGGTAGTGAAGTTAGCGATCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGTCGAAGGTTCCCCCCCACGCCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGTC‑AACAAGACGGTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGAc                                                                                                                                                                                            <  1:76018/251‑1 (MQ=255)
                                             aaaCACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGAC‑tttt                                                                                                                                                                                       <  1:423180/251‑1 (MQ=255)
                                                caccacCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGAC‑TTTTgcg                                                                                                                                                                                    <  2:104116/251‑1 (MQ=255)
                                                   caccaTGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGAC‑TTTTGCGCGt                                                                                                                                                                                 <  1:420247/251‑1 (MQ=255)
                                                      caTGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGAC‑TTTTGCGCGTAgg                                                                                                                                                                              <  1:474590/251‑1 (MQ=255)
                                                              agaTTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGAGACGAAGGGTACCACCCACACCCACCGCGGCCCTTCCCGGGAGGGGGGGGAGGGATAAATTTTCGCC‑AACAAGACAGTAATCATCGGCCAAAAGAGACAATCCAGCCGCGATGACCAGCCGACCCTAAAACCGACGAGCGGGGGAATCGCC‑TTTTgcaccgagggaacagcc                                                                                                                                                                      >  2:546497/1‑236 (MQ=255)
                                                                           tACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGAC‑TTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCg                                                                                                                                                         >  1:484335/1‑251 (MQ=255)
                                                                               aaCAGGGTATTGAAGTTATCCACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCCCCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCAAC‑TTTTGCGCGTAGGTATAGGCCTAACCGGCCTCCGGaaa                                                                                                                                                     <  2:266372/251‑1 (MQ=255)
                                                                                  aGGGTATTGACGTTAGCGCCCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAAAGGCGACGCAGGGTACCACACACACCCCCGGCGGGCCTTCCGGGAAGAGGGCGGCGGGATAGATTTTCGGC‑CACAAGAAGTTAATCATCGGCAAAAAGAGACAAACCAGCAAGGATGACACGCCCACCATAAATCCGACGTGCCGGGTAAAC‑ACAGTTTTAGCGTAGGTATAAGCCAAACCGGCCCCCGGGAACTg                                                                                                                                                  >  2:213054/1‑251 (MQ=255)
                                                                                  aGGGTATTGAATTTAGCGACCAGGTTGACGCTCTTCAGAGTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCAGCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGAC‑TTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTg                                                                                                                                                  <  1:95828/251‑1 (MQ=255)
                                                                                       aTTGAAGTTAGCGACCAGGTTGACGCTCGTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGAC‑TTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGaa                                                                                                                                             <  2:471533/251‑1 (MQ=255)
                                                                                          gAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGAC‑TTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCa                                                                                                                                          <  1:442872/251‑1 (MQ=255)
                                                                                            cgTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGGGACGAAGGGTACCCCCCCCCACCACGGCGGCACATCCAGGGAGAGGGCGGGGAGATAGATTTTACCC‑AACACCACGGTAATCATCGGCAAAAAACGAACATCCCGCAGCGATGACCAGCCGACCCTAAATCCGACGTGCAGGTTAATCGGC‑CTTTGCGGCTAGGTATAGGCCGAACCCGCCCCCGGAAACTGGCGacacagt                                                                                                                                        >  2:439587/2‑245 (MQ=255)
                                                                                               tAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGAC‑TTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGGGAACCAGttt                                                                                                                                      <  2:381259/250‑1 (MQ=255)
                                                                                                      cAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGGAGGGTAACACCCACCCCCACGGCGGCAATTCCGGGAAGAGGGCGGGGAGATAGATTTTCGCCCAACA‑GACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCCTAAATCCGACGTGCGGGTTAATCGAC‑TTTTGcgcg                                                                                                                                                                                  >  1:375678/1‑199 (MQ=255)
                                                                                                      cAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGAC‑TTTTGcgcg                                                                                                                                                                                  <  2:375678/199‑1 (MQ=255)
                                                                                                         gTTGACGCTCTTCAGATTCGCGGCGGATAAAATGGCGACGAAGGGTACCACCCCCCCCCACGGCGGCACTTCCGGGAAGAGGGCGGGGAGATAGATTTTCGCC‑AAAAAGACGTTAAACATCGGCCAAAAGAGACAATACCGCAGCGGTGACCAGCCGACCATAAATCCCACGGGCGGGGTTATAGTC‑TTTTGCGCGTAGGTATAAGCCGAACCGGCCCCCCGAAACCTGCGAAACCGTTTGCCGtcgccg                                                                                                                            >  2:583534/1‑245 (MQ=255)
                                                                                                                     cagcttCGCGGCTGTTAAAATGGCGACGAAGGTTAACACCCACACACACGGCGGCACTTCCGGGAAGAGGGGGGAGAGATAGGTTTTCGGC‑AACAAGACGTTAATCATCGGCAAAAAGAGAAAAACCAGCCGCGATGAACAGCCGACCCTAAATCCCCAGGACGGGTTAAACGGA‑TTTTTCGCGGAGGTATAGGCCGAAACGGCCTCCCGAAAATAGCGAACCCGTTTTCCGTACGTGTTAGCGGTaaaa                                                                                                                >  2:280025/5‑249 (MQ=255)
                                                                                                                                             gacgaAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGAC‑TTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAAt                                                                                       >  2:119679/1‑251 (MQ=255)
                                                                                                                                               cgaAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGAC‑TTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAAt                                                                                       >  2:183254/1‑249 (MQ=255)
                                                                                                                                                      aaccaccCCCCCCCACGGAGGCACCGCGGGGAAGAGGGAGGAGAGATAGATTGCCGCC‑ATCAAGACGTTAATCATAGGCAAAAAGAGCAAATCCATCAGGGAGGACTGGCCGACCAAAAATCCGACGTGGGGGTTAAGCGCC‑TTTGGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGGTTGCCGTAGCTGATAGCGTTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGcc                                                                               <  2:462056/249‑1 (MQ=255)
                                                                                                                                                       accaccCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGAC‑TTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGCTAGCGGTAACCAGCACACCCGCCAGCGCCAGCAATCAGGCCGCCgg                                                                             >  2:507231/1‑251 (MQ=255)
                                                                                                                                                               cacccacGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGAC‑TTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGAc                                                                     >  1:176746/1‑251 (MQ=255)
                                                                                                                                                                      ggcggcACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGAC‑TTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCt                                                              >  1:90402/1‑251 (MQ=255)
                                                                                                                                                                           cACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGAC‑TTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGc                                                         >  2:140301/1‑251 (MQ=255)
                                                                                                                                                                            aCTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGAC‑TTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCcg                                                                                                      >  2:109380/1‑205 (MQ=255)
                                                                                                                                                                            aCTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGAC‑TTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCcg                                                                                                      <  1:109380/205‑1 (MQ=255)
                                                                                                                                                                                  ggtaaGAGGGCAGAGATATAGAGTTGCGTC‑AACAAGACGTAAACCACCGGCAAAAAGATACAATCCAGGAGCGATGCCCAGCCGCCCATAAATCCGACGTGCGGGTTAACCGTC‑TTTTGCGCGTAGGGATAGGCCGAACCGGCCTCCGGAAACTGGCGACCCAGTTTGCCGTTGCTGAGAGCGTGAAACATCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGAc                                                  <  1:69856/248‑1 (MQ=255)
                                                                                                                                                                                           gCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGAC‑TTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCaa                                          >  1:5515/1‑250 (MQ=255)
                                                                                                                                                                                               agagaAAGTTTTGCGCC‑AACAAGACGTTAAGCATCGGCAAAAAGAGACAAGCCAGCAGCGATGACCAGCCGACCAAAAATCCGACGGGCGGGTTAATCGAC‑TTTTGCGCGTAGGTATAGGCCGAACCGGCCCCCGGAAACGGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTa                                     <  1:368932/251‑1 (MQ=255)
                                                                                                                                                                                                    tagatTTTCGCC‑AACAAGACGTTAATCATCGGCACAAAGAGACAATCCAGCAGCGATGACCAGCCGCCCATAAATCCGACGTGCGGGTTAATCGAC‑TTTTGCGCGTCGGTATAGGCCGAACCGGCCTCCGGAACCTGGAGAACCAGTTTGCCGTAGATGATAGCGGTAACCAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGCCACTACAAAAAGTATCaa                                 >  2:540487/1‑250 (MQ=255)
                                                                                                                                                                                                              cc‑AACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGAC‑TTTTGCGCGTCGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATc                      >  1:49716/1‑251 (MQ=255)
                                                                                                                                                                                                                aaacaaGACGTTAATCAACGGAAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGGGCGGGTTAATCGAC‑TTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCg                     <  2:5515/250‑1 (MQ=255)
                                                                                                                                                                                                                  acaaGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGAC‑TTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGc                   >  1:107914/1‑251 (MQ=255)
                                                                                                                                                                                                                            aaTCATCGGCAAAAAGAGACAATCAAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATGGAC‑TTTTGCGCGTAGGTATAGGCCGAACCGGCCGCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAgg         <  2:568841/251‑1 (MQ=255)
                                                                                                                                                                                                                              acagcGGCAAAAAGAGACAACCCAGCATCGATGACCAGCCGACCATAAACCCGACGTGCGGGTTAACCGAC‑TTGTGCGCGGAGGTATAGGCCGAACCGGCCTCCGGAAACGGGCGAACCAGTTTGCCGTAGCTGAGAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGTACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGcc       <  1:533694/247‑1 (MQ=255)
                                                                                                                                                                                                                                   ggCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGAC‑TTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGCTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGAcc  >  1:555747/1‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                               
GCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCC‑AACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGAC‑TTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACC  >  NC_000913/1358621‑1359097

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: ‑