Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,041,216 T→C 7.9% Q73R (CAG→CGG) ‡ 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,2160TC 7.9% 771.4 / 38.8 240Q73R (CAG→CGG) ‡mobBmolybdopterin‑guanine dinucleotide biosynthesis protein B
Reads supporting (aligned to +/- strand):  ref base T (118/102);  new base C (5/14);  total (124/116)
Fisher's exact test for biased strand distribution p-value = 3.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

Only 100 of 247 total aligned reads displayed.

TAATATCCAGTAATGCGACATCAAGATTAAGCGGCACATCACTGGCTACAGCAATAACATGCCTGTCTATCACTAATTCTTCAGGTCGATGTCCGGCTCCATCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTCAAAAGCGTAGTTTTTCCGGTACCACTCCACGCGGCAAAGGCGAGTAACGGTATCATCGTTTTTCCTGCCATCGGGCAAGCTCCTCTGGCGTATTCACGTTAACAAATGCATCTTTATGA  >  NC_000913/4040965‑4041481
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                         
             tGCGACATCAAGATTAAGCGGCACATCACTGGCTACAGCAATAACATGCCTGTCTATCACTAATTCTTCAGGTCGATGTCCGGCTCCATCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgc                                                                                                                                                                                                                                       >  2:268723/1‑275 (MQ=255)
                                  cacaTCACTGGCTACAGCAATAACATGCCTGTCTATCACTAATTCTTCAGGTCGATGTCCGGCTCCATCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTTCCGCGCCAGACTTCCGCAGCTCATAGCTATCtttcc                                                                                                                                                                                                                  >  4:81638/1‑273 (MQ=255)
                                        acgggCTACAGCAATAACCGTCCGGTCTATCACTAATTGTGCAGGTCGATGTACGGCTCACTCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGGTTTGAGGTAGCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGc                                                                                                                                                                                                             <  4:92995/271‑1 (MQ=255)
                                                cAGCAATAACATGCCTGTCGATCACTAATTCTTCAGGTCGATGTCCGGCTCCATCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAAc                                                                                                                                                                                                     <  4:268723/274‑1 (MQ=255)
                                                                ggcTATAACTAATTCTTCAGGTTGATTTCCGGCTCCATCGCGAAACAGCACAACCTTTGCGATCGCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCg                                                                                                                                                                                    <  4:159032/273‑1 (MQ=255)
                                                                gTCTATCACTAATTCTTCAGGTCGATGTCCGGCTCCATCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTATTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTGGCGCACCTCATAGCTATCTTTACCTGGCTTATCAACATCCATATTCTGGTGc                                                                                                                                                                                     >  1:40950/1‑274 (MQ=255)
                                                                   tATCACTAATTCTTCAGGTCGATGTCCGGCTCCATCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCt                                                                                                                                                                                                            >  1:23502/1‑248 (MQ=255)
                                                                   tATCACGCATTCTTCAGGTCGATGTACGGCTCCATCGCGAAACAGCACAAGCTTTGCGCTCCCTTCATGCTTCAGCCCTTCCACCCGCATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCTAGAAACTGTAGATCCAGCTCTTCTTCGTCAGGTGATTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCCACGATGGTTTGTGCCTCCCCAGCCTTGCGCAGCTCAGTGCTTACTTTGCCTGGCTTAACAACATCCaacacactggcac                                                                                                                                                                                    >  4:204415/1‑260 (MQ=255)
                                                                     tCACTAATTCGTCAGGTCGATGTCAGGCTCCATCGCGAAACAGCACAATCTTGGCGATCTCTTCATGCTTAACCCCTTCTGCCAGAATCAAATCCAGCCTTGAGGTATCCATTCGACTTGCGAGAAACGGTAGATCCAGCTTTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGc                                                                                                                                                                               <  2:447/275‑1 (MQ=255)
                                                                                   ggTCGATGTCCGGCTCCATCGCGAACCAGCACAATCTGTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTgg                                                                                                                                                                 <  2:20182/275‑1 (MQ=255)
                                                                                             cGGCTCCATCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATc                                                                                                                                                                                                        <  4:40264/226‑1 (MQ=255)
                                                                                                       gcgAAACAGCAGAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTCGTGTTTCCGTCACCAAGGCCCATCGTTgct                                                                                                                                                                                                                                                                           <  8:1220572/149‑1 (MQ=255)
                                                                                                       gcgAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTgct                                                                                                                                                                                                                                                                           <  8:1064003/149‑1 (MQ=255)
                                                                                                                cacaATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCaa                                                                                                                                                                                                                                                                  <  6:910910/149‑1 (MQ=255)
                                                                                                                  caATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACg                                                                                                                                                                                                                                                                >  5:1058759/1‑149 (MQ=255)
                                                                                                                   aaTCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGa                                                                                                                                                                                                                                                               <  5:44777/149‑1 (MQ=255)
                                                                                                                     tCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCCAGCCAGCAACGATg                                                                                                                                                                                                                                                             <  7:827157/149‑1 (MQ=255)
                                                                                                                      cTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCCAGCCAGCAACGATgg                                                                                                                                                                                                                                                            <  6:394856/149‑1 (MQ=255)
                                                                                                                      cTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCCAGCCAGCAACGATgg                                                                                                                                                                                                                                                            <  8:814858/149‑1 (MQ=255)
                                                                                                                        ttGCGTTCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATCCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGtt                                                                                                                                                                                                                                                          <  5:555418/149‑1 (MQ=255)
                                                                                                                         tGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGttt                                                                                                                                                                                                                                                         <  6:1396676/149‑1 (MQ=255)
                                                                                                                            gATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTtgtg                                                                                                                                                                                                                                                      >  6:736576/1‑149 (MQ=255)
                                                                                                                             aTCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGc                                                                                                                                                                                                                                                     >  6:1236856/1‑149 (MQ=255)
                                                                                                                               ctctTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGc                                                                                                                                                                                                                           >  2:380977/1‑173 (MQ=255)
                                                                                                                                tctTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGNTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTCTCAACATCCCTATCATGGTTCGTANGCTTAATCAGCCCTGGACGGTTCCCTCTGGNNCATAATGCCGGGATCAGTTTTTTCAAAAg                                                                                                                    >  3:48216/1‑275 (MQ=255)
                                                                                                                                tctTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCcgc                                                                                                                                                                                                                                                  <  7:1421270/149‑1 (MQ=255)
                                                                                                                                     aTGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAg                                                                                                                                                                                                                                             >  8:320932/1‑149 (MQ=255)
                                                                                                                                         ttAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCtt                                                                                                                                                                                                                                         >  5:1079004/1‑149 (MQ=255)
                                                                                                                                           aaaCCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATg                                                                                                                                                                                                                                                             <  6:45487/127‑1 (MQ=255)
                                                                                                                                           aaaCCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgc                                                                                                                                                                                                                                       <  7:320932/149‑1 (MQ=255)
                                                                                                                                            aaCCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgcg                                                                                                                                                                                                                                      <  5:1324181/149‑1 (MQ=255)
                                                                                                                                             aCCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCCAGCCAGCAACGAt                                                                                                                                                                                                                                                              >  6:868352/1‑124 (MQ=255)
                                                                                                                                             aCCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCCAGCCAGCAACGAt                                                                                                                                                                                                                                                              <  5:868352/124‑1 (MQ=255)
                                                                                                                                              cccTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCa                                                                                                                                                                                                                                    >  8:945307/1‑149 (MQ=255)
                                                                                                                                               ccTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCATCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAg                                                                                                                                                                                                                                   >  7:562283/1‑149 (MQ=255)
                                                                                                                                                 ttCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCt                                                                                                                                                                                                                                 >  5:1248595/1‑149 (MQ=255)
                                                                                                                                                 ttCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCt                                                                                                                                                                                                                                 >  7:138815/1‑149 (MQ=255)
                                                                                                                                                     aCCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATa                                                                                                                                                                                                                             <  8:562283/149‑1 (MQ=255)
                                                                                                                                                       cAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGc                                                                                                                                                                                                                           >  5:99252/1‑149 (MQ=255)
                                                                                                                                                        aGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCt                                                                                                                                                                                                                          <  8:1347981/149‑1 (MQ=255)
                                                                                                                                                            tCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCCAGCCAGCAACGATg                                                                                                                                                                                                                                                             <  6:399345/110‑1 (MQ=255)
                                                                                                                                                             cAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCtt                                                                                                                                                                                                                     <  6:1196094/149‑1 (MQ=255)
                                                                                                                                                              aaaTCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCttt                                                                                                                                                                                                                    >  5:855841/1‑149 (MQ=255)
                                                                                                                                                               aacccAGCTTTGAGTGATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCGTCTTCGTCTGGGGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGAAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATTCCTCTGGCGCATAATGCCGGGATCAGTTTTTTCAAAAGCGTAGTTTTTCCGGTACCACTCCACGCGGCa                                                                                     <  2:81638/272‑1 (MQ=255)
                                                                                                                                                               aaTCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTg                                                                                                                                                                                                                   <  5:224104/149‑1 (MQ=255)
                                                                                                                                                                 tCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGcc                                                                                                                                                                                                                 >  5:1437801/1‑149 (MQ=255)
                                                                                                                                                                  ccgtcTTTGAGGTTTCCATTCGACTGGCTGGAAACGTTAGATCCAGCGCTTCTTGTTCTGGTGTTTCCGGCATCAAGGCCCATCGTTGCTGGCGGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTCAAAAGCGTAGTTTTTCCGGTACCACTCCACGCGGCAAAg                                                                                  <  4:12487/271‑1 (MQ=255)
                                                                                                                                                                  ccAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCt                                                                                                                                                                                                                >  8:1105916/1‑149 (MQ=255)
                                                                                                                                                                    aGCTTTGAGGTATCCATTCGACTTGCGAGGAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTTCTTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTNTCCACATCCATATCATgg                                                                                                                                                                                        >  4:352432/1‑171 (MQ=255)
                                                                                                                                                                       tttGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCtt                                                                                                                                                                                                           >  6:514875/1‑149 (MQ=255)
                                                                                                                                                                       attGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCtt                                                                                                                                                                                                           <  8:939105/148‑1 (MQ=255)
                                                                                                                                                                         tGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTAt                                                                                                                                                                                                         <  8:458686/149‑1 (MQ=255)
                                                                                                                                                                          gAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATc                                                                                                                                                                                                        <  6:1339229/149‑1 (MQ=255)
                                                                                                                                                                          gAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATc                                                                                                                                                                                                        <  7:1064137/149‑1 (MQ=255)
                                                                                                                                                                             gTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAAc                                                                                                                                                                                                     >  6:262533/1‑149 (MQ=255)
                                                                                                                                                                                tCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCTATATc                                                                                                                                                                                                  >  8:1261378/1‑149 (MQ=255)
                                                                                                                                                                                   aTTCGACTTGCGAGAAACTTTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGc                                                                                                                                                                                                                                  <  6:250393/114‑1 (MQ=255)
                                                                                                                                                                                    ttCNACTNGCGAGAANCTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCt                                                                                                                                                                                                            <  3:20724/135‑1 (MQ=255)
                                                                                                                                                                                    ttCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCt                                                                                                                                                                                                                      >  6:1382012/1‑125 (MQ=255)
                                                                                                                                                                                    ttCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCata                                                                                                                                                                                              >  7:468544/1‑149 (MQ=255)
                                                                                                                                                                                          ttGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCa                                                                                                                                                                                                                                    >  1:15298/1‑105 (MQ=255)
                                                                                                                                                                                          ttGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTc                                                                                                                                                                                                                                >  4:308924/1‑109 (MQ=255)
                                                                                                                                                                                           tGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGt                                                                                                                                                                                       <  8:248969/149‑1 (MQ=255)
                                                                                                                                                                                            gCGAGAAACTGTAGCTCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCACATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTTCGTATGCTTCATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGAACAGTTTTTTCAAAAGCGTAGTTTTTCCGGTACCACTTCACGCGGCCAAGGCGAGTAACGGTATCATCGTTTTTc                                                         >  2:124012/1‑274 (MQ=255)
                                                                                                                                                                                             cGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGc                                                                                                                                                                                     <  7:716930/149‑1 (MQ=255)
                                                                                                                                                                                             cGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGc                                                                                                                                                                                     <  8:1436078/149‑1 (MQ=255)
                                                                                                                                                                                              gagaAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAg                                                                                                                                                                                                                            <  2:308924/109‑1 (MQ=255)
                                                                                                                                                                                                gaAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTa                                                                                                                                                                                  >  6:796676/1‑149 (MQ=255)
                                                                                                                                                                                                   aCTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGc                                                                                                                                                                               <  6:1097015/149‑1 (MQ=255)
                                                                                                                                                                                                   aCTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGc                                                                                                                                                                               <  8:697103/149‑1 (MQ=255)
                                                                                                                                                                                                     tGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCtt                                                                                                                                                                             <  5:703738/149‑1 (MQ=255)
                                                                                                                                                                                                       tAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTaa                                                                                                                                                                           <  7:95485/149‑1 (MQ=255)
                                                                                                                                                                                                         gATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGc                                                                                                                                                                               <  6:1072035/143‑1 (MQ=255)
                                                                                                                                                                                                             cAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGcc                                                                                                                                                                     <  5:1235387/149‑1 (MQ=255)
                                                                                                                                                                                                               gCTCTTCTGCGTCTAGTGTTTGCGTCATCAACGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCGCATAGCTATCTTTGCCTGTCTTTTCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCt                                                                                                                                                                   <  7:902450/149‑1 (MQ=255)
                                                                                                                                                                                                                  cttcttCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGa                                                                                                                                                                <  5:262533/149‑1 (MQ=255)
                                                                                                                                                                                                                      ttcGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGa                                                                                                                                                            >  5:935134/1‑149 (MQ=255)
                                                                                                                                                                                                                         gTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGc                                                                                                                                                                               <  6:1454529/127‑1 (MQ=255)
                                                                                                                                                                                                                         gTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGc                                                                                                                                                                               >  5:1454529/1‑127 (MQ=255)
                                                                                                                                                                                                                          tCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATccc                                                                                                                                                        >  8:407429/1‑149 (MQ=255)
                                                                                                                                                                                                                            tGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCctc                                                                                                                                                      <  6:950948/149‑1 (MQ=255)
                                                                                                                                                                                                                                gtTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGgc                                                                                                                                                  >  7:132066/1‑149 (MQ=255)
                                                                                                                                                                                                                                  ttCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGgcgc                                                                                                                                                >  5:1374518/1‑149 (MQ=255)
                                                                                                                                                                                                                                   tCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCa                                                                                                                                               >  5:1392555/1‑149 (MQ=255)
                                                                                                                                                                                                                                      gTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATaa                                                                                                                                            >  6:1114417/1‑149 (MQ=255)
                                                                                                                                                                                                                                       tcatcaAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATcc                                                                                                                                                         >  8:261229/1‑135 (MQ=255)
                                                                                                                                                                                                                                         atcaAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATc                                                                                                                                                          <  5:175208/132‑1 (MQ=255)
                                                                                                                                                                                                                                            aaGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCgg                                                                                                                                      <  5:514875/149‑1 (MQ=255)
                                                                                                                                                                                                                                               aaccATCGTTGCCAGCCAGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGAt                                                                                                                                   >  7:45208/3‑149 (MQ=255)
                                                                                                                                                                                                                                                cccATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATc                                                                                                                                  <  7:210523/149‑1 (MQ=255)
                                                                                                                                                                                                                                                accATCGTTGCGAGCCAGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATc                                                                                                                                  >  5:1098089/2‑149 (MQ=255)
                                                                                                                                                                                                                                                 ccATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCa                                                                                                                                 <  7:377821/149‑1 (MQ=255)
                                                                                                                                                                                                                                                  cATCGTTGCCAGCCAGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATGTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAg                                                                                                                                >  8:765251/1‑149 (MQ=255)
                                                                                                                                                                                                                                                  cATCGTTGCCAGCCAGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTaa                                                                                                                                                                           <  8:872816/106‑1 (MQ=255)
                                                                                                                                                                                                                                                   aTCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgcgc                                                                                                                                                                                                                                     >  8:626455/1‑47 (MQ=255)
                                                                                                                                                                                                                                                   aTCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTgcgc                                                                                                                                                                                                                                     <  7:626455/47‑1 (MQ=255)
                                                                                                                                                                                                                                                    tCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGtt                                                                                                                              <  8:468544/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      gTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGtttt                                                                                                                            >  7:1257974/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        tGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGtttttt                                                                                                                          <  8:11917/149‑1 (MQ=255)
                                                                                                                                                                                                                                                         gctggctgGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATg                                                                                                                                                                                         <  2:350875/85‑1 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                         
TAATATCCAGTAATGCGACATCAAGATTAAGCGGCACATCACTGGCTACAGCAATAACATGCCTGTCTATCACTAATTCTTCAGGTCGATGTCCGGCTCCATCGCGAAACAGCACAATCTTTGCGATCTCTTCATGCTTAAACCCTTCGACCAGAATCAAATCCAGCTTTGAGGTATCCATTCGACTTGCGAGAAACTGTAGATCCAGCTCTTCTTCGTCTGGTGTTTCCGTCATCAAGGCCCATCGTTGCTGGCTGGCAACGATGGTTTGTGCCGCGCCAGCCTTGCGCAGCTCATAGCTATCTTTGCCTGGCTTATCAACATCCATATCATGGTGCGTATGCTTAATCAGCCCTGGACGGATCCCTCTGGCGCATAATGCCGGGATCAGTTTTTTCAAAAGCGTAGTTTTTCCGGTACCACTCCACGCGGCAAAGGCGAGTAACGGTATCATCGTTTTTCCTGCCATCGGGCAAGCTCCTCTGGCGTATTCACGTTAACAAATGCATCTTTATGA  >  NC_000913/4040965‑4041481

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