Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,041,221 G→A 8.3% A71A (GCC→GCT mobB ← molybdopterin‑guanine dinucleotide biosynthesis protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,041,2210GA 8.3% 736.3 / 41.1 240A71A (GCC→GCTmobBmolybdopterin‑guanine dinucleotide biosynthesis protein B
Reads supporting (aligned to +/- strand):  ref base G (118/102);  new base A (6/14);  total (124/116)
Fisher's exact test for biased strand distribution p-value = 6.01e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.76e-02

Only 100 of 250 total aligned reads displayed.

TAATATCCAGTAATGCGACATCAAGATTAAGCGGCACATCACTGGCTACAGCAATAACATGCCTGTCTATCACTAATTCTTCAGGTCGATGTCCGGCTCCATCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTCAAAAGCGTAGTTTTTCCGGTACCACTCCACGCGGCAAAGGCGAGTAACGGTATCATCGTTTTTCCTGCCATCGGGCAAGCTCCTCTGGCGTATTCACGTTAACAAATGCATCTTTATGA  >  NC_000913/4040965‑4041481
                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                    
                                  cacaTCACTGGCTACAGCAATAACATGCCTGTCTATCACTAATTCTTCAGGTCGATGTCCGGCTCCATCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTTCCGCGCCAGACTTCCGCAGCTCATAGCTATCtttcc                                                                                                                                                                                                                  >  4:81638/1‑273 (MQ=255)
                                        acgggCTACAGCAATAACCGTCCGGTCTATCACTAATTGTGCAGGTCGATGTACGGCTCACTCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGGTTTGAGGTAGCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGc                                                                                                                                                                                                             <  4:92995/271‑1 (MQ=255)
                                                cAGCAATAACATGCCTGTCGATCACTAATTCTTCAGGTCGATGTCCGGCTCCATCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAAc                                                                                                                                                                                                     <  4:268723/274‑1 (MQ=255)
                                                 aGCAATAACATGCCGGTCTATAACTAATTCTCCAGGTCGATTTCCGGCTCCATCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACAt                                                                                                                                                                                                   <  4:117417/275‑1 (MQ=255)
                                                      tAACATGCCTGTCTATCACTAATTCTTCAGGTCGATGTCCGGCTCCATCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTACTTTGGCCTGGCTTTTCAccatccatt                                                                                                                                                                                              >  4:20182/1‑274 (MQ=255)
                                                                ggcTATAACTAATTCTTCAGGTTGATTTCCGGCTCCATCGCGAAACAGCACAACCTTTGCGATCGCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCg                                                                                                                                                                                    <  4:159032/273‑1 (MQ=255)
                                                                   tATCACTAATTCTTCAGGTCGATGTCCGGCTCCATCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCt                                                                                                                                                                                                            >  1:23502/1‑248 (MQ=255)
                                                                                   ggTCGATGTCCGGCTCCATCGCGAACCAGCACAATCTGTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTgg                                                                                                                                                                 <  2:20182/275‑1 (MQ=255)
                                                                                         tGTCCGGCTCCATCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCt                                                                                                                                                                                                            >  2:40264/1‑226 (MQ=255)
                                                                                             cGGCTCCATCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATc                                                                                                                                                                                                        <  4:40264/226‑1 (MQ=255)
                                                                                                                 acaATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAAc                                                                                                                                                                                                                                                                 >  6:1207447/1‑149 (MQ=255)
                                                                                                                  caATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACg                                                                                                                                                                                                                                                                >  5:1058759/1‑149 (MQ=255)
                                                                                                                   aaTCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGa                                                                                                                                                                                                                                                               <  5:44777/149‑1 (MQ=255)
                                                                                                                   aaTCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGa                                                                                                                                                                                                                                                               >  5:850428/1‑149 (MQ=255)
                                                                                                                    aTCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCTATCGTTGCCAGCCAGCAACGAt                                                                                                                                                                                                                                                              <  6:707226/149‑1 (MQ=255)
                                                                                                                     tCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGCTCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCCAGCCAGCAACGATg                                                                                                                                                                                                                                                             <  6:946872/149‑1 (MQ=255)
                                                                                                                     tCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCCAGCCAGCAACGATg                                                                                                                                                                                                                                                             <  6:1105593/149‑1 (MQ=255)
                                                                                                                     tCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCCAGCCAGCAACGATg                                                                                                                                                                                                                                                             <  7:827157/149‑1 (MQ=255)
                                                                                                                      cTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCCAGCCAGCAACGATgg                                                                                                                                                                                                                                                            <  6:394856/149‑1 (MQ=255)
                                                                                                                      cTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCCAGCCAGCAACGATgg                                                                                                                                                                                                                                                            <  8:281411/149‑1 (MQ=255)
                                                                                                                       tttGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCCAGCCAGCAACGATggg                                                                                                                                                                                                                                                           <  6:1371865/149‑2 (MQ=255)
                                                                                                                       tttGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAAGCCCATCGTTGCTGGCTGGCAACGATGGt                                                                                                                                                                                                                                                           <  6:896134/149‑1 (MQ=255)
                                                                                                                        ttGCGTTCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATCCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGtt                                                                                                                                                                                                                                                          <  5:555418/149‑1 (MQ=255)
                                                                                                                         tGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGttt                                                                                                                                                                                                                                                         <  6:1396676/149‑1 (MQ=255)
                                                                                                                            gATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTtgtg                                                                                                                                                                                                                                                      >  6:77678/1‑149 (MQ=255)
                                                                                                                               ctctTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGc                                                                                                                                                                                                                           >  2:380977/1‑173 (MQ=255)
                                                                                                                                tctTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGNTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTCTCAACATCCCTATCATGGTTCGTANGCTTAATCAGCCCTGGACGGTTCCCTCTGGNNCATAATGCCGGGATCAGTTTTTTCAAAAg                                                                                                                    >  3:48216/1‑275 (MQ=255)
                                                                                                                                    cATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCa                                                                                                                                                                                                                                              >  7:697103/1‑149 (MQ=255)
                                                                                                                                    cATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCa                                                                                                                                                                                                                                              >  5:522836/1‑149 (MQ=255)
                                                                                                                                       gCTTAAACCCTTCTACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGcc                                                                                                                                                                                                                                           >  8:365687/1‑149 (MQ=255)
                                                                                                                                         ttAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCtt                                                                                                                                                                                                                                         >  5:1079004/1‑149 (MQ=255)
                                                                                                                                         ttAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTCCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCtt                                                                                                                                                                                                                                         >  7:1059864/1‑149 (MQ=255)
                                                                                                                                          tAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCCCATACTGCCGGGATCAGTTTTTTCAAAATCTTAGtttt                                                                                                           >  2:256994/1‑274 (MQ=255)
                                                                                                                                           aaaCCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATg                                                                                                                                                                                                                                                             <  6:45487/127‑1 (MQ=255)
                                                                                                                                           aaaCCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgc                                                                                                                                                                                                                                       <  7:430977/149‑1 (MQ=255)
                                                                                                                                           aaaCCCTTCGACCAGAATCAAAGCCAGCTTTGAGGTATCCATTCGACTTCCGAGAAACTGTAGTTCCAGCTCTTTTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTCAAAAGCGTAGTTTTTc                                                                                                         <  1:48216/275‑1 (MQ=255)
                                                                                                                                            aaCCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgcg                                                                                                                                                                                                                                      <  5:1324181/149‑1 (MQ=255)
                                                                                                                                             aCCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCCAGCCAGCAACGAt                                                                                                                                                                                                                                                              <  5:868352/124‑1 (MQ=255)
                                                                                                                                                cTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGc                                                                                                                                                                                                                                  >  8:541604/1‑149 (MQ=255)
                                                                                                                                                 ttCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCt                                                                                                                                                                                                                                 >  5:1248595/1‑149 (MQ=255)
                                                                                                                                                       cAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGc                                                                                                                                                                                                                           >  5:99252/1‑149 (MQ=255)
                                                                                                                                                            tCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCt                                                                                                                                                                                                                      >  7:1176167/1‑149 (MQ=255)
                                                                                                                                                            tCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCCAGCCAGCAACGATg                                                                                                                                                                                                                                                             <  6:399345/110‑1 (MQ=255)
                                                                                                                                                            tCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCCAGCCAGCAACGATg                                                                                                                                                                                                                                                             >  5:399345/1‑110 (MQ=255)
                                                                                                                                                             cAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCtt                                                                                                                                                                                                                     <  6:1196094/149‑1 (MQ=255)
                                                                                                                                                               aacccAGCTTTGAGTGATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCGTCTTCGTCTGGGGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGAAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATTCCTCTGGCGCATAATGCCGGGATCAGTTTTTTCAAAAGCGTAGTTTTTCCGGTACCACTCCACGCGGCa                                                                                     <  2:81638/272‑1 (MQ=255)
                                                                                                                                                               aaTCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTg                                                                                                                                                                                                                   <  5:224104/149‑1 (MQ=255)
                                                                                                                                                                 tCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGcc                                                                                                                                                                                                                 >  5:1437801/1‑149 (MQ=255)
                                                                                                                                                                  ccAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCt                                                                                                                                                                                                                <  8:817996/149‑1 (MQ=255)
                                                                                                                                                                    aGCTTTGAGGTATCCATTCGACTTGCGAGGAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTTCTTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTNTCCACATCCATATCATgg                                                                                                                                                                                        >  4:352432/1‑171 (MQ=255)
                                                                                                                                                                       tttGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCtt                                                                                                                                                                                                           >  6:514875/1‑149 (MQ=255)
                                                                                                                                                                          gAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATc                                                                                                                                                                                                        <  6:1339229/149‑1 (MQ=255)
                                                                                                                                                                             gTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAAc                                                                                                                                                                                                     >  6:262533/1‑149 (MQ=255)
                                                                                                                                                                                   aTTCGACTTGCGAGAAACTTTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGc                                                                                                                                                                                                                                  <  6:250393/114‑1 (MQ=255)
                                                                                                                                                                                   aTTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCat                                                                                                                                                                                               >  8:786722/1‑149 (MQ=255)
                                                                                                                                                                                    ttCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCata                                                                                                                                                                                              >  8:1163636/1‑149 (MQ=255)
                                                                                                                                                                                    ttCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCata                                                                                                                                                                                              >  7:468544/1‑149 (MQ=255)
                                                                                                                                                                                          ttGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCa                                                                                                                                                                                                                                    >  1:15298/1‑105 (MQ=255)
                                                                                                                                                                                          ttGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTc                                                                                                                                                                                                                                >  4:308924/1‑109 (MQ=255)
                                                                                                                                                                                          ttGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATgg                                                                                                                                                                                        >  7:1369948/1‑149 (MQ=255)
                                                                                                                                                                                          ttGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATgg                                                                                                                                                                                        >  8:95485/1‑149 (MQ=255)
                                                                                                                                                                                            gCGAGAAACTGTAGCTCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCACATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTTCGTATGCTTCATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGAACAGTTTTTTCAAAAGCGTAGTTTTTCCGGTACCACTTCACGCGGCCAAGGCGAGTAACGGTATCATCGTTTTTc                                                         >  2:124012/1‑274 (MQ=255)
                                                                                                                                                                                              gagaAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTc                                                                                                                                                                                                                                <  3:15298/105‑1 (MQ=255)
                                                                                                                                                                                              gagaAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAg                                                                                                                                                                                                                            <  2:308924/109‑1 (MQ=255)
                                                                                                                                                                                              gagaAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCg                                                                                                                                                                                    >  8:468129/1‑149 (MQ=255)
                                                                                                                                                                                                   aCTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGc                                                                                                                                                                               <  6:1097015/149‑1 (MQ=255)
                                                                                                                                                                                                   aCTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGc                                                                                                                                                                               <  7:547926/149‑1 (MQ=255)
                                                                                                                                                                                                     tGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCtt                                                                                                                                                                             <  5:703738/149‑1 (MQ=255)
                                                                                                                                                                                                         gATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATc                                                                                                                                                                         <  6:1046149/149‑1 (MQ=255)
                                                                                                                                                                                                           tCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATTAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAg                                                                                                                                                                       <  8:873948/149‑1 (MQ=255)
                                                                                                                                                                                                           tCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAg                                                                                                                                                                       >  5:924303/1‑149 (MQ=255)
                                                                                                                                                                                                             cAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGcc                                                                                                                                                                     <  5:1235387/149‑1 (MQ=255)
                                                                                                                                                                                                             cAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGcc                                                                                                                                                                     <  8:138815/149‑1 (MQ=255)
                                                                                                                                                                                                               gCTCTTCTGCGTCTAGTGTTTGCGTCATCAACGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCGCATAGCTATCTTTGCCTGTCTTTTCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCt                                                                                                                                                                   <  7:902450/149‑1 (MQ=255)
                                                                                                                                                                                                                      ttcGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATg                                                                                                                                                                                         >  8:414941/1‑120 (MQ=255)
                                                                                                                                                                                                                      ttcGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGa                                                                                                                                                            <  8:763994/149‑1 (MQ=255)
                                                                                                                                                                                                                       tcGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGAt                                                                                                                                                           <  7:1402009/149‑1 (MQ=255)
                                                                                                                                                                                                                         gTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGc                                                                                                                                                                               >  5:1454529/1‑127 (MQ=255)
                                                                                                                                                                                                                         gTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGc                                                                                                                                                                               <  6:1454529/127‑1 (MQ=255)
                                                                                                                                                                                                                            tGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTg                                                                                                                                                                  >  2:344298/1‑137 (MQ=255)
                                                                                                                                                                                                                                  ttCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGgcgc                                                                                                                                                >  5:1374518/1‑149 (MQ=255)
                                                                                                                                                                                                                                   tCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCa                                                                                                                                               >  5:1392555/1‑149 (MQ=255)
                                                                                                                                                                                                                                   cccGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCa                                                                                                                                               <  8:1059864/148‑1 (MQ=255)
                                                                                                                                                                                                                                       tcatcaAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATcc                                                                                                                                                         <  7:261229/135‑1 (MQ=255)
                                                                                                                                                                                                                                         atcaAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATc                                                                                                                                                          <  5:175208/132‑1 (MQ=255)
                                                                                                                                                                                                                                         atcaAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGc                                                                                                                                         >  8:629373/1‑149 (MQ=255)
                                                                                                                                                                                                                                         atcaAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGc                                                                                                                                         >  7:1258117/1‑149 (MQ=255)
                                                                                                                                                                                                                                            aaGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCgg                                                                                                                                      <  6:850428/149‑1 (MQ=255)
                                                                                                                                                                                                                                              ggCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGa                                                                                                                                    <  8:60269/149‑1 (MQ=255)
                                                                                                                                                                                                                                               gCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGAt                                                                                                                                   <  8:151883/149‑1 (MQ=255)
                                                                                                                                                                                                                                                cccATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATc                                                                                                                                  <  7:210523/149‑1 (MQ=255)
                                                                                                                                                                                                                                                accATCGTTGCGAGCCAGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATc                                                                                                                                  >  5:1098089/2‑149 (MQ=255)
                                                                                                                                                                                                                                                accATCGTTGCCAGCCAGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGa                                                                                                                                                                <  8:255797/118‑1 (MQ=255)
                                                                                                                                                                                                                                                accATCGTTGCCAGCCAGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATc                                                                                                                                  >  8:939736/2‑149 (MQ=255)
                                                                                                                                                                                                                                                 ccATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCa                                                                                                                                 <  7:377821/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   aTCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgcgc                                                                                                                                                                                                                                     <  7:626455/47‑1 (MQ=255)
                                                                                                                                                                                                                                                       ttGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGttttt                                                                                                                           <  8:1257974/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        tGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGtttttt                                                                                                                          <  8:11917/149‑1 (MQ=255)
                                                                                                                                                                                                                                                             gctggcAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTCaaca                                                                                                                     <  8:1331917/149‑3 (MQ=255)
                                                                                                                                                                                                                                                               tggcAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTCAAAAGc                                                                                                                   >  8:1020295/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                    
TAATATCCAGTAATGCGACATCAAGATTAAGCGGCACATCACTGGCTACAGCAATAACATGCCTGTCTATCACTAATTCTTCAGGTCGATGTCCGGCTCCATCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTCAAAAGCGTAGTTTTTCCGGTACCACTCCACGCGGCAAAGGCGAGTAACGGTATCATCGTTTTTCCTGCCATCGGGCAAGCTCCTCTGGCGTATTCACGTTAACAAATGCATCTTTATGA  >  NC_000913/4040965‑4041481

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