Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,009,732 C→G 5.3% intergenic (+20/‑89) yedF → / → yedK putative TusA family sulfurtransferase/DUF159 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,009,7320CG 5.3% 264.7 / 6.0 76intergenic (+20/‑89)yedF/yedKputative TusA family sulfurtransferase/DUF159 family protein
Reads supporting (aligned to +/- strand):  ref base C (40/32);  new base G (2/2);  total (42/34)
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.93e-01

TCCTGGAAGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTCGCGGA  >  NC_000913/2009585‑2009879
                                                                                                                                                   |                                                                                                                                                   
tCCTGGAAGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGAc                                                                                                                                                     <  6:115028/148‑1 (MQ=255)
           gtgAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGa                                                                                                                                          <  2:362425/148‑1 (MQ=255)
             gAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGaa                                                                                                                                         >  3:348711/1‑147 (MQ=255)
                   cTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGtt                                                                                                                                  <  7:123420/148‑1 (MQ=255)
                       ccGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCt                                                                                                                              >  1:169347/1‑148 (MQ=255)
                        cGCAGTCGATCAATAATATTCCACTGGATGCTCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCtt                                                                                                                             <  8:185878/148‑1 (MQ=255)
                        cGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCtt                                                                                                                             <  8:162530/148‑1 (MQ=255)
                          cAGTCGATCAATAATATTCCACTGGATGCGCTTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTg                                                                                                                           <  1:62795/148‑1 (MQ=255)
                              cGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTAt                                                                                                                        >  3:80460/1‑147 (MQ=255)
                                aTCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGGTGTGACCCACAGCGGGACCATCGTTTGATTTGTATca                                                                                                                      >  2:126645/1‑147 (MQ=255)
                                aTCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTCGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCCTTACCTTTGTATca                                                                                                                      >  8:216709/1‑147 (MQ=255)
                                      aataTTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGAt                                                                                                               >  6:45875/1‑148 (MQ=255)
                                                   aTGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGt                                                                                                                          >  2:193294/1‑124 (MQ=255)
                                                   aTGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGt                                                                                                                          <  1:193294/124‑1 (MQ=255)
                                                   aTGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGccccc                                                                                                  <  2:169347/148‑1 (MQ=255)
                                                   aTGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGccccc                                                                                                  <  8:32499/148‑1 (MQ=255)
                                                     gcgcGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCAc                                                                                                >  2:363475/1‑148 (MQ=255)
                                                      cgcgTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACt                                                                                               <  6:288900/148‑1 (MQ=255)
                                                             tCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTc                                                                                        >  7:371640/1‑148 (MQ=255)
                                                             tCACGGCGTTACGGTGCTGGATATTCAGCAAGAGGGGCCGACCATTCGTTATTTAATTCAGAAGTACTCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTAGTATCACATAGAGTTTACCCGCCCCCa                                                                                                 <  6:360104/139‑1 (MQ=255)
                                                              cACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAAATCAGAAGTAATCACTAATTCCCGCTGTGACCCACATCGGGAAAATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACtagatagg                                                                                        >  5:360104/1‑140 (MQ=255)
                                                                gggCTTTACGGTGCTTGATATTCAGGATGACGGGCCGTCGATTCGTTATTTAATTCAGAAGTAATCCCTCATTCGCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATGACATAGATTTTACCCGCCCCCACTCCTATTCttt                                                                                     <  5:377554/147‑1 (MQ=255)
                                                                cGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAAATCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTATGTATCACATACATTTAACCCGCCCCCACACCTATACttt                                                                                     >  6:377554/1‑148 (MQ=255)
                                                                  gCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCCTTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTc                                                                                   <  4:294061/148‑1 (MQ=255)
                                                                  gCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCtttt                                                                                    >  7:240749/1‑147 (MQ=255)
                                                                  gCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCtttt                                                                                    >  5:74736/1‑147 (MQ=255)
                                                                      taCGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGt                                                                                <  4:80460/147‑1 (MQ=255)
                                                                       aCGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTAt                                                                              >  6:119750/1‑148 (MQ=255)
                                                                         ggTGCTGGATATTCAGCAAGACGGTCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGa                                                                            >  3:173650/1‑148 (MQ=255)
                                                                            gCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGAtct                                                                         >  3:276220/1‑148 (MQ=255)
                                                                              tGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCCCTCCTATTCTTTTCAGTATGAtct                                                                         >  3:237732/1‑146 (MQ=255)
                                                                                gATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCag                                                                      >  4:81238/1‑147 (MQ=255)
                                                                                 atatTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCag                                                                      >  3:316221/1‑146 (MQ=255)
                                                                                  tatTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTACCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCaggag                                                                   >  8:359205/1‑148 (MQ=255)
                                                                                        gCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACTTCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATATCAGTAGGTTTc                                                              >  5:103659/1‑147 (MQ=255)
                                                                                        gCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCt                                                             <  8:371640/148‑1 (MQ=255)
                                                                                         cAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTa                                                            <  4:316221/148‑1 (MQ=255)
                                                                                          aaGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAt                                                           >  3:107241/1‑148 (MQ=255)
                                                                                             aCGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAtgtg                                                        <  8:240749/148‑1 (MQ=255)
                                                                                              cGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAtgtgt                                                       <  5:45875/148‑1 (MQ=255)
                                                                                                ggCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTgg                                                     <  5:119750/148‑1 (MQ=255)
                                                                                                   cGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTgg                                                     >  1:34362/1‑145 (MQ=255)
                                                                                                   cGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTgg                                                     <  2:34362/145‑1 (MQ=255)
                                                                                                       cATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCttt                                              <  4:348711/148‑1 (MQ=255)
                                                                                                          tCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGtt                                                                <  7:43133/127‑1 (MQ=255)
                                                                                                          tCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGtt                                                                >  8:43133/1‑127 (MQ=255)
                                                                                                          tCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGcc                                           >  7:383021/1‑148 (MQ=255)
                                                                                                          tCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGcc                                           >  8:33336/1‑148 (MQ=255)
                                                                                                          tCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGcc                                           >  2:230451/1‑148 (MQ=255)
                                                                                                                tttAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATcc                                     >  1:145761/1‑148 (MQ=255)
                                                                                                                 ttAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATcc                                     >  8:3754/1‑147 (MQ=255)
                                                                                                                 ttAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATcc                                     >  8:9101/1‑147 (MQ=255)
                                                                                                                      tCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTa                                                                               <  1:193348/100‑1 (MQ=255)
                                                                                                                      tCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTa                                                                               >  2:193348/1‑100 (MQ=255)
                                                                                                                      tCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAAcg                               >  8:240091/1‑148 (MQ=255)
                                                                                                                        aGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAAcgcg                             <  6:29768/148‑1 (MQ=255)
                                                                                                                        aGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAAcgcg                             >  2:309829/1‑148 (MQ=255)
                                                                                                                           aGTAATCCCTCATTCCCGCTTTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCGAATCCCAAACGCGTGa                          <  2:40952/148‑1 (MQ=255)
                                                                                                                            gTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATATATTTTACCCGCCCCAACTCCTATTCTTTTAAGTATTATCTCAGGAGGTTTATATGTGTGGACGCTTTGCCCAATCACAAACGCGTg                           >  4:132667/1‑146 (MQ=255)
                                                                                                                                 cccTCATTCACGCTGTGACCCACAGCGGGAACAACGTAAGCTTTGTATCACATAGATTTTACCCGCCCCCAGACCTATGATTTTCAGTATAATCTCAGCAGGTATCTATGTGGGGGCCCTTTGCACAATCCCAAACGCGTGACGATTa                    >  4:262450/1‑148 (MQ=255)
                                                                                                                                  ccTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTa                    <  1:230451/147‑1 (MQ=255)
                                                                                                                                   cTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTAcc                  <  4:237732/148‑1 (MQ=255)
                                                                                                                                     tgttCCCGCTGTGGGTCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGaa                         <  5:118056/137‑1 (MQ=255)
                                                                                                                                     tgttCCCGCTGTGGGTCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGaa                         >  6:118056/3‑139 (MQ=255)
                                                                                                                                     cATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGTAGGTTTCTATGTGTGGACGCTTTGGCAAATCCCAAACGCGTGAAGATTACCtt                >  8:279217/1‑148 (MQ=255)
                                                                                                                                     cATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCtt                >  8:16918/1‑148 (MQ=255)
                                                                                                                                        tCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCaggagg                                                                  >  2:118912/1‑95 (MQ=255)
                                                                                                                                        tCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCaggagg                                                                  <  1:118912/95‑1 (MQ=255)
                                                                                                                                        tCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGAAGCTTTGCACAATCCCAAACGCGTGAAGATTACCTTgcg             >  6:166797/1‑148 (MQ=255)
                                                                                                                                         cccGCTGTGGGTCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTa                                                            >  7:39708/1‑100 (MQ=255)
                                                                                                                                         cccGCTGTGGGTCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTa                                                            <  8:39708/100‑1 (MQ=255)
                                                                                                                                         cccGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTgcg             <  7:359205/147‑1 (MQ=255)
                                                                                                                                           cGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTCCCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCtt          >  5:105685/1‑148 (MQ=255)
                                                                                                                                                tgACCCACAGAGGGAACATCGTTACCTTTGTATCACATATATTTTCCCCGCCCCCCCTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATTTGTGGACGCTTTGCCAAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTcgc     <  3:132667/148‑1 (MQ=255)
                                                                                                                                                ggACCCACAGCGGGAACATCGTTACCTTTGTCTCACATCGATTTTACCCGCCCCCACTCCTCTTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTcg      <  6:103659/146‑1 (MQ=255)
                                                                                                                                                  aCCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTcgcg    <  4:107241/147‑1 (MQ=255)
                                                                                                                                                  aCCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTCGCgg   <  1:309829/148‑1 (MQ=255)
                                                                                                                                                   cccACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTCGCGGa  <  6:74736/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCCTGGAAGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTCGCGGA  >  NC_000913/2009585‑2009879

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