Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,370,155 A→G intergenic (+152/‑61) pspE → / → ycjM thiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,370,1550AG99.3% 595.0 / NA 151intergenic (+152/‑61)pspE/ycjMthiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (94/56);  minor base T (1/0);  total (95/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 151 total aligned reads displayed.

AATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAACGCGTCT  >  NC_000913/1370007‑1370301
                                                                                                                                                     |                                                                                                                                                  
aaTTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCa                                                                                                                                                  <  3:525484/151‑1 (MQ=255)
    cACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGtg                                                                                                                                               <  4:338219/151‑2 (MQ=255)
    cACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGttt                                                                                                                                              >  4:816149/1‑151 (MQ=255)
    cACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGtt                                                                                                                                               <  1:966764/150‑1 (MQ=255)
       cATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTAtt                                                                                                                                           >  4:1036514/1‑151 (MQ=255)
          gAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTAtttt                                                                                                                                         <  2:27254/150‑1 (MQ=255)
          gAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTAtttt                                                                                                                                         <  1:815794/150‑1 (MQ=255)
           aGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTAtttt                                                                                                                                         >  1:1799018/1‑150 (MQ=255)
           aGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCt                                                                                                                                       <  4:749398/151‑1 (MQ=255)
           aGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCt                                                                                                                                       >  1:180429/1‑151 (MQ=255)
              ggcgCTTTTGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTg                                                                                                                                      <  4:1219747/149‑1 (MQ=255)
              ggcgCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCg                                                                                                                                    >  2:275273/1‑151 (MQ=255)
               gcgcTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCgg                                                                                                                                   >  1:786653/1‑151 (MQ=255)
                cgcTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCgg                                                                                                                                   >  1:17180/1‑150 (MQ=255)
                    tATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAATCCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGagta                                                                                                                               <  1:853082/150‑1 (MQ=255)
                         gcCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGtctct                                                                                                                         >  1:634192/1‑151 (MQ=255)
                          ccGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGtctctc                                                                                                                        >  1:700478/1‑151 (MQ=255)
                                  tttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCa                                                                                                                  <  4:1566724/150‑1 (MQ=255)
                                   ttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATg                                                                                                                >  4:1525611/1‑151 (MQ=255)
                                     ttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATggg                                                                                                              >  4:74539/1‑151 (MQ=255)
                                        tGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGAcc                                                                                                           >  1:1743324/1‑151 (MQ=255)
                                         gTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGa                                                                                                             <  1:492179/148‑1 (MQ=255)
                                         gTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGAcc                                                                                                           <  4:1047852/150‑1 (MQ=255)
                                          tCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGc                                                                                                         >  4:1721301/1‑151 (MQ=255)
                                             cacCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTAc                                                                                                      >  1:1571984/1‑151 (MQ=255)
                                                cTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCAc                                                                                                   >  1:97231/1‑151 (MQ=255)
                                                           ccccATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCaa                                                                                         <  1:1066200/150‑1 (MQ=255)
                                                           ccccATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAAc                                                                                        >  2:530538/1‑151 (MQ=255)
                                                            cccATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAAc                                                                                        >  4:1505255/1‑150 (MQ=255)
                                                            cccATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAAc                                                                                        >  4:1670966/1‑150 (MQ=255)
                                                            cccATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACa                                                                                       <  2:21031/151‑1 (MQ=255)
                                                                tAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGa                                                                                    >  2:58505/1‑150 (MQ=255)
                                                                tAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGa                                                                                    >  4:878031/1‑150 (MQ=255)
                                                                tAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGa                                                                                    >  2:226918/1‑150 (MQ=255)
                                                                 aGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGaaa                                                                                  <  4:1617354/151‑1 (MQ=255)
                                                                   attattTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAAc                                                                                 >  4:21597/1‑150 (MQ=255)
                                                                    ttattTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAAc                                                                                 >  1:57296/1‑149 (MQ=255)
                                                                     tattTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAAc                                                                                 >  4:101676/1‑148 (MQ=255)
                                                                      attTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGaa                                                                             <  4:995576/151‑1 (MQ=255)
                                                                      attTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGa                                                                              >  2:436759/1‑150 (MQ=255)
                                                                       tttGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGaa                                                                             >  3:399810/1‑150 (MQ=255)
                                                                        ttGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGaaa                                                                            >  1:1732241/1‑150 (MQ=255)
                                                                        ttGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGaaa                                                                            >  2:431976/1‑150 (MQ=255)
                                                                        ttGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGaaa                                                                            >  1:929630/1‑150 (MQ=255)
                                                                         tGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGaaaa                                                                           <  1:1063967/150‑1 (MQ=255)
                                                                          gCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAAtt                                                                         <  1:649644/151‑1 (MQ=255)
                                                                           cGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAAtt                                                                         >  4:1007850/1‑150 (MQ=255)
                                                                             tCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACg                                                                      >  4:1461573/1‑151 (MQ=255)
                                                                              cAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACgg                                                                     >  1:312745/1‑151 (MQ=255)
                                                                              cAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACgg                                                                     <  4:1338951/151‑1 (MQ=255)
                                                                              cAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACg                                                                      <  3:485340/150‑1 (MQ=255)
                                                                               aGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACgg                                                                     >  4:540281/1‑150 (MQ=255)
                                                                                gCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAATCCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGAt                                                                   <  4:1862532/151‑1 (MQ=255)
                                                                                 cTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGAtt                                                                  >  4:1766961/1‑151 (MQ=255)
                                                                                   cacaAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTAc                                                                >  1:1282062/1‑151 (MQ=255)
                                                                                     caAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCt                                                              >  3:75383/1‑151 (MQ=255)
                                                                                     caAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCt                                                              >  1:123060/1‑151 (MQ=255)
                                                                                      aaaTACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTg                                                             >  1:768876/1‑151 (MQ=255)
                                                                                       aaTACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTgg                                                            >  1:1198870/1‑151 (MQ=255)
                                                                                        aTACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGa                                                           <  4:1132731/151‑1 (MQ=255)
                                                                                          aCGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACg                                                         <  2:138874/151‑1 (MQ=255)
                                                                                            gCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGa                                                        >  4:701847/1‑150 (MQ=255)
                                                                                                tttCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCt                                                   >  4:638561/1‑151 (MQ=255)
                                                                                                  tCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTAc                                                 >  4:684731/1‑151 (MQ=255)
                                                                                                  tCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTAc                                                 >  4:35849/1‑151 (MQ=255)
                                                                                                   cccTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTAc                                                 <  4:1141707/150‑1 (MQ=255)
                                                                                                   cccTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACg                                                >  4:395180/1‑151 (MQ=255)
                                                                                                    ccTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACgg                                               >  1:1221363/1‑151 (MQ=255)
                                                                                                    ccTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACgg                                               >  4:152971/1‑151 (MQ=255)
                                                                                                     cTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACgg                                               >  1:41339/1‑150 (MQ=255)
                                                                                                     cTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGt                                              >  4:654907/1‑151 (MQ=255)
                                                                                                      tGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTg                                             >  1:1264435/1‑151 (MQ=255)
                                                                                                       ggTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAGCAGAAAATTACGGATTACCTGGACGAAATCTACGGTg                                             <  4:438481/150‑1 (MQ=255)
                                                                                                            aaaaTGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACAt                                       >  3:492844/1‑151 (MQ=255)
                                                                                                             aaaTGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACAtt                                      >  2:361149/1‑151 (MQ=255)
                                                                                                                 gATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTAcc                                  >  1:1135225/1‑151 (MQ=255)
                                                                                                                  aTTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCg                                 >  4:883728/1‑151 (MQ=255)
                                                                                                                    ttCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGc                                >  4:1494791/1‑150 (MQ=255)
                                                                                                                    ttCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGc                                >  4:817099/1‑150 (MQ=255)
                                                                                                                    ttCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGc                                >  4:120439/1‑150 (MQ=255)
                                                                                                                      ccTGCGTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAAc                             >  4:1703965/1‑151 (MQ=255)
                                                                                                                           gTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCAtt                        <  4:49588/151‑1 (MQ=255)
                                                                                                                           gTGACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCAt                         <  1:1841945/150‑1 (MQ=255)
                                                                                                                             gACTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCAttt                       >  2:327934/1‑150 (MQ=255)
                                                                                                                               cTAAAACCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTAc                     >  1:963096/1‑150 (MQ=255)
                                                                                                                                 aaaaCCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAg                   <  1:1367689/150‑1 (MQ=255)
                                                                                                                                  aaaCCCTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGaa                 <  4:144440/151‑1 (MQ=255)
                                                                                                                                     cccTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACt              >  4:503029/1‑151 (MQ=255)
                                                                                                                                      ccTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACtt             <  3:101585/151‑1 (MQ=255)
                                                                                                                                       cTTGTGCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTg            >  2:140652/1‑151 (MQ=255)
                                                                                                                                         tgtgCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAATCTTGt           <  2:18210/150‑1 (MQ=255)
                                                                                                                                          gtgCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTaa         <  1:540920/151‑1 (MQ=255)
                                                                                                                                          gtgCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTa          >  4:553485/1‑150 (MQ=255)
                                                                                                                                          gtgCTCAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTa          >  1:1449283/1‑150 (MQ=255)
                                                                                                                                               cAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAACGCGt    >  2:92295/1‑151 (MQ=255)
                                                                                                                                               cAATTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAACGCGt    <  4:1330570/151‑1 (MQ=255)
                                                                                                                                                 aTTGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAACGCGTCt  <  1:70235/151‑1 (MQ=255)
                                                                                                                                                  ttGTCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAACGCGACt  >  3:371886/1‑150 (MQ=255)
                                                                                                                                                  ttGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAACGCGt    >  3:173787/1‑148 (MQ=255)
                                                                                                                                                  ttGGCAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAACGCGt    <  2:173794/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGTCTCTCGTTTCATGGGACCGCTACCACGGAAAGGCAACATGAAACAGAAAATTACGGATTACCTGGACGAAATCTACGGTGGAACATTTACCGCAACTCATTTACAGAAACTTGTAACGCGTCT  >  NC_000913/1370007‑1370301

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: