Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,134,222 G→A 6.4% G141S (GGC→AGC)  flgF → flagellar component of cell‑proximal portion of basal‑body rod

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,134,2220GA 6.4% 192.4 / 3.2 62G141S (GGC→AGC) flgFflagellar component of cell‑proximal portion of basal‑body rod
Reads supporting (aligned to +/- strand):  ref base G (33/25);  new base A (2/2);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

GCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCAGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGG  >  NC_000913/1134078‑1134370
                                                                                                                                                |                                                                                                                                                    
gCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCAGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGcac                                                                                                                                                  <  6:205375/149‑1 (MQ=255)
    aCGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCAGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAAt                                                                                                                                               <  8:322753/148‑1 (MQ=255)
            gAAGGGTATACGCGTAATGGCAGCATTCAGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCt                                                                                                                                      <  5:231982/149‑1 (MQ=255)
                   ataCGCGTAATGGCAGCATTCAGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCTAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCg                                                                                                                               >  1:139451/1‑149 (MQ=255)
                      cgcgTAATGGCAGCATTCAGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATTGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGc                                                                                                                            >  1:132707/1‑149 (MQ=255)
                         gTAATGGCAGCATTCAGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCTGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGAt                                                                                                                         >  7:237483/1‑149 (MQ=255)
                               gcagcaTTCAGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGc                                                                                                                    >  6:328386/1‑148 (MQ=255)
                               gcagcaTTCAGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCAGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGc                                                                                                                    >  2:351134/1‑148 (MQ=255)
                                    aTTCAGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATAc                                                                                                              >  7:325254/1‑149 (MQ=255)
                                     ttCAGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCTAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACg                                                                                                             <  2:139451/149‑1 (MQ=255)
                                        aGGGTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGTGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGATATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGtt                                                                                                          <  4:287718/149‑1 (MQ=255)
                                          gTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTgc                                                                                                        >  8:177086/1‑149 (MQ=255)
                                          gTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTgc                                                                                                        >  6:8215/1‑149 (MQ=255)
                                          gTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTgc                                                                                                        <  3:349289/149‑1 (MQ=255)
                                               tCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAg                                                                                                   >  6:251326/1‑149 (MQ=255)
                                                  cACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAg                                                                                                >  4:21519/1‑149 (MQ=255)
                                                       ggCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGt                                                                                           >  5:307013/1‑149 (MQ=255)
                                                        gCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGTCGCTCAATCCGGGCTATCCGGCAAATACGCTTGCGCCAGAACGGctac                                                                                          >  8:166673/1‑146 (MQ=255)
                                                             tGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTaaa                                                                                     >  5:131666/1‑149 (MQ=255)
                                                             tGACAATTCACGGGAATCCGGTGATAGGCGAGGATGGGCCAAATGATTTGCAGGAAGGGGAGGAAATCACTATTGCTGCCGATGGCACAAACTCTGCGCTCAATCCGGGCGATCCGGAAAATACGGTTGCGCCAGAAGGGCTTCTaaaa                                                                                     >  8:336201/1‑149 (MQ=255)
                                                               aCAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAAc                                                                                    <  5:251326/148‑1 (MQ=255)
                                                               aCAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACt                                                                                   >  3:136567/1‑149 (MQ=255)
                                                               aCAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACt                                                                                   >  7:163666/1‑149 (MQ=255)
                                                               aCAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCCGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTaaa                                                                                     <  1:351134/147‑1 (MQ=255)
                                                                 aaTTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCGAGTAGGGCGGCTTAAACTgg                                                                                 >  4:215269/1‑149 (MQ=255)
                                                                 aaTTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTgg                                                                                 <  3:21519/149‑1 (MQ=255)
                                                                 aaTTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCAAGTAGGGCGTCTTAAACTgg                                                                                 >  1:125634/1‑149 (MQ=255)
                                                                  aTTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAAtctc                                                                                                                                            >  1:338019/1‑89 (MQ=255)
                                                                  aTTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAAtctc                                                                                                                                            <  2:338019/89‑1 (MQ=255)
                                                                   ttCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGCAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGCGCGATCCGTCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGt                                                                                >  2:331367/1‑148 (MQ=255)
                                                                             tCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGCTCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACgg                                                                     >  2:35944/1‑149 (MQ=255)
                                                                             tCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACgg                                                                     >  6:234812/1‑149 (MQ=255)
                                                                                 gTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCggg                                                                                                                             >  8:71465/1‑89 (MQ=255)
                                                                                 gTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCggg                                                                                                                             <  7:71465/89‑1 (MQ=255)
                                                                                 gTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCggg                                                                                                                             >  6:110709/1‑89 (MQ=255)
                                                                                 gTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCggg                                                                                                                             <  5:110709/89‑1 (MQ=255)
                                                                                     tAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGaa                                                             <  1:48557/149‑1 (MQ=255)
                                                                                          gAGGCTGGGCCAATTGCTGTGCCGTAAGGGGCGGAAATCCCTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGGTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCa                                                        <  2:67795/149‑1 (MQ=255)
                                                                                                gggCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCgg                                                  >  6:313572/1‑149 (MQ=255)
                                                                                                 ggCCAATTTCTGTGCCGGAAGGGGAGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCTGGCGATCCGGCAAATACGGTTGCGCCATTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGcgcg                                                   <  6:128992/147‑1 (MQ=255)
                                                                                                 ggCCAATTGCTGTGCCGGAAGGGGCGGAACTCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAgc                                                      <  5:59547/144‑1 (MQ=255)
                                                                                                 ggCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAgc                                                      >  6:59547/1‑144 (MQ=255)
                                                                                                 ggCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGcgcg                                                   <  6:131666/147‑1 (MQ=255)
                                                                                                    cAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTgac                                              <  8:310040/149‑1 (MQ=255)
                                                                                                          cTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATAGCTCAATGTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTTAACTGGt                                                                                <  3:2998/109‑1 (MQ=255)
                                                                                                          cTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATAGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGt                                                                                <  3:7839/109‑1 (MQ=255)
                                                                                                          cTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATAGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGt                                                                                >  4:7839/1‑109 (MQ=255)
                                                                                                          cTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATAGCACAATCTCGGCGCTCAATCCCGGCGATCCGGCAAATACGGATGCGCCAGTAGCGCGTCTTAAACTGGt                                                                                >  4:2998/1‑109 (MQ=255)
                                                                                                           tgtgCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAAACCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCa                                       >  4:262359/1‑149 (MQ=255)
                                                                                                            gtgCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCgg                                                                                                                              <  3:34322/61‑1 (MQ=255)
                                                                                                            gtgCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCgg                                                                                                                              >  4:34322/1‑61 (MQ=255)
                                                                                                               ccGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCAttt                                    <  5:328386/148‑1 (MQ=255)
                                                                                                                 ggAAGGGGCGGTAATCACTATTGCGGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTGTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTc                                 <  1:331367/149‑1 (MQ=255)
                                                                                                                 ggAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTc                                 <  6:307013/149‑1 (MQ=255)
                                                                                                                   aaGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCCCGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGt                               >  7:108593/1‑149 (MQ=255)
                                                                                                                       ggCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGGCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTaa                           <  2:132707/149‑1 (MQ=255)
                                                                                                                            aaaTCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAgcg                        >  1:300703/1‑147 (MQ=255)
                                                                                                                                  cTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAAcc                >  7:234899/1‑149 (MQ=255)
                                                                                                                                     ttGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAg             <  5:234812/149‑1 (MQ=255)
                                                                                                                                     ttGCTGACGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGACAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGAAGCGAAGTGCAGCGCGGTGACGAAGGCATTTTTCGTTTAAGCGCAGAAAccc               <  7:166673/147‑1 (MQ=255)
                                                                                                                                         tgcCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCa         >  7:95863/1‑149 (MQ=255)
                                                                                                                                         tgcCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAAAACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCa         >  1:115645/1‑149 (MQ=255)
                                                                                                                                                ggCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTgg  >  8:374572/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCAGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAGTGCAGCGCGGTGACGACGGCATTTTTCGTTTAAGCGCAGAAACCCAGGCCACGCGTGG  >  NC_000913/1134078‑1134370

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