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.

AAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGA  >  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)
         gCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGt                                                                                                                                             <  4:278059/149‑1 (MQ=255)
           aaTGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttt                                                                                                                                           <  2:53902/149‑1 (MQ=255)
           aaTGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttt                                                                                                                                           >  3:911653/1‑149 (MQ=255)
           aaTGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttt                                                                                                                                           >  3:954299/1‑149 (MQ=255)
            aTGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGtttt                                                                                                                                          <  3:568971/149‑1 (MQ=255)
            aTGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctca                                                                                                                                           <  3:769686/149‑5 (MQ=255)
             tGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttt                                                                                                                                         <  4:200858/149‑1 (MQ=255)
              gCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGtttttt                                                                                                                                        <  1:287483/149‑1 (MQ=255)
              gCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGtttttt                                                                                                                                        >  3:314180/1‑149 (MQ=255)
              gCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGtttttt                                                                                                                                        <  1:36157/149‑1 (MQ=255)
               ccGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatgg                                                                                                                                        <  4:408364/149‑8 (MQ=255)
                cGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGtttttttt                                                                                                                                      <  3:876843/149‑1 (MQ=255)
                  gACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAGCGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTc                                                                                                                                    <  4:858874/149‑1 (MQ=255)
                  gACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTc                                                                                                                                    <  1:904173/149‑1 (MQ=255)
                   aCGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCt                                                                                                                                   <  4:702778/149‑1 (MQ=255)
                     gCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGt                                                                                                                                <  2:702243/149‑1 (MQ=255)
                       cAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTca                                                                                                                              >  3:389636/1‑149 (MQ=255)
                       cAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGCTGAGAACGCCGGGGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGGTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTGTTCTGTca                                                                                                                              >  2:498689/1‑149 (MQ=255)
                                caagcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTAtt                                                                                                                     <  3:611604/149‑1 (MQ=255)
                                 aagcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTAttt                                                                                                                    >  1:885453/1‑149 (MQ=255)
                                 aagcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTAttt                                                                                                                    >  1:53707/1‑149 (MQ=255)
                                 aagcaaAAGAGAGCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCAGGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTAtt                                                                                                                     <  1:996742/149‑1 (MQ=255)
                                         gagaTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATcccc                                                                                                            <  4:293905/149‑1 (MQ=255)
                                           gaTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCAt                                                                                                          <  4:389636/149‑1 (MQ=255)
                                            aTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACTATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATa                                                                                                         <  4:257062/149‑1 (MQ=255)
                                             tCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAg                                                                                                        >  1:603182/1‑149 (MQ=255)
                                                ttAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTAt                                                                                                                                                     >  1:1054562/1‑103 (MQ=255)
                                                ttAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTAt                                                                                                                                                     <  2:1054562/103‑1 (MQ=255)
                                                ttAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATTGatt                                                                                                     >  2:300436/1‑149 (MQ=255)
                                                    cGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGattatt                                                                                                  <  1:1053798/149‑1 (MQ=255)
                                                    cGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTAttt                                                                                                 <  3:180934/149‑1 (MQ=255)
                                                    cGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTAttt                                                                                                 <  2:781195/149‑1 (MQ=255)
                                                        aTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCtt                                                                                                                                                       <  4:181884/93‑1 (MQ=255)
                                                          gggATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCa                                                                                           >  1:779069/1‑149 (MQ=255)
                                                           ggATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAg                                                                                          <  2:932411/149‑1 (MQ=255)
                                                            gATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGc                                                                                         <  3:780170/149‑1 (MQ=255)
                                                                     aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATa                                                                               <  4:1034259/149‑1 (MQ=255)
                                                                     aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATa                                                                               >  3:695589/1‑149 (MQ=255)
                                                                     aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATa                                                                               >  2:805575/1‑149 (MQ=255)
                                                                            gaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCtttt                                                                        >  2:956211/1‑149 (MQ=255)
                                                                              aCGCCGTTGGCTTGAACGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGCTCGTCCCTTATGCGCCGTTTTTTTTTCTATCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACACATACGCCTTTTc                                                                      <  2:1051709/149‑1 (MQ=255)
                                                                                gCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTAtt                                                                                                                     <  1:22277/101‑1 (MQ=255)
                                                                                gCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTAtt                                                                                                                     >  2:22277/1‑101 (MQ=255)
                                                                                 ccGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCt                                                                   <  4:722276/149‑1 (MQ=255)
                                                                                  cGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTg                                                                  <  1:459322/149‑1 (MQ=255)
                                                                                      ggCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaa                                                              >  1:248003/1‑149 (MQ=255)
                                                                                         cTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaaaa                                                           >  3:366574/1‑149 (MQ=255)
                                                                                          tGAAAGACATCGCGATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaaaaa                                                          <  4:273619/149‑1 (MQ=255)
                                                                                          tGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaaaaa                                                          <  1:331720/149‑1 (MQ=255)
                                                                                               gACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAtt                                                     >  4:481563/1‑149 (MQ=255)
                                                                                                 cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTc                                                   <  1:830161/149‑1 (MQ=255)
                                                                                                  aTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGTGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTcc                                                  >  4:691939/1‑149 (MQ=255)
                                                                                                  aTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTc                                                   <  1:486510/148‑1 (MQ=255)
                                                                                                  aTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTc                                                   >  2:486510/1‑148 (MQ=255)
                                                                                                   tCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCt                                                 >  3:397129/1‑149 (MQ=255)
                                                                                                   tCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCt                                                 <  2:885453/149‑1 (MQ=255)
                                                                                                      cAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCg                                              <  4:386022/149‑1 (MQ=255)
                                                                                                         tGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGa                                           <  3:749918/149‑1 (MQ=255)
                                                                                                          gCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTTGGCCACACCTTATTAACACCCCCTTAGATTATTGGTGTCGCCCACAAAAAACCTTTTTTCCCTGTAAAAAATTGTTTTCCTGcttac                                           >  4:492634/1‑145 (MQ=255)
                                                                                                            cGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTa                                        >  2:902629/1‑149 (MQ=255)
                                                                                                             gAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTaa                                       >  4:532058/1‑149 (MQ=255)
                                                                                                                   cAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCt                                 >  2:751744/1‑149 (MQ=255)
                                                                                                                    aaaGGTTAAAAGAATTGACCATGAGCGGCGCTTATGCGCCGTTTTTTTGTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTg                                                                  <  3:654302/115‑1 (MQ=255)
                                                                                                                    aaaGGTTAAAAGAATTCCCCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCtt                                <  3:822131/149‑1 (MQ=255)
                                                                                                                    aaaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCtt                                >  1:582067/1‑149 (MQ=255)
                                                                                                                    aaaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCtt                                <  4:933844/149‑1 (MQ=255)
                                                                                                                    aaaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCtt                                <  1:751744/149‑1 (MQ=255)
                                                                                                                    aaaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCACCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTg                                                                  >  4:654302/1‑115 (MQ=255)
                                                                                                                        gTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTaaaa                                     >  4:1032426/1‑140 (MQ=255)
                                                                                                                        gTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTaaaa                                     <  3:1032426/140‑1 (MQ=255)
                                                                                                                         ttAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCt                           >  3:497063/1‑149 (MQ=255)
                                                                                                                           aaaaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTtgt                              <  1:1006848/144‑1 (MQ=255)
                                                                                                                           aaaaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCa                         >  3:112759/1‑149 (MQ=255)
                                                                                                                           aaaaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTAATTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCtttt                              >  2:1006848/1‑142 (MQ=255)
                                                                                                                            aaaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCaa                        >  1:634159/1‑149 (MQ=255)
                                                                                                                             aaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAt                       <  2:1063953/149‑1 (MQ=255)
                                                                                                                             aaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAt                       <  1:158633/149‑1 (MQ=255)
                                                                                                                             aaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAt                       >  1:431598/1‑149 (MQ=255)
                                                                                                                               gAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTg                     >  1:72843/1‑149 (MQ=255)
                                                                                                                               gAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCacac                                                                                                                           <  1:173714/48‑1 (MQ=255)
                                                                                                                               gAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTg                     >  4:23273/1‑149 (MQ=255)
                                                                                                                                        aTGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAt            >  1:309677/1‑149 (MQ=255)
                                                                                                                                        aTGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAt            >  1:715395/1‑149 (MQ=255)
                                                                                                                                         tGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAtt           >  1:53231/1‑149 (MQ=255)
                                                                                                                                         tGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAtt           <  1:956211/149‑1 (MQ=255)
                                                                                                                                           aacggcgCATAAGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAtttt         >  2:397742/3‑149 (MQ=255)
                                                                                                                                           aacggcgCATAAGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAtttt         >  2:792676/3‑149 (MQ=255)
                                                                                                                                           aacggcgCATAAGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAtttt         >  3:926308/3‑149 (MQ=255)
                                                                                                                                            gcggcgCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTc        <  2:103156/149‑1 (MQ=255)
                                                                                                                                               gcgcTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGc     >  2:347114/1‑149 (MQ=255)
                                                                                                                                                 gcTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCgg   >  3:1004682/1‑149 (MQ=255)
                                                                                                                                                  cTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGa  <  2:72843/149‑1 (MQ=255)
                                                                                                                                                  cTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGa  <  3:481563/149‑1 (MQ=255)
                                                                                                                                                  cTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGa  <  2:715395/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
AAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGA  >  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: