Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,123,626 A→C 5.0% C299G (TGC→GGC)  cpsB ← mannose‑1‑phosphate guanyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,123,6260AC 5.0% 284.3 / 5.9 80C299G (TGC→GGC) cpsBmannose‑1‑phosphate guanyltransferase
Reads supporting (aligned to +/- strand):  ref base A (32/44);  new base C (2/2);  total (34/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.39e-01

GGTAGCGGTCGCCCGCGTCGATAGAGTCATATTTGCCCCACGGACGATACACTTCGCGATGCACCCGATGCTCATGGCGACCATCGGCTTTGATCTGCTCGACCACTTTTTTCACATCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTT‑TGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACGCTTCTTCATCCACGCGAATAAAATTGAGATCCGGATCGACGGCGCTCATCGCTTTTTCACAGGCATCGAGGATATCCGGGCGATATTTTTTCAGTTCTTCGAGATAGCGT  >  NC_000913/2123365‑2123888
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                      
ggTAGCGGTCGCCCGCGTCGATAGAGTCATATTTGCCCCACGGACGATACACTTCGCGATGCACCCGATGCTCATGGCGACCATCGGCTTTGATCTGCTCGACCACTTTTTTCACATCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTT‑TGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCCTACACATAGCTGTTTTCAGTTTTGTGATTAATCCCATCGCCGTGGCAAACGTTGccc                                                                                                                                                                                                                                                              >  2:371387/1‑272 (MQ=255)
                                                     tCGCGATGCACCCGATGCTCATGGCGACCATCGGCTTTGATCTGCTCGACCACTTTTTTCACATCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTT‑TGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTTGGCGCTCATTCCCCATAATGAAGACCAGGAGCCAACATCGCTCCa                                                                                                                                                                                                      >  2:169336/1‑275 (MQ=255)
                                                     tCGCGATGCACCCGATGCTCATGGCGACCATCGGCTTTGATCTGCTCGACCACTTTTTTCACATCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTT‑TGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCa                                                                                                                                                                                                      >  2:224649/1‑275 (MQ=255)
                                                           tGCACCCGATGCTCATGGCGACCATCGGCTTTGATCTGCTCGACCACTTTTTTCACATCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTT‑TGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCTCTCCAGCCCGc                                                                                                                                                                                                >  1:19665/1‑275 (MQ=255)
                                                                     gCTCATGGCGACCATCGGCTTTGATCTGCTCGACCACTTTTTTCACATCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTT‑TGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTTGGCGGTGTGGGCGCTGACCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCgg                                                                                                                                                                                       <  2:99811/274‑1 (MQ=255)
                                                                              gACCATCGGCTTTGATCTGCTCGACCACTTTTTTCACATCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTT‑TGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCacaaca                                                                                                                                                                             >  2:499765/1‑275 (MQ=255)
                                                                                           gttCTGCTCGACCCCTTTTTTGCCATCCTGTACCGCTTTAGGGTTGGAAATTAGCCCCGCATCTTGTGGCCT‑‑CACCACTTCCAGATCTTTCCCCCCGGCGGGGGTGACCAGGGCGGATTCAGCATACACATAGCTGTTTTCAGGTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGTCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGGGCCACCATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTa                                                                                                                                                                 <  4:499760/272‑1 (MQ=255)
                                                                                                 cTCGACCACTTTTTTCACATCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTT‑TGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCCTAATGAAGACCAGGGGCCAACATCGCTCCAGCCCGCATCCCTCGGCCACACCACAGCATCTgtcccttaaaca                                                                                                                                                           >  1:135454/1‑263 (MQ=255)
                                                                                                   gcccctcTTTTTTC‑CCTCCGTTACCCCTTTACGTGCGGCAATCAGACCCGCTACTT‑TGGTCTG‑CACCACTACCACATCTTTCACCCCGACGGTGATGACCAGGCCAGATTCAGCATACACATAGCTGTGTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGa                                                                                                                                                       <  4:224645/269‑1 (MQ=255)
                                                                                                     gccACTTTTTTCACATCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTT‑TGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTa                                                                                                                                                                 <  2:21405/263‑1 (MQ=255)
                                                                                                                    tCCTTTACCGCGTTACGGTCGGCAATCAGCACCGCATCTT‑TGGTCTG‑CACCACTACCGGATCTTTCACCCCGGCGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGa                                                                                                                                       <  2:401240/275‑1 (MQ=255)
                                                                                                                    tCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTT‑TGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGATGATCTCCCATAATGAAGACCAGGAGCCAACATCGCCCCAGCCCGCATCCCACCGCACCACAACAGCATCTGCCTTAACGTCCATGACCGCGGAATCccccaa                                                                                                                                       >  4:390889/1‑273 (MQ=255)
                                                                                                                     ccTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTT‑TGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGATTTCAGCTTTGTGATTAATCACATCGCCGTGGCAAACg                                                                                                                                                                                                                                                                    >  5:26262/1‑149 (MQ=255)
                                                                                                                       tGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTT‑TGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGtt                                                                                                                                                                                                                                                                  <  8:430842/149‑1 (MQ=255)
                                                                                                                        gTACCGCGTTACGGTCGGCAATCAGCACCGCATCTT‑TGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTg                                                                                                                                                                                                                                                                 <  7:429162/149‑1 (MQ=255)
                                                                                                                             gcgTTACGGTCGGCAATCAGCACCGCATCTT‑TGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAg                                                                                                                                                                                                                            >  4:52797/1‑181 (MQ=255)
                                                                                                                             gcgTTACGGTCGGCAATCAGCACCGCATCTT‑TGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCg                                                                                                                               >  2:347092/1‑274 (MQ=255)
                                                                                                                              cgTTACGGTCGGCAACCAGAACGCCATCTT‑TGGCGTGACCCCAC‑ACCAGTTCTTTCACCCCGACGTGGGTGACCAGGCCGGATTCGTCATACCAATAGCTGGTTTCAGTTGTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCGCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCggg                                                                                                                             <  4:184348/275‑1 (MQ=255)
                                                                                                                                ttACGGTCGGCAATCAGCACCGCATCTT‑TGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGACGTGGCAAACGTTGACCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCACTCCAGCCCGCATCCATCGGCCCCCCCACCGCATCTGCGCGCCGTTCCTTTACCCCGAGACCCCACAACCCTTCCGGGCa                                                                                                                           >  4:495171/1‑275 (MQ=255)
                                                                                                                                 tACGGTCGGCAATCAGCACCGCATCTT‑TGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCa                                                                                                                                                                                                                        <  2:52798/181‑1 (MQ=255)
                                                                                                                                 tACGGTCGGCAATCAGCACCGCATCTT‑TGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACa                                                                                                                                                                                                              <  2:228823/191‑1 (MQ=255)
                                                                                                                                        ggCAATCAGCACCGCATCTT‑TGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCgaga                                                                                                                    >  2:128781/1‑274 (MQ=255)
                                                                                                                                        ggCAATCAGCACCGCATCTT‑TGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGaa                                                                                                                   <  2:390894/275‑1 (MQ=255)
                                                                                                                                            aTCAGCACCGCATCTT‑TGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCt                                                                                                                                                                                                                                             <  8:118159/149‑1 (MQ=255)
                                                                                                                                             tCAGCACCGCATCTT‑TGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTg                                                                                                                                                                                                                                            >  7:290387/1‑149 (MQ=255)
                                                                                                                                                gCACCGCATCTT‑TGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGAtc                                                                                                                                                                                                                                         <  5:114979/149‑1 (MQ=255)
                                                                                                                                                  aCCGCATCTT‑TGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGAtctc                                                                                                                                                                                                                                       <  5:403465/149‑1 (MQ=255)
                                                                                                                                                      cATCTT‑TGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCAt                                                                                                                                                                                                                                   <  8:290387/149‑1 (MQ=255)
                                                                                                                                                         cTT‑TGGTCTG‑CACCACTACCACATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACGCTTCTTCATCCAc                                                                                                   <  2:495176/274‑1 (MQ=255)
                                                                                                                                                             tGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGaa                                                                                                                                                                                                                             <  8:93352/149‑1 (MQ=255)
                                                                                                                                                             tGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGaa                                                                                                                                                                                                                             <  5:295161/149‑1 (MQ=255)
                                                                                                                                                             tGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGaa                                                                                                                                                                                                                             <  7:97429/149‑1 (MQ=255)
                                                                                                                                                             tGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAAAGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGaa                                                                                                                                                                                                                             <  8:276617/149‑1 (MQ=255)
                                                                                                                                                                tCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGAc                                                                                                                                                                                                                          <  6:26262/149‑1 (MQ=255)
                                                                                                                                                                  tG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGc                                                                                                                                                                     >  4:192772/1‑200 (MQ=255)
                                                                                                                                                                       ccacTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTa                                                                                                                                                                 <  2:192774/200‑1 (MQ=255)
                                                                                                                                                                           tACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCaa                                                                                                                                                                                                                <  8:39598/149‑1 (MQ=255)
                                                                                                                                                                           tACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACACCGCCGAGGGCAACGTTTGCCAcggcgat                                                                                                                                                                                                                                                      <  8:442261/111‑3 (MQ=37)
                                                                                                                                                                             ccAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTccc                                                                                                                                                                                                                                     <  8:245586/126‑1 (MQ=255)
                                                                                                                                                                             ccAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTccc                                                                                                                                                                                                                                     >  7:245586/1‑126 (MQ=255)
                                                                                                                                                                             ccAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACa                                                                                                                                                                                                              >  7:282364/1‑149 (MQ=255)
                                                                                                                                                                             ccAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACACCGCCGAGGGCAACGTTTGCCAcggcgat                                                                                                                                                                                                                                                      >  8:46427/1‑107 (MQ=37)
                                                                                                                                                                              cAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACAt                                                                                                                                                                                                             >  6:219764/1‑149 (MQ=255)
                                                                                                                                                                              cAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACAt                                                                                                                                                                                                             <  5:190118/149‑1 (MQ=255)
                                                                                                                                                                               aGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCa                                                                                                                                                                                                                 <  6:142385/144‑1 (MQ=255)
                                                                                                                                                                               aGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCa                                                                                                                                                                                                                 >  5:142385/1‑144 (MQ=255)
                                                                                                                                                                               aGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTa                                                                                                                                                                 <  1:143699/192‑1 (MQ=255)
                                                                                                                                                                                 aTCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGc                                                                                                                                                                                                          >  7:319480/1‑149 (MQ=255)
                                                                                                                                                                                   cTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTc                                                                                                                                                                                                        <  7:184544/149‑1 (MQ=255)
                                                                                                                                                                                     ttCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCt                                                                                                                                                                                                                                             >  2:281937/1‑110 (MQ=255)
                                                                                                                                                                                     ttCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAGACGCTTCTTCATCCACGCGAATAAAATTGAGATCCGGATCGAc                                                                        >  2:72413/1‑275 (MQ=255)
                                                                                                                                                                                         ccccGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGAtc                                                                                                                                                                                                                                         <  4:281933/110‑1 (MQ=255)
                                                                                                                                                                                           ccGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACACCGCCGAGGGCAACGTTTGCCACGGCGAtgtgattaa                                                                                                                                                                                                                                              >  8:78586/1‑98 (MQ=21)
                                                                                                                                                                                           ccGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACACCGCCGAGGGCAACGTTTGCCACGGCGAtgtgattaa                                                                                                                                                                                                                                              <  7:78586/103‑6 (MQ=18)
                                                                                                                                                                                             gACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCgaga                                                                                                                    >  4:233077/1‑223 (MQ=255)
                                                                                                                                                                                                 gtggtgACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACg                                                                                                                <  2:233081/223‑1 (MQ=255)
                                                                                                                                                                                                  tggtgACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCggg                                                                                                                             <  1:44058/209‑1 (MQ=255)
                                                                                                                                                                                                       aCCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGgtg                                                                                                                                                                                                                                                     >  2:507341/1‑84 (MQ=255)
                                                                                                                                                                                                        ccAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGa                                                                                                                                       <  2:308831/193‑1 (MQ=255)
                                                                                                                                                                                                          aGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACGCTTCTTCATCCACGCGAATAAAATTGAGATCCGGATCGACGGCGCTCATCGCTTTTTcaca                                                   <  2:445344/275‑1 (MQ=255)
                                                                                                                                                                                                          aGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCACATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACGCTTCTTCATCCACGCGAATAAAATTGAGATCCGGATCGACGGCGCTCATCGCTTTTTcaca                                                   <  2:304102/275‑1 (MQ=255)
                                                                                                                                                                                                           ggCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTggg                                                                                                                                                                                                                                                 <  4:507336/84‑1 (MQ=255)
                                                                                                                                                                                                              cAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCacaaca                                                                                                                                                                             <  5:219764/149‑1 (MQ=255)
                                                                                                                                                                                                                  ttCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCAt                                                                                                                                                                         >  5:370387/1‑149 (MQ=255)
                                                                                                                                                                                                                        aTACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTc                                                                                                                                                                                                        >  7:348320/1‑112 (MQ=255)
                                                                                                                                                                                                                        aTACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTc                                                                                                                                                                                                        <  8:348320/112‑1 (MQ=255)
                                                                                                                                                                                                                            acaTAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACg                                                                                                                                                               >  7:363979/1‑149 (MQ=255)
                                                                                                                                                                                                                                    gTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTa                                                                                                                                                                 <  2:306750/139‑1 (MQ=255)
                                                                                                                                                                                                                                     ttttCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGcatc                                                                                                                                                                                             >  2:167196/1‑110 (MQ=255)
                                                                                                                                                                                                                                     ttttCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCacaaca                                                                                                                                                                             <  5:467121/126‑1 (MQ=255)
                                                                                                                                                                                                                                     ttttCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCacaaca                                                                                                                                                                             >  6:467121/1‑126 (MQ=255)
                                                                                                                                                                                                                                         cAGTTTTGAGATTAAGCACATAGCCGTGGCAAACGTTGGCCGCGGCGGTGTGGGGGCTGATCTCCCATAATGAAGACCAGGAGGCAACATCGATCAAGCCCGaatccatc                                                                                                                                                                                         <  4:167194/110‑1 (MQ=255)
                                                                                                                                                                                                                                          aGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGt                                                                                                                                                 <  6:187178/149‑1 (MQ=255)
                                                                                                                                                                                                                                                 tgATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGAcc                                                                                                                                                     >  8:453187/1‑138 (MQ=255)
                                                                                                                                                                                                                                                 tgATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGAcc                                                                                                                                                     <  7:453187/138‑1 (MQ=255)
                                                                                                                                                                                                                                                 tgATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCAc                                                                                                                                          >  6:35118/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    ttAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCGAGCCAGCATCCATCGGCACCACAACAGCATCTGCCGTAAGTTCCATGACCGCGTAATACACCGa                                                                                                                                       >  8:240142/1‑149 (MQ=255)
                                                                                                                                                                                                                                                           caTCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACGCTTCTTCATCCACGCGAATAAAATTGAGATCCGGATCGACGGCGCTCATCGCTTTTTCACAGGCATCGAGGATATCCGGGCGATATTTTTTCAGTTCTTCGAGATAGCGt  >  2:75946/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                  gTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCGCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACGCTTCTTCATCCACGCGAATAAAATTGAGATCCGGATCGACGGCGCTCATCGCTTTTTCACAGGCa                                               <  1:46609/223‑1 (MQ=255)
                                                                                                                                                                                                                                                                   tGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGc                                                                                                                        <  5:35118/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                      
GGTAGCGGTCGCCCGCGTCGATAGAGTCATATTTGCCCCACGGACGATACACTTCGCGATGCACCCGATGCTCATGGCGACCATCGGCTTTGATCTGCTCGACCACTTTTTTCACATCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTT‑TGGTCTG‑CACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACGCTTCTTCATCCACGCGAATAAAATTGAGATCCGGATCGACGGCGCTCATCGCTTTTTCACAGGCATCGAGGATATCCGGGCGATATTTTTTCAGTTCTTCGAGATAGCGT  >  NC_000913/2123365‑2123888

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑