Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,727,780 G→T 5.3% intergenic (+24/‑57) nemA → / → gloA chromate reductase, quinone reductase, FMN‑linked; N‑Ethylmaleimide reductase; old yellow enzyme/glyoxalase I, Ni‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,727,7800GT 5.3% 317.9 / 5.3 94intergenic (+24/‑57)nemA/gloAchromate reductase, quinone reductase, FMN‑linked; N‑Ethylmaleimide reductase; old yellow enzyme/glyoxalase I, Ni‑dependent
Reads supporting (aligned to +/- strand):  ref base G (52/37);  new base T (2/3);  total (54/40)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

GATTGCGAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCAAAGTGCTGGGCATGAAACTGCTGCGTAC  >  NC_000913/1727639‑1727928
                                                                                                                                             |                                                                                                                                                    
gATTGCGAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAgccgcc                                                                                                                                                 <  3:399358/147‑1 (MQ=255)
  ttGCGAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCCCAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCta                                                                                                                                             >  5:261651/1‑149 (MQ=255)
   tcatgcccagcactttggtataaaaatcgatggagcgttgcaaatcgccaacgcgcagcatggtatgaagaagacgcatctttttatcctcaaaatggctaacagattcaaaatgttgttttagtatAGCGGCGTAAAGCCGCCGCta                                                                                                                                             <  2:119966‑M1/21‑1 (MQ=255)
     cGAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCta                                                                                                                                             <  3:278266/146‑1 (MQ=255)
     cGAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCta                                                                                                                                             >  4:278266/1‑146 (MQ=255)
     cGAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATAc                                                                                                                                          >  3:435108/1‑149 (MQ=255)
     cGAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATAc                                                                                                                                          >  7:392669/1‑149 (MQ=255)
     atgcccagcactttggtataaaaatcgatggagcgttgcaaatcgccaacgcgcagcatggtatgaagaagacgcatctttttatcctcaaaatggctaacagattcaaaatgttgttttagtatAGCGGCGTAAAGCCGCCGCTATAc                                                                                                                                          >  1:119966‑M1/126‑149 (MQ=255)
      gAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATAc                                                                                                                                          >  4:60434/1‑148 (MQ=255)
           cGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTaaaa                                                                                                                                     >  7:98298/1‑148 (MQ=255)
             gATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCTTGCAGAAAGATTCTACGGTGGAGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATACAACATTGCGAGCGGAGTACAGCCGCCGCTATACAAAaaca                                                                                                                                   >  8:301377/1‑148 (MQ=255)
                  ggTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCACCGCTATACTAAAACATCAttt                                                                                                                             >  3:454531/1‑149 (MQ=255)
                   gTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGTCGTAAAGCCACCGCTATACTAAAACAACAtttt                                                                                                                            <  4:454531/149‑1 (MQ=255)
                   gTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACAtttt                                                                                                                            <  6:98970/149‑1 (MQ=255)
                   gTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGGGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAacaaca                                                                                                                                <  6:261651/145‑1 (MQ=255)
                     cGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGa                                                                                                                          <  8:392669/149‑1 (MQ=255)
                                  cgcAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCAt                                                                                                             <  3:169552/149‑1 (MQ=255)
                                   gcAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCAtt                                                                                                            >  7:308638/1‑149 (MQ=255)
                                       gttgcaaatcgccaacgcgcagcatggtatgaagaagacgcatctttttatcctcaaaatggctaacagattcaaaatgttgttttagtatAGCGGCGGCTTTACGCCGCTATACTAAAACAACATTTTGAATCTGTTa                                                                                                                  >  2:152295‑M1/92‑139 (MQ=255)
                                       gttgcaaatcgccaacgcgcagcatggtatgaagaagacgcatctttttatcctcaaaatggctaacagattcaaaatgttgttttagtatAGCGGCGGCTTTACGCCGCTATACTAAAACAACATTTTGAATCTGTTa                                                                                                                  >  2:158942‑M1/92‑139 (MQ=255)
                                       gttgcaaatcgccaacgcgcagcatggtatgaagaagacgcatctttttatcctcaaaatggctaacagattcaaaatgttgttttagtatAGCGGCGGCTTTACGCCGCTATACTAAAACAACATTTTGAATCTGTTa                                                                                                                  <  1:158942‑M1/48‑1 (MQ=255)
                                       gttgcaaatcgccaacgcgcagcatggtatgaagaagacgcatctttttatcctcaaaatggctaacagattcaaaatgttgttttagtatAGCGGCGGCTTTACGCCGCTATACTAAAACAACATTTTGAATCTGTTa                                                                                                                  <  1:152295‑M1/48‑1 (MQ=255)
                                        gCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTCGCCATTTTGAg                                                                                                       >  6:320501/1‑149 (MQ=255)
                                         cTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAgg                                                                                                      <  6:486614/149‑1 (MQ=255)
                                          tGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGa                                                                                                     >  7:372945/1‑149 (MQ=255)
                                           gAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGa                                                                                                     >  5:91396/1‑148 (MQ=255)
                                            aGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATa                                                                                                   <  2:328840/149‑1 (MQ=255)
                                                   cccACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGAt                                                                                            <  4:389078/149‑1 (MQ=255)
                                                   cccACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGAt                                                                                            <  5:440099/149‑1 (MQ=255)
                                                   cccACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGATGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGAt                                                                                            <  6:290856/149‑1 (MQ=255)
                                                    ccACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACAttt                                                                                                                             >  6:306201/1‑115 (MQ=255)
                                                    ccACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACAttt                                                                                                                             <  5:306201/115‑1 (MQ=255)
                                                    ccACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATg                                                                                           >  3:386213/1‑149 (MQ=255)
                                                      acaGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAacaaca                                                                                                                                <  6:35218/110‑1 (MQ=255)
                                                      acaGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAacaaca                                                                                                                                >  5:35218/1‑110 (MQ=255)
                                                       caGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGt                                                                                        >  4:328474/1‑149 (MQ=255)
                                                        aGCGGGCCGAACGTTTCTACGGTGGCGGCGCGGCAGGCTATCCCGATTACCCGACGGTGTCATCCAACCTTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGtc                                                                                       <  4:381839/149‑1 (MQ=255)
                                                           gagCCGAAAGTTTCTACGCTGGCGGCGCGGAAGGCTAAACCCATTACCCGACGTTGTAAACCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGtcta                                                                                     >  4:187394/3‑147 (MQ=255)
                                                               cGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAg                                                                                                       <  1:454294/126‑1 (MQ=255)
                                                               cGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAg                                                                                                       >  2:454294/1‑126 (MQ=255)
                                                                gAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAacaaca                                                                                                                                >  1:442253/1‑100 (MQ=255)
                                                                gAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAacaaca                                                                                                                                <  2:442253/100‑1 (MQ=255)
                                                                gAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCa                                                                                >  6:39475/1‑148 (MQ=255)
                                                                gAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCa                                                                                >  7:451879/1‑148 (MQ=255)
                                                                gAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCa                                                                                >  4:418445/1‑148 (MQ=255)
                                                                gAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCa                                                                                >  8:344764/1‑148 (MQ=255)
                                                                gAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCa                                                                                >  1:14944/1‑148 (MQ=255)
                                                                gAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCAt                                                                               <  3:60434/149‑1 (MQ=255)
                                                                 aaaGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGAt                                                                                            <  6:456422/135‑1 (MQ=255)
                                                                 aaaGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGAt                                                                                            >  5:456422/1‑135 (MQ=255)
                                                                 aaaGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATa                                                                              >  4:14706/1‑149 (MQ=255)
                                                                        cTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCAtgct                                                                       >  2:313118/1‑149 (MQ=255)
                                                                            ggTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTgcgc                                                                   <  8:451879/149‑1 (MQ=255)
                                                                             gTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCCACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATTCTGcgcg                                                                  >  8:447697/1‑149 (MQ=255)
                                                                                gcggcgCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTg                                                               >  5:209685/1‑149 (MQ=255)
                                                                                  ggcgCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGc                                                             <  3:14706/149‑1 (MQ=255)
                                                                                     gcgGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGAt                                                          <  8:308638/149‑1 (MQ=255)
                                                                                       ggAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGAttt                                                        <  3:328474/149‑1 (MQ=255)
                                                                                        gAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTg                                                       <  4:386213/149‑1 (MQ=255)
                                                                                         aaGGCTATACCGATTACCCGACGTTGTAATCCAACATTGTGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGAt                                                          >  1:340752/1‑145 (MQ=255)
                                                                                         aaGGCTATACCGATTACCCGACGTTGTAATCCAACATTGTGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGAt                                                          <  2:340752/145‑1 (MQ=255)
                                                                                              tataCCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGc                                                 >  2:444140/1‑149 (MQ=255)
                                                                                              tataCCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGACTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGc                                                 >  4:80899/1‑149 (MQ=255)
                                                                                                 aCCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTc                                               >  4:110911/1‑148 (MQ=255)
                                                                                                  ccGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTcc                                              <  2:14944/148‑1 (MQ=255)
                                                                                                     aTTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGc                                                                                          >  8:28313/1‑101 (MQ=255)
                                                                                                     aTTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGc                                                                                          <  7:28313/101‑1 (MQ=255)
                                                                                                      ttACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGATCCATCGa                                         >  7:500622/1‑149 (MQ=255)
                                                                                                         cccGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAATCAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTAATACCCTCCTGCGCGTTGGCGAGTTGCAACGCTCAATCGAttt                                      >  6:364694/1‑149 (MQ=255)
                                                                                                         cccGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACCACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATTCCATGCTGCGCGTTGGCGATTTGAAACGCTCCATCGAttt                                      >  8:223910/1‑149 (MQ=255)
                                                                                                         cccGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGAttt                                      <  3:418445/149‑1 (MQ=255)
                                                                                                         cccGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGAttt                                      >  5:268007/1‑149 (MQ=255)
                                                                                                            gACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTta                                   >  5:237290/1‑149 (MQ=255)
                                                                                                            gACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTta                                   >  1:86433/1‑149 (MQ=255)
                                                                                                            attcaaaatgttgttttagtatAGCGGCGGCTTTACGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTta                                   <  8:59934‑M1/127‑1 (MQ=255)
                                                                                                             aCGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTtat                                  <  6:237290/149‑1 (MQ=255)
                                                                                                              cggtGTAATCCACCATTGCGAGCCGCGGAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTACCCATTTTGAGGATAAAAAGCTGCGTCTTCTTCATACCAAGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTtata                                 <  6:15142/146‑1 (MQ=255)
                                                                                                               gTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATAc                                >  1:371256/1‑149 (MQ=255)
                                                                                                                ttGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATAc                                >  4:267829/1‑148 (MQ=255)
                                                                                                                   tAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCaaa                            <  1:313118/149‑1 (MQ=255)
                                                                                                                    aaTCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCAAAg                           >  8:102521/1‑149 (MQ=255)
                                                                                                                     aTCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCAAAGt                          >  3:295733/1‑149 (MQ=255)
                                                                                                                     aTCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCAAAGt                          <  8:372945/149‑1 (MQ=255)
                                                                                                                      tCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCAAAGTg                         >  4:146622/1‑149 (MQ=255)
                                                                                                                        cAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCAAAGTGCt                       <  3:187394/149‑1 (MQ=255)
                                                                                                                          aCATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCAAAGTGCTgg                     <  3:110911/149‑1 (MQ=255)
                                                                                                                           cATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCAAAGTGCTggg                    >  1:244840/1‑149 (MQ=255)
                                                                                                                             ttGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTta                                   <  2:179941/132‑1 (MQ=255)
                                                                                                                             ttGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTta                                   >  1:179941/1‑132 (MQ=255)
                                                                                                                                 gAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCAAAGTGCTGGGCATGa               >  5:35872/1‑148 (MQ=255)
                                                                                                                                   GCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCAAAGTGCTGGGCATGAAAc            >  7:102521‑M1/1‑149 (MQ=255)
                                                                                                                                       CGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCAAAGTGCTGGGCATGAAActgct        <  1:484772‑M1/149‑1 (MQ=255)
                                                                                                                                        GTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCAAAGTGCTGGGCATGAAActgctg       >  5:39475‑M1/1‑149 (MQ=255)
                                                                                                                                        GTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCAAAGTGCTGGGCATGAAActgctg       <  2:424798‑M1/149‑1 (MQ=255)
                                                                                                                                          aaaGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCAAAGTGCTGGGCATGAAACTGCTGCg     >  8:156880/1‑149 (MQ=255)
                                                                                                                                             gccgccGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCAAAGTGCTGGGCATGAAACTGCTGCGTAc  <  3:146622/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
GATTGCGAACCCGGATCTGGTCGCCCGCTTGCAGCGCAAAGCTGAGCTTAACCCACAGCGTGCCGAAAGTTTCTACGGTGGCGGCGCGGAAGGCTATACCGATTACCCGACGTTGTAATCCAACATTGCGAGCGGCGTAAAGCCGCCGCTATACTAAAACAACATTTTGAATCTGTTAGCCATTTTGAGGATAAAAAGATGCGTCTTCTTCATACCATGCTGCGCGTTGGCGATTTGCAACGCTCCATCGATTTTTATACCAAAGTGCTGGGCATGAAACTGCTGCGTAC  >  NC_000913/1727639‑1727928

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