Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,123,634 C→G 5.7% G296A (GGC→GCC) ‡ cpsB ← mannose‑1‑phosphate guanyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,123,6340CG 5.7% 253.1 / 3.1 77G296A (GGC→GCC) ‡cpsBmannose‑1‑phosphate guanyltransferase
Reads supporting (aligned to +/- strand):  ref base C (29/42);  new base G (2/3);  total (32/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.55e-01

GGTAGCGGTCGCCCGCGTCGATAGAGTCATATTTGCCCCACGGACGATACACTTCGCGATGCACCCGATGCTCATGGCGACCATCGGCTTTGATCTGCTCGACCACTTTTTTCACATCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACGCTTCTTCATCCACGCGAATAAAATTGAGATCCGGATCGACGGCGCTCATCGCTTTTTCACAGGCATCGAGGATATCCGGGCGATATTTTTTCAGTTCTTCGAGATAGCGT  >  NC_000913/2123365‑2123888
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                              
ggTAGCGGTCGCCCGCGTCGATAGAGTCATATTTGCCCCACGGACGATACACTTCGCGATGCACCCGATGCTCATGGCGACCATCGGCTTTGATCTGCTCGACCACTTTTTTCACATCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCCTACACATAGCTGTTTTCAGTTTTGTGATTAATCCCATCGCCGTGGCAAACGTTGccc                                                                                                                                                                                                                                                              >  2:371387/1‑272 (MQ=255)
                                                     tCGCGATGCACCCGATGCTCATGGCGACCATCGGCTTTGATCTGCTCGACCACTTTTTTCACATCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTTGGCGCTCATTCCCCATAATGAAGACCAGGAGCCAACATCGCTCCa                                                                                                                                                                                                      >  2:169336/1‑275 (MQ=255)
                                                     tCGCGATGCACCCGATGCTCATGGCGACCATCGGCTTTGATCTGCTCGACCACTTTTTTCACATCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCa                                                                                                                                                                                                      >  2:224649/1‑275 (MQ=255)
                                                           tGCACCCGATGCTCATGGCGACCATCGGCTTTGATCTGCTCGACCACTTTTTTCACATCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCTCTCCAGCCCGc                                                                                                                                                                                                >  1:19665/1‑275 (MQ=255)
                                                                     gCTCATGGCGACCATCGGCTTTGATCTGCTCGACCACTTTTTTCACATCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTTGGCGGTGTGGGCGCTGACCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCgg                                                                                                                                                                                       <  2:99811/274‑1 (MQ=255)
                                                                              gACCATCGGCTTTGATCTGCTCGACCACTTTTTTCACATCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCacaaca                                                                                                                                                                             >  2:499765/1‑275 (MQ=255)
                                                                                           gttCTGCTCGACCCCTTTTTTGCCATCCTGTACCGCTTTAGGGTTGGAAATTAGCCCCGCATCTTGTGGCCTCACCACTTCCAGATCTTTCCCCCCGGCGGGGGTGACCAGGGCGGATTCAGCATACACATAGCTGTTTTCAGGTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGTCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGGGCCACCATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTa                                                                                                                                                                 <  4:499760/272‑1 (MQ=255)
                                                                                                 cTCGACCACTTTTTTCACATCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCCTAATGAAGACCAGGGGCCAACATCGCTCCAGCCCGCATCCCTCGGCCACACCACAGCATCTgtcccttaaaca                                                                                                                                                           >  1:135454/1‑263 (MQ=255)
                                                                                                   gcccctcTTTTTTCCCTCCGTTACCCCTTTACGTGCGGCAATCAGACCCGCTACTTTGGTCTGCACCACTACCACATCTTTCACCCCGACGGTGATGACCAGGCCAGATTCAGCATACACATAGCTGTGTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGa                                                                                                                                                       <  4:224645/269‑1 (MQ=255)
                                                                                                     gccACTTTTTTCACATCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTa                                                                                                                                                                 <  2:21405/263‑1 (MQ=255)
                                                                                                                    tCCTTTACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCGGATCTTTCACCCCGGCGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGa                                                                                                                                       <  2:401240/275‑1 (MQ=255)
                                                                                                                    tCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGATGATCTCCCATAATGAAGACCAGGAGCCAACATCGCCCCAGCCCGCATCCCACCGCACCACAACAGCATCTGCCTTAACGTCCATGACCGCGGAATCccccaa                                                                                                                                       >  4:390889/1‑273 (MQ=255)
                                                                                                                             gcgTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAg                                                                                                                                                                                                                            >  4:52797/1‑181 (MQ=255)
                                                                                                                             gcgTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCg                                                                                                                               >  2:347092/1‑274 (MQ=255)
                                                                                                                              cgTTACGGTCGGCAACCAGAACGCCATCTTTGGCGTGACCCCACACCAGTTCTTTCACCCCGACGTGGGTGACCAGGCCGGATTCGTCATACCAATAGCTGGTTTCAGTTGTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCGCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCggg                                                                                                                             <  4:184348/275‑1 (MQ=255)
                                                                                                                                ttACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGACGTGGCAAACGTTGACCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCACTCCAGCCCGCATCCATCGGCCCCCCCACCGCATCTGCGCGCCGTTCCTTTACCCCGAGACCCCACAACCCTTCCGGGCa                                                                                                                           >  4:495171/1‑275 (MQ=255)
                                                                                                                                 tACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCa                                                                                                                                                                                                                        <  2:52798/181‑1 (MQ=255)
                                                                                                                                 tACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACa                                                                                                                                                                                                              <  2:228823/191‑1 (MQ=255)
                                                                                                                                        ggCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCgaga                                                                                                                    >  2:128781/1‑274 (MQ=255)
                                                                                                                                        ggCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGaa                                                                                                                   <  2:390894/275‑1 (MQ=255)
                                                                                                                                            aTCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCt                                                                                                                                                                                                                                             <  8:118159/149‑1 (MQ=255)
                                                                                                                                             tCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTg                                                                                                                                                                                                                                            >  7:290387/1‑149 (MQ=255)
                                                                                                                                                gCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGAtc                                                                                                                                                                                                                                         <  5:114979/149‑1 (MQ=255)
                                                                                                                                                  aCCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGAtctc                                                                                                                                                                                                                                       <  5:403465/149‑1 (MQ=255)
                                                                                                                                                      cATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCAt                                                                                                                                                                                                                                   <  8:290387/149‑1 (MQ=255)
                                                                                                                                                         cTTTGGTCTGCACCACTACCACATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACGCTTCTTCATCCAc                                                                                                   <  2:495176/274‑1 (MQ=255)
                                                                                                                                                             tGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGaa                                                                                                                                                                                                                             <  5:295161/149‑1 (MQ=255)
                                                                                                                                                             tGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGaa                                                                                                                                                                                                                             <  7:97429/149‑1 (MQ=255)
                                                                                                                                                             tGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGaa                                                                                                                                                                                                                             <  8:93352/149‑1 (MQ=255)
                                                                                                                                                             tGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAAAGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGaa                                                                                                                                                                                                                             <  8:276617/149‑1 (MQ=255)
                                                                                                                                                                tCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGAc                                                                                                                                                                                                                          <  6:26262/149‑1 (MQ=255)
                                                                                                                                                                  tGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGc                                                                                                                                                                     >  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                                                                                                                                                     <  7:453187/138‑1 (MQ=255)
                                                                                                                                                                                                                                                 tgATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGAcc                                                                                                                                                     >  8:453187/1‑138 (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)
                                                                                                                                                                                                                                                                              gCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACGctt                                                                                                             <  6:370387/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                              
GGTAGCGGTCGCCCGCGTCGATAGAGTCATATTTGCCCCACGGACGATACACTTCGCGATGCACCCGATGCTCATGGCGACCATCGGCTTTGATCTGCTCGACCACTTTTTTCACATCCTGTACCGCGTTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACGCTTCTTCATCCACGCGAATAAAATTGAGATCCGGATCGACGGCGCTCATCGCTTTTTCACAGGCATCGAGGATATCCGGGCGATATTTTTTCAGTTCTTCGAGATAGCGT  >  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: