Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,370,029 T→A 6.5% intergenic (+26/‑187) 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,0290TA 6.5% 411.7 / 13.9 122intergenic (+26/‑187)pspE/ycjMthiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein
Reads supporting (aligned to +/- strand):  ref base T (53/61);  new base A (4/4);  total (57/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

Only 100 of 125 total aligned reads displayed.

GTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTA  >  NC_000913/1369881‑1370176
                                                                                                                                                    |                                                                                                                                                      
gtgtATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTAt                                                                                                                                                        <  3:640642/149‑1 (MQ=255)
 tgtATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATg                                                                                                                                                       >  4:210709/1‑149 (MQ=255)
  gtATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATgc                                                                                                                                                      <  1:1010188/149‑1 (MQ=255)
       gCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑t                                                                                                                                                >  2:88763/1‑149 (MQ=255)
       gCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑t                                                                                                                                                <  4:278059/149‑1 (MQ=255)
         aaTGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttt                                                                                                                                              >  3:911653/1‑149 (MQ=255)
         aaTGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttt                                                                                                                                              <  2:53902/149‑1 (MQ=255)
          aTGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctca                                                                                                                                              <  3:769686/149‑5 (MQ=255)
            gCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑tttttt                                                                                                                                           <  1:287483/149‑1 (MQ=255)
            gCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑tttttt                                                                                                                                           >  3:314180/1‑149 (MQ=255)
             ccGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatgg                                                                                                                                           <  4:237649/149‑8 (MQ=255)
             ccGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatgg                                                                                                                                           <  4:408364/149‑8 (MQ=255)
             ccGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatgg                                                                                                                                           <  2:561367/149‑8 (MQ=255)
              cGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑tttttttt                                                                                                                                         <  3:876843/149‑1 (MQ=255)
               ggACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttttttttt                                                                                                                                        <  2:987922/149‑1 (MQ=255)
                gACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAGCGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTc                                                                                                                                       <  4:858874/149‑1 (MQ=255)
                 aCGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCt                                                                                                                                      <  4:702778/149‑1 (MQ=255)
                   gCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑Gt                                                                                                                                   <  2:702243/149‑1 (MQ=255)
                     cAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGCTGAGAACGCCGGGGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGGTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTGTTCT‑GTca                                                                                                                                 >  2:498689/1‑149 (MQ=255)
                              caagcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTAtt                                                                                                                        <  4:314180/149‑1 (MQ=255)
                              caagcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTAtt                                                                                                                        <  3:611604/149‑1 (MQ=255)
                               aagcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTAttt                                                                                                                       >  4:180934/1‑149 (MQ=255)
                               aagcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTAttt                                                                                                                       >  1:53707/1‑149 (MQ=255)
                               aagcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTAttt                                                                                                                       >  1:885453/1‑149 (MQ=255)
                               aagcaaAAGAGAGCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCAGGAGCGGCGCTTATGCGCCGTTTTTTTTTTCT‑GTCACACCTTAtt                                                                                                                        <  1:996742/149‑1 (MQ=255)
                                    aaaGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATc                                                                                                                  >  3:722276/1‑149 (MQ=255)
                                       gagaTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATcccc                                                                                                               <  4:293905/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                                                                                                    <  2:781195/149‑1 (MQ=255)
                                                  cGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTAttt                                                                                                    <  3:180934/149‑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)
                                                                  cACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAAt                                                                                   >  4:1032303/1‑149 (MQ=255)
                                                                   gcGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGTCGAAGGTCAAAGGTTAAAAGAATTCACCATGGGCGGCTCTAATGCACCG‑TTTTTTTTTCT‑GTCACCCATAATTTACATCCCAATAGAGTATTGGTTTCAGCCCAC‑AAATa                                                                                  <  2:1028388/148‑1 (MQ=255)
                                                                   aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATa                                                                                  <  4:1034259/149‑1 (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‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTccc                                                                       <  2:779069/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)
                                                                                        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)
                                                                                                 tCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCt                                                    <  2:885453/149‑1 (MQ=255)
                                                                                                 tCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCt                                                    >  3:397129/1‑149 (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)
                                                                                                              ggTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAAc                                       <  1:88763/149‑1 (MQ=255)
                                                                                                                 cAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCt                                    >  2:751744/1‑149 (MQ=255)
                                                                                                                  aaaGGTTAAAAGAATTGACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTGTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTg                                                                     <  3:654302/115‑1 (MQ=255)
                                                                                                                  aaaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCtt                                   <  4:933844/149‑1 (MQ=255)
                                                                                                                  aaaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCtt                                   >  1:582067/1‑149 (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)
                                                                                                                       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)
                                                                                                                         aaaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGAGTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCa                            >  4:559689/1‑149 (MQ=255)
                                                                                                                          aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCaa                           <  3:583823/149‑1 (MQ=255)
                                                                                                                          aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCaa                           >  1:634159/1‑149 (MQ=255)
                                                                                                                           aaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAt                          >  1:431598/1‑149 (MQ=255)
                                                                                                                           aaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAt                          <  2:1063953/149‑1 (MQ=255)
                                                                                                                           aaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAt                          <  1:158633/149‑1 (MQ=255)
                                                                                                                             gAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGT‑AAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTg                        >  1:72843/1‑149 (MQ=255)
                                                                                                                             gAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCacac                                                                                                                              >  2:173714/1‑48 (MQ=255)
                                                                                                                             gAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCacac                                                                                                                              <  1:173714/48‑1 (MQ=255)
                                                                                                                             gAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTg                        >  4:23273/1‑149 (MQ=255)
                                                                                                                               aTTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGAc                      >  3:859052/1‑149 (MQ=255)
                                                                                                                                     cATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTa                >  2:425835/1‑149 (MQ=255)
                                                                                                                                     cATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTa                >  4:42570/1‑149 (MQ=255)
                                                                                                                                     cATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTa                >  1:365118/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            >  3:926308/3‑149 (MQ=255)
                                                                                                                                          acggcgCATAAGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTc           >  4:768372/2‑149 (MQ=255)
                                                                                                                                             gcgcTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGc        >  2:347114/1‑149 (MQ=255)
                                                                                                                                               gcTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCgg      >  3:1004682/1‑149 (MQ=255)
                                                                                                                                               gcTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCgg      >  3:403895/1‑149 (MQ=255)
                                                                                                                                                cTTATGCGCCGTTTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGT‑AAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGa     <  2:72843/149‑1 (MQ=255)
                                                                                                                                                cTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGa     <  3:481563/149‑1 (MQ=255)
                                                                                                                                                  tATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGagt   <  1:300436/149‑1 (MQ=255)
                                                                                                                                                   aTGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGagta  <  2:846912/149‑1 (MQ=255)
                                                                                                                                                   aTGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGagta  <  4:954299/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAC‑AAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTA  >  NC_000913/1369881‑1370176

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