Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,370,043 (T)9→7 intergenic (+40/‑172) pspE → / → ycjM thiosulfate sulfurtransferase PspE/glucosylglycerate phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,370,0360T.100.0% 431.3 / NA 102intergenic (+33/‑180)pspE/ycjMthiosulfate sulfurtransferase PspE/glucosylglycerate phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (51/51);  total (51/51)
*NC_0009131,370,0370T.99.0% 422.4 / ‑3.1 102intergenic (+34/‑179)pspE/ycjMthiosulfate sulfurtransferase PspE/glucosylglycerate phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (50/51);  minor base A (1/0);  total (51/51)
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 121 total aligned reads displayed.

GACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTA  >  NC_000913/1369897‑1370176
                                                                                                                                           ||                                                                                                                                           
gACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGt                                                                                                                                              >  1:939143/1‑140 (MQ=255)
      aGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttt                                                                                                                                        <  1:1111282/140‑1 (MQ=255)
       gTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGtttttttc                                                                                                                                       >  1:2996670/1‑139 (MQ=255)
       gTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGtttttttc                                                                                                                                       >  2:86417/1‑139 (MQ=255)
       gTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttt                                                                                                                                        <  2:397608/139‑1 (MQ=255)
       gTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttt                                                                                                                                        >  1:397608/1‑139 (MQ=255)
          aGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGtttttttctgt                                                                                                                                    >  2:2170765/1‑136 (MQ=255)
           gggCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGtttttttctgtc                                                                                                                                   <  2:3066550/140‑6 (MQ=255)
           gggCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGtttttttctgtc                                                                                                                                   <  1:2402700/140‑6 (MQ=255)
            ggCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGtttttttctgtca                                                                                                                                  <  1:1540921/140‑7 (MQ=255)
            ggCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGtttttttctgtca                                                                                                                                  <  2:1153888/140‑7 (MQ=255)
            ggCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGtttttttctgtca                                                                                                                                  >  1:1464143/1‑134 (MQ=255)
            ggCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGtttttttctgtca                                                                                                                                  <  2:1833876/140‑7 (MQ=255)
            ggCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGtttttttctgtca                                                                                                                                  <  2:119635/140‑7 (MQ=255)
            ggCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGtttttttctgtca                                                                                                                                  <  2:2996670/140‑7 (MQ=255)
             gcaagcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGtttttttctgtcac                                                                                                                                 >  2:1619241/1‑133 (MQ=255)
              caagcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTcaca                                                                                                                              <  2:1464143/140‑1 (MQ=255)
               aagcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCacac                                                                                                                             >  2:2509931/1‑140 (MQ=255)
                       gagaTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTa                                                                                                                     <  1:3039242/140‑1 (MQ=255)
                         gaTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACa                                                                                                                   <  1:2288685/140‑1 (MQ=255)
                         gaTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACa                                                                                                                   <  2:1807778/140‑1 (MQ=255)
                         gaTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACa                                                                                                                   <  2:1794607/140‑1 (MQ=255)
                         gaTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACa                                                                                                                   <  1:1206912/140‑1 (MQ=255)
                         gaTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACa                                                                                                                   <  1:3388647/140‑1 (MQ=255)
                          aTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACAt                                                                                                                  >  2:1428675/1‑140 (MQ=255)
                            ccTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTcaca                                                                                                                              <  2:3000531/126‑1 (MQ=255)
                               tAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCa                                                                                                             <  2:939143/140‑1 (MQ=255)
                               tAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCa                                                                                                             <  1:2509931/140‑1 (MQ=255)
                                aGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCAt                                                                                                            >  1:853212/1‑140 (MQ=255)
                                   gagaTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGa                                                                                                         <  2:2612452/140‑1 (MQ=255)
                                    agaTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGat                                                                                                        >  2:80780/1‑140 (MQ=255)
                                       tGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCATAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTAc                                                                                                                    >  2:2020020/1‑125 (MQ=255)
                                       tGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGattat                                                                                                     >  1:1755164/1‑140 (MQ=255)
                                             ataCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGt                                                                                               <  2:2814603/140‑1 (MQ=255)
                                             ataCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGt                                                                                               <  2:1755164/140‑1 (MQ=255)
                                             ataCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGt                                                                                               <  2:853212/140‑1 (MQ=255)
                                             ataCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGt                                                                                               <  2:3357327/140‑1 (MQ=255)
                                             ataCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGt                                                                                               <  2:270971/140‑1 (MQ=255)
                                             ataCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGt                                                                                               <  1:3129599/140‑1 (MQ=255)
                                             ataCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGt                                                                                               <  1:2995694/140‑1 (MQ=255)
                                                  cACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCt                                                                                          <  2:1443105/140‑1 (MQ=255)
                                                   aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTc                                                                                         >  1:1533854/1‑140 (MQ=255)
                                                   aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTc                                                                                         >  1:2811029/1‑140 (MQ=255)
                                                   aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTc                                                                                         >  2:571564/1‑140 (MQ=255)
                                                   aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTc                                                                                         <  1:170654/140‑1 (MQ=255)
                                                    cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTca                                                                                        >  1:3330726/1‑140 (MQ=255)
                                                    cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTca                                                                                        >  1:67182/1‑140 (MQ=255)
                                                    cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGCCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTca                                                                                        >  2:2119375/1‑140 (MQ=255)
                                                        gagaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGatta                                                                                                      <  1:2897647/122‑1 (MQ=255)
                                                        gagaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGatta                                                                                                      >  2:2897647/1‑122 (MQ=255)
                                                         agaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAAt                                                                                   >  1:2112271/1‑140 (MQ=255)
                                                          gaACGCCGGTGGCCTGAAAGACATCGCAATGCCGGAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATa                                                                                  <  1:2017060/140‑1 (MQ=255)
                                                          gaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATa                                                                                  <  2:2072859/140‑1 (MQ=255)
                                                           aaCGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATAc                                                                                 >  2:2332450/1‑140 (MQ=255)
                                                           aaCGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATAc                                                                                 >  1:3035267/1‑140 (MQ=255)
                                                               ccGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCtt                                                                             >  2:1062566/1‑140 (MQ=255)
                                                                    ggCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTcc                                                                        <  2:3035267/140‑1 (MQ=255)
                                                                     gCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTccc                                                                       <  1:2332450/140‑1 (MQ=255)
                                                                        tGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTgg                                                                    >  1:860292/1‑140 (MQ=255)
                                                                          aaaGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTa                                                                  >  1:1274043/1‑140 (MQ=255)
                                                                          aaaGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTa                                                                  >  2:1544927/1‑140 (MQ=255)
                                                                            aGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaa                                                                <  1:1428675/140‑1 (MQ=255)
                                                                               cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCAAACCtt                                                                                                                          <  2:484986/79‑1 (MQ=255)
                                                                               cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCAAACCtt                                                                                                                          >  1:484986/1‑79 (MQ=255)
                                                                                    cAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAtt                                                        <  2:3161381/140‑1 (MQ=255)
                                                                                          cGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGc                                                  <  1:2018399/140‑1 (MQ=255)
                                                                                           gAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCg                                                 <  1:1062566/140‑1 (MQ=255)
                                                                                           gAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCg                                                 <  1:813980/140‑1 (MQ=255)
                                                                                            aaGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCgg                                                >  1:703381/1‑139 (MQ=255)
                                                                                               gTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGAc                                             <  2:1544985/140‑1 (MQ=255)
                                                                                               gTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGAc                                             >  2:528035/1‑140 (MQ=255)
                                                                                                  aaaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTaa                                          <  1:528035/140‑1 (MQ=255)
                                                                                                   aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTaaa                                         >  1:496971/1‑140 (MQ=255)
                                                                                                   aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTaaa                                         >  1:410867/1‑140 (MQ=255)
                                                                                                   aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTaaa                                         >  2:3460147/1‑140 (MQ=255)
                                                                                                         aaaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCtt                                   >  2:2478569/1‑140 (MQ=255)
                                                                                                         aaaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCtt                                   <  2:496971/140‑1 (MQ=255)
                                                                                                         aaaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTAACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCtt                                   >  1:2203002/1‑140 (MQ=255)
                                                                                                          aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTtg                                  >  1:2265592/1‑140 (MQ=255)
                                                                                                           aaGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTtgt                                 <  2:1940800/140‑1 (MQ=255)
                                                                                                           aaGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTtgt                                 <  1:1904160/140‑1 (MQ=255)
                                                                                                                ttCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTaa                                          <  1:49009/126‑1 (MQ=255)
                                                                                                                ttCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCa                            >  1:2698848/1‑140 (MQ=255)
                                                                                                                ttCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCa                            >  2:31920/1‑140 (MQ=255)
                                                                                                                ttCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCa                            >  1:2856181/1‑140 (MQ=255)
                                                                                                                  cACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAt                          <  1:2478569/140‑1 (MQ=255)
                                                                                                                    ccATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTg                        <  2:1276816/140‑1 (MQ=255)
                                                                                                                    ccATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTg                        <  2:2856181/140‑1 (MQ=255)
                                                                                                                    ccATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTg                        <  1:1544927/140‑1 (MQ=255)
                                                                                                                      aTGAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGAc                      <  1:993551/140‑1 (MQ=255)
                                                                                                                        gAGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAg                    >  2:556998/1‑140 (MQ=255)
                                                                                                                         aGCGGCGCTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGt                   >  1:915356/1‑140 (MQ=255)
                                                                                                                             gcgcTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAt               >  2:333353/1‑140 (MQ=255)
                                                                                                                              cgcTTATGCGCCGATTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAt               >  2:64839/1‑140 (MQ=255)
                                                                                                                              cgcTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAtt              <  2:1274043/140‑1 (MQ=255)
                                                                                                                               gcTTATGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAttt             >  2:153968/1‑140 (MQ=255)
                                                                                                                                   aTGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTc           <  2:1028267/138‑1 (MQ=255)
                                                                                                                                   aTGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTc           >  1:1028267/1‑138 (MQ=255)
                                                                                                                                   aTGCGCCG‑‑TTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTg         >  2:304617/1‑140 (MQ=255)
                                                                                                                                            ttttttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGagta  <  1:64839/140‑1 (MQ=255)
                                                                                                                                           ||                                                                                                                                           
GACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTA  >  NC_000913/1369897‑1370176

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: