Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,009,732 C→G 6.7% 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 6.7% 291.0 / 9.2 86intergenic (+20/‑89)yedF/yedKputative TusA family sulfurtransferase/DUF159 family protein
Reads supporting (aligned to +/- strand):  ref base C (42/38);  new base G (2/4);  total (44/42)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTCGCGG  >  NC_000913/2009592‑2009878
                                                                                                                                            |                                                                                                                                                  
aggggtgAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCg                                                                                                                                            <  1:329526/146‑1 (MQ=255)
     tgAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGaa                                                                                                                                        >  7:198065/1‑148 (MQ=255)
      gAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACa                                                                                                                                      >  3:59482/1‑149 (MQ=255)
       aGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCGGGaatga                                                                                                                                     <  1:84827/149‑4 (MQ=255)
        gCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGGGTCACAGCGGGaatgag                                                                                                                                    <  5:378025/149‑5 (MQ=255)
        gCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATc                                                                                                                                    <  5:58432/149‑1 (MQ=255)
           aCTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGt                                                                                                                                  >  4:97729/1‑148 (MQ=255)
                 cGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCttt                                                                                                                           <  4:305630/149‑1 (MQ=255)
                 cGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCttt                                                                                                                           <  3:325952/149‑1 (MQ=255)
                     gTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTAt                                                                                                                       <  7:327502/149‑1 (MQ=255)
                      tCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTTGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATc                                                                                                                      <  6:275982/149‑1 (MQ=255)
                        gATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGaa                                                                                                                                        <  2:336042/129‑1 (MQ=255)
                        gATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGaa                                                                                                                                        >  1:336042/1‑129 (MQ=255)
                        gATCAATAATATTCCACTGGATGCGCGTAATCACAGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTAt                                                                                                                       >  3:374165/1‑146 (MQ=255)
                        gATCAATAATATTCCACTGGATGCGCGTAATAACGGCTATACGGTGCTGGATATTCAGCAAGACGGGTCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATca                                                                                                                     >  7:132983/1‑148 (MQ=255)
                         aTCAATAATATTCCACTGGATGCGCGTAATAACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATcaca                                                                                                                   >  6:285761/1‑149 (MQ=255)
                            aataatATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATcaca                                                                                                                   >  4:114749/1‑146 (MQ=255)
                             ataataTTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATcaca                                                                                                                   >  7:416895/1‑145 (MQ=255)
                               aataTTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGAtt                                                                                                             >  1:254476/1‑149 (MQ=255)
                               aataTTCAACTGCATGAGCGAAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCTCTGTGACCCACAGCGTGAACATCGATACATTTGTATCACATAGAtt                                                                                                             >  8:254764/1‑149 (MQ=255)
                                 tatTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGAtttt                                                                                                           <  8:198065/149‑1 (MQ=255)
                                    tCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTAcc                                                                                                        >  4:351401/1‑149 (MQ=255)
                                      cACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCg                                                                                                      >  5:128184/1‑149 (MQ=255)
                                           gATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGccccc                                                                                                 >  8:112266/1‑149 (MQ=255)
                                            aTGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTGGGTATTTATTTCAGAAGTAAGCCATCATTCCCGCGGTGCCCCACCGCTGGAACTTCGTTACCTTTGGATCACATATAGTTTACCGCCCCCCa                                                                                                >  4:75597/1‑149 (MQ=255)
                                            aTGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCa                                                                                                >  4:75596/1‑149 (MQ=255)
                                             tGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCATAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCAc                                                                                               <  6:99566/149‑1 (MQ=255)
                                              gcgcGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCCCt                                                                                              >  4:285325/1‑149 (MQ=255)
                                               cgcgTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGCCCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTc                                                                                             >  3:66642/1‑149 (MQ=255)
                                               cgcgTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTc                                                                                             <  8:161748/149‑1 (MQ=255)
                                                    aaTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTCATCCCTCATTCCCGCTGTGACCCCCAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTagt                                                                                        >  6:86704/1‑147 (MQ=255)
                                                    aaTCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTAt                                                                                         >  8:159888/1‑148 (MQ=255)
                                                     aTCACGGCTATACGGTTCTGGATATTCAGCAAGACGGGCCGACGATTCGTTATTTAATTCAGAAGTAATCCCTCATTGCCGCTGTGAGCCACAGCGGGAACATCGTTAACTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTc                                                                                       <  8:236085/149‑1 (MQ=255)
                                                      tCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCGCAGCGGGAACATCGTTAcc                                                                                                                              <  8:405976/109‑1 (MQ=255)
                                                       cACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCt                                                                                      <  4:161611/148‑1 (MQ=255)
                                                       cACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACACGCCCCCACTCCTATTCtt                                                                                     >  3:161611/1‑149 (MQ=255)
                                                        aCGGCTATACGGTGCTGGATATTCAGCAAGACGGGTCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCttt                                                                                    <  8:132983/149‑1 (MQ=255)
                                                        aCGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCttt                                                                                    >  7:138388/1‑149 (MQ=255)
                                                         cGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCtttt                                                                                   <  3:351401/149‑1 (MQ=255)
                                                         cGGATATACGGTGCGGGATATTCCCCAACACGGGCCGACCATTCTTTATTTAATTCAGAATTAATCCCTAATTCCCGCTGTGACCCACAGCGGTAACATCGTTACCTTTGTATCACATATATTTTACCCGCCCCCACTCCTATTAtttt                                                                                   <  6:128184/149‑1 (MQ=255)
                                                          ggCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTAt                                                                                                                       >  3:67283/1‑112 (MQ=255)
                                                          ggCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTAt                                                                                                                       <  4:67283/112‑1 (MQ=255)
                                                             tataAGGTGCTGGATATTCAGCAAGAAGGGCCGAACAATCGTTATTTAATTCAGAAGTAATCCCTCATTACCGCTGTGACCCACAGCGGGAACATCGTTACATTTGTATCACATAGATTTTACCCGCCCCCCCTCCTATACTATTCAg                                                                                >  3:72274/1‑148 (MQ=255)
                                                               taCGGTGCTGGATATTCAGCAAGACTGTCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACTTAGATTTTACCCGCCCCCAGTCCTATTCTTTTCAGTAt                                                                             <  5:177229/149‑1 (MQ=255)
                                                               taCGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGt                                                                               <  4:59482/147‑1 (MQ=255)
                                                               taCGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTAt                                                                             <  8:416895/149‑1 (MQ=255)
                                                                       tGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAGCATCGTGACCGTTGTAGCACATAGATTTTACCCGCCCCCGCTCCTATTCTTTTCAGTAGGATCTCag                                                                     >  6:75505/1‑149 (MQ=255)
                                                                         gATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCagg                                                                    <  7:159888/148‑1 (MQ=255)
                                                                         gATATTCAGCAAGACGGGCCGACCATTCGGTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGAATGATCTCagga                                                                   <  8:226006/149‑1 (MQ=255)
                                                                           tatTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGCTTTTACCCGCCCCCACTCCTAGTCTTTTCAGTATTATCTCaggagg                                                                 >  2:379471/1‑149 (MQ=255)
                                                                             ttCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGtt                                                               >  7:4904/1‑149 (MQ=255)
                                                                                aGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCt                                                            <  1:297564/149‑1 (MQ=255)
                                                                                  cAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCaggagg                                                                 <  8:88909/142‑1 (MQ=255)
                                                                                  cAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCaggagg                                                                 >  7:88909/1‑142 (MQ=255)
                                                                                  cAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAt                                                          >  5:123057/1‑149 (MQ=255)
                                                                                    aGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAtgt                                                        <  8:74462/149‑1 (MQ=255)
                                                                                    aGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAtgg                                                        <  5:285761/149‑2 (MQ=255)
                                                                                      aCGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAtgtgt                                                      >  5:243934/1‑149 (MQ=255)
                                                                                       cGGGCCGACCATTCTTTATTTATTTCATATGTATTCACTCATTCCCGCTGTGACCGACTGCGGGAACATCGTTCCCTTTGTATCTCATAGATTTTACACGCCCGCTCTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAtgtgtg                                                     <  4:374165/149‑1 (MQ=255)
                                                                                        gggCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTAtgtgtg                                                     >  5:174440/1‑148 (MQ=255)
                                                                                        gggCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTgg                                                    <  2:254476/149‑1 (MQ=255)
                                                                                          gCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGAc                                                  >  7:37355/1‑149 (MQ=255)
                                                                                                cATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTg                                            >  8:27874/1‑149 (MQ=255)
                                                                                                  ttCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGcc                                          <  3:97729/149‑1 (MQ=255)
                                                                                                     gTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCaa                                       <  3:114749/149‑1 (MQ=255)
                                                                                                         tttAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGAtct                                                                        <  7:124265/112‑1 (MQ=255)
                                                                                                         tttAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGAtct                                                                        >  8:124265/1‑112 (MQ=255)
                                                                                                         tttAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATccc                                   <  7:193700/149‑1 (MQ=255)
                                                                                                            aaTTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCaa                                 <  3:369694/148‑1 (MQ=255)
                                                                                                             aTTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATACTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCacac                               >  7:197273/1‑149 (MQ=255)
                                                                                                                 aGAAGTAATCCCTCATTCCCGCTGTGCCCCACAGCGGGAACATCGTTACCTTTGTATCACATAGAGTTTCCCCGCCCCCCCTCCTATGCTTTTCAGTATGATCTCAGGAGGTTTCGATGTGTGGACGCTTTGCCCAATCCCAAACGCGt                           >  7:392831/1‑149 (MQ=255)
                                                                                                                      tAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGa                      <  7:286159/149‑1 (MQ=255)
                                                                                                                        aTCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGAt                     >  4:202416/1‑148 (MQ=255)
                                                                                                                        aTCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGAt                     >  6:137539/1‑148 (MQ=255)
                                                                                                                           ccTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTAcc                 <  7:112266/149‑1 (MQ=255)
                                                                                                                           ccTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTATTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCGAATCCCAAACGCGTGAAGATTAcc                 >  6:99357/1‑149 (MQ=255)
                                                                                                                               aTTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCtt               <  8:37355/147‑1 (MQ=255)
                                                                                                                               aTTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTgc             <  3:202416/149‑1 (MQ=255)
                                                                                                                                ttCCCGCTGTGGGTCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCag                                                                     <  2:235922/92‑1 (MQ=255)
                                                                                                                                ttCCCGCTGTGGGTCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCag                                                                     >  1:235922/1‑92 (MQ=255)
                                                                                                                                ttCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTgcg            >  2:1655/1‑149 (MQ=255)
                                                                                                                                  cccGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTgcgc           <  5:137539/148‑1 (MQ=255)
                                                                                                                                   ccGCTGTGGGTCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCtt         >  7:28081/1‑149 (MQ=255)
                                                                                                                                    cGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTAGTCTTTTCAg                                                                                >  1:243648/1‑77 (MQ=255)
                                                                                                                                         tgACCAACAGCGGTAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCGCCACTCCTATTCTTTTCAGTATGATCTCAGTAGGTTTCTATGTGTGGACGCTTTGCCAAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTcgc    <  7:254764/148‑1 (MQ=255)
                                                                                                                                           aCCCACAGCGGGACCATCGTTACCTTTGTAGCACATAGCTTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGCCGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTCGCgg  >  6:117119/1‑148 (MQ=255)
                                                                                                                                            cccACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGAGCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTCGCgg  >  7:421698/1‑147 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
AGTGGTGAGCGACTGTCCGCAGTCGATCAATAATATTCCACTGGATGCGCGTAATCACGGCTATACGGTGCTGGATATTCAGCAAGACGGGCCGACCATTCGTTATTTAATTCAGAAGTAATCCCTCATTCCCGCTGTGACCCACAGCGGGAACATCGTTACCTTTGTATCACATAGATTTTACCCGCCCCCACTCCTATTCTTTTCAGTATGATCTCAGGAGGTTTCTATGTGTGGACGCTTTGCCCAATCCCAAACGCGTGAAGATTACCTTGCGCTTCTCGCGG  >  NC_000913/2009592‑2009878

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