Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,370,026 T→A 6.2% intergenic (+23/‑190) 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,0260TA 6.2% 408.5 / 10.6 122intergenic (+23/‑190)pspE/ycjMthiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein
Reads supporting (aligned to +/- strand):  ref base T (54/60);  new base A (4/4);  total (58/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-01

Only 100 of 127 total aligned reads displayed.

AAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGA  >  NC_000913/1369879‑1370173
                                                                                                                                                   |                                                                                                                                                      
aaGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCtt                                                                                                                                                       >  4:498102/1‑149 (MQ=255)
aaGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCtt                                                                                                                                                       <  4:790918/149‑1 (MQ=255)
aaGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTGAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCtt                                                                                                                                                       <  4:141828/149‑1 (MQ=255)
  gtgtATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTAt                                                                                                                                                     <  3:640642/149‑1 (MQ=255)
   tgtATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATg                                                                                                                                                    >  4:210709/1‑149 (MQ=255)
    gtATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATgc                                                                                                                                                   <  1:1010188/149‑1 (MQ=255)
         gCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑t                                                                                                                                             <  4:278059/149‑1 (MQ=255)
           aaTGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttt                                                                                                                                           <  2:53902/149‑1 (MQ=255)
           aaTGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttt                                                                                                                                           >  3:911653/1‑149 (MQ=255)
           aaTGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttt                                                                                                                                           >  3:954299/1‑149 (MQ=255)
            aTGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑tttt                                                                                                                                          <  3:568971/149‑1 (MQ=255)
            aTGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctca                                                                                                                                           <  3:769686/149‑5 (MQ=255)
             tGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttttt                                                                                                                                         <  4:200858/149‑1 (MQ=255)
              gCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑tttttt                                                                                                                                        <  1:287483/149‑1 (MQ=255)
              gCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑tttttt                                                                                                                                        >  3:314180/1‑149 (MQ=255)
              gCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑tttttt                                                                                                                                        <  1:36157/149‑1 (MQ=255)
               ccGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatgg                                                                                                                                        <  4:408364/149‑8 (MQ=255)
                cGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑tttttttt                                                                                                                                      <  3:876843/149‑1 (MQ=255)
                  gACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAGCGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTc                                                                                                                                    <  4:858874/149‑1 (MQ=255)
                  gACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTc                                                                                                                                    <  1:904173/149‑1 (MQ=255)
                   aCGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCt                                                                                                                                   <  4:702778/149‑1 (MQ=255)
                     gCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑Gt                                                                                                                                <  2:702243/149‑1 (MQ=255)
                       cAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTca                                                                                                                              >  3:389636/1‑149 (MQ=255)
                       cAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGCTGAGAACGCCGGGGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGGTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTGTTCT‑GTca                                                                                                                              >  2:498689/1‑149 (MQ=255)
                                caagcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTAtt                                                                                                                     <  3:611604/149‑1 (MQ=255)
                                 aagcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTAttt                                                                                                                    >  1:885453/1‑149 (MQ=255)
                                 aagcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTAttt                                                                                                                    >  1:53707/1‑149 (MQ=255)
                                 aagcaaAAGAGAGCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCAGGAGCGGCGCTTATGCGCCGTTTTTTTTTTCT‑GTCACACCTTAtt                                                                                                                     <  1:996742/149‑1 (MQ=255)
                                         gagaTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATcccc                                                                                                            <  4:293905/149‑1 (MQ=255)
                                           gaTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCAt                                                                                                          <  4:389636/149‑1 (MQ=255)
                                            aTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACTATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATa                                                                                                         <  4:257062/149‑1 (MQ=255)
                                             tCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAg                                                                                                        >  1:603182/1‑149 (MQ=255)
                                                ttAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTAt                                                                                                                                                     >  1:1054562/1‑103 (MQ=255)
                                                ttAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTAt                                                                                                                                                     <  2:1054562/103‑1 (MQ=255)
                                                ttAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATTGatt                                                                                                     >  2:300436/1‑149 (MQ=255)
                                                    cGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGattatt                                                                                                  <  1:1053798/149‑1 (MQ=255)
                                                    cGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTAttt                                                                                                 <  3:180934/149‑1 (MQ=255)
                                                    cGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTAttt                                                                                                 <  2:781195/149‑1 (MQ=255)
                                                        aTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCtt                                                                                                                                                       <  4:181884/93‑1 (MQ=255)
                                                          gggATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCa                                                                                           >  1:779069/1‑149 (MQ=255)
                                                           ggATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAg                                                                                          <  2:932411/149‑1 (MQ=255)
                                                            gATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGc                                                                                         <  3:780170/149‑1 (MQ=255)
                                                                     aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATa                                                                               <  4:1034259/149‑1 (MQ=255)
                                                                     aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATa                                                                               >  3:695589/1‑149 (MQ=255)
                                                                     aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATa                                                                               >  2:805575/1‑149 (MQ=255)
                                                                            gaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCtttt                                                                        >  2:956211/1‑149 (MQ=255)
                                                                              aCGCCGTTGGCTTGAACGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGCTCGTCCCTTATGCGCCG‑TTTTTTTTTCT‑ATCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑ACATACGCCTTTTc                                                                      <  2:1051709/149‑1 (MQ=255)
                                                                                gCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTAtt                                                                                                                     <  1:22277/101‑1 (MQ=255)
                                                                                gCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTAtt                                                                                                                     >  2:22277/1‑101 (MQ=255)
                                                                                 ccGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCt                                                                   <  4:722276/149‑1 (MQ=255)
                                                                                  cGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTg                                                                  <  1:459322/149‑1 (MQ=255)
                                                                                      ggCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTaa                                                              >  1:248003/1‑149 (MQ=255)
                                                                                         cTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTaaaaa                                                           >  3:366574/1‑149 (MQ=255)
                                                                                          tGAAAGACATCGCGATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTaaaaaa                                                          <  4:273619/149‑1 (MQ=255)
                                                                                          tGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTaaaaaa                                                          <  1:331720/149‑1 (MQ=255)
                                                                                               gACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGAtt                                                     >  4:481563/1‑149 (MQ=255)
                                                                                                 cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTc                                                   <  1:830161/149‑1 (MQ=255)
                                                                                                  aTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGTGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTcc                                                  >  4:691939/1‑149 (MQ=255)
                                                                                                  aTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTc                                                   <  1:486510/148‑1 (MQ=255)
                                                                                                  aTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTc                                                   >  2:486510/1‑148 (MQ=255)
                                                                                                   tCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCt                                                 >  3:397129/1‑149 (MQ=255)
                                                                                                   tCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCt                                                 <  2:885453/149‑1 (MQ=255)
                                                                                                      cAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCg                                              <  4:386022/149‑1 (MQ=255)
                                                                                                         tGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGa                                           <  3:749918/149‑1 (MQ=255)
                                                                                                          gCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTTTGGCCACACCTTATTAACACCCCCTTAGATTATTGGTGTC‑GCCCACAAAAAACCTTTTTTCCCTGTAAAAAATTGTTTTCCTGcttac                                           >  4:492634/1‑145 (MQ=255)
                                                                                                            cGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTa                                        >  2:902629/1‑149 (MQ=255)
                                                                                                             gAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTaa                                       >  4:532058/1‑149 (MQ=255)
                                                                                                                   cAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCt                                 >  2:751744/1‑149 (MQ=255)
                                                                                                                    aaaGGTTAAAAGAATTGACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTGTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTg                                                                  <  3:654302/115‑1 (MQ=255)
                                                                                                                    aaaGGTTAAAAGAATTCCCCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCtt                                <  3:822131/149‑1 (MQ=255)
                                                                                                                    aaaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCtt                                >  1:582067/1‑149 (MQ=255)
                                                                                                                    aaaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCtt                                <  4:933844/149‑1 (MQ=255)
                                                                                                                    aaaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCtt                                <  1:751744/149‑1 (MQ=255)
                                                                                                                    aaaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCACCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTg                                                                  >  4:654302/1‑115 (MQ=255)
                                                                                                                        gTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTaaaa                                     >  4:1032426/1‑140 (MQ=255)
                                                                                                                        gTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTaaaa                                     <  3:1032426/140‑1 (MQ=255)
                                                                                                                         ttAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCt                           >  3:497063/1‑149 (MQ=255)
                                                                                                                           aaaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTtgt                              <  1:1006848/144‑1 (MQ=255)
                                                                                                                           aaaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCa                         >  3:112759/1‑149 (MQ=255)
                                                                                                                           aaaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTAATTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCtttt                              >  2:1006848/1‑142 (MQ=255)
                                                                                                                            aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCaa                        >  1:634159/1‑149 (MQ=255)
                                                                                                                             aaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAt                       <  2:1063953/149‑1 (MQ=255)
                                                                                                                             aaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAt                       <  1:158633/149‑1 (MQ=255)
                                                                                                                             aaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAt                       >  1:431598/1‑149 (MQ=255)
                                                                                                                               gAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGT‑AAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTg                     >  1:72843/1‑149 (MQ=255)
                                                                                                                               gAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCacac                                                                                                                           <  1:173714/48‑1 (MQ=255)
                                                                                                                               gAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTg                     >  4:23273/1‑149 (MQ=255)
                                                                                                                                        aTGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAt            >  1:309677/1‑149 (MQ=255)
                                                                                                                                        aTGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAt            >  1:715395/1‑149 (MQ=255)
                                                                                                                                         tGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAtt           >  1:53231/1‑149 (MQ=255)
                                                                                                                                         tGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAtt           <  1:956211/149‑1 (MQ=255)
                                                                                                                                           aacggcgCATAAGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAtttt         >  2:397742/3‑149 (MQ=255)
                                                                                                                                           aacggcgCATAAGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAtttt         >  2:792676/3‑149 (MQ=255)
                                                                                                                                           aacggcgCATAAGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAtttt         >  3:926308/3‑149 (MQ=255)
                                                                                                                                            gcggcgCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTc        <  2:103156/149‑1 (MQ=255)
                                                                                                                                               gcgcTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGc     >  2:347114/1‑149 (MQ=255)
                                                                                                                                                 gcTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCgg   >  3:1004682/1‑149 (MQ=255)
                                                                                                                                                  cTTATGCGCCGTTTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGT‑AAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGa  <  2:72843/149‑1 (MQ=255)
                                                                                                                                                  cTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGa  <  3:481563/149‑1 (MQ=255)
                                                                                                                                                  cTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGa  <  2:715395/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
AAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGA  >  NC_000913/1369879‑1370173

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