Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,370,026 T→A 6.1% 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.1% 503.6 / 17.7 147intergenic (+23/‑190)pspE/ycjMthiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein
Reads supporting (aligned to +/- strand):  ref base T (74/63);  new base A (6/3);  total (81/66)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

Only 100 of 151 total aligned reads displayed.

AAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAG  >  NC_000913/1369879‑1370174
                                                                                                                                                   |                                                                                                                                                      
aaGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCtt                                                                                                                                                       >  1:308231/1‑149 (MQ=255)
aaGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCtt                                                                                                                                                       <  4:1392297/149‑1 (MQ=255)
 aGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTa                                                                                                                                                      <  3:98700/149‑1 (MQ=255)
  gtgtATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTAt                                                                                                                                                     <  4:729355/149‑1 (MQ=255)
  gtgtATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTAt                                                                                                                                                     <  4:565479/149‑1 (MQ=255)
   tgtATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATg                                                                                                                                                    >  4:720828/1‑149 (MQ=255)
    gtATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATgc                                                                                                                                                   <  3:864078/149‑1 (MQ=255)
      aTTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATgcgc                                                                                                                                                 >  2:825897/1‑149 (MQ=255)
       ttGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGcc                                                                                                                                                >  4:472372/1‑149 (MQ=255)
         gCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑t                                                                                                                                            >  1:1315245/1‑149 (MQ=255)
          cAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑tt                                                                                                                                           <  3:1186112/149‑1 (MQ=255)
           aaTGCCGGACGCCAGATAGGGCAAGCAGAAGATATCCTTAGCGAGATGGGCTATACCCACGATGAGAACGCCGGTGGCCTGAAAAACATCGCAATGCCGAAGGTCAAAGGCTAAAAGAATTCACCATGAGCGGCGCGTACGAGCCG‑‑ttt                                                                                                                                          >  2:1314579/1‑149 (MQ=255)
             tGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑ttttt                                                                                                                                        >  1:359181/1‑149 (MQ=255)
               ccGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑ttttttt                                                                                                                                      <  3:90151/149‑1 (MQ=255)
                 ggACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑ttttttttt                                                                                                                                    >  3:473908/1‑149 (MQ=255)
                          tCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCacac                                                                                                                           <  2:575267/149‑1 (MQ=255)
                           cAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACAcc                                                                                                                          <  4:288880/149‑1 (MQ=255)
                              ggCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTa                                                                                                                       >  1:372551/1‑149 (MQ=255)
                               gcaagcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTAt                                                                                                                      >  3:48821/1‑149 (MQ=255)
                               gcaagcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTAt                                                                                                                      >  4:805953/1‑149 (MQ=255)
                                 aagcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTAttt                                                                                                                    >  3:169001/1‑149 (MQ=255)
                                   gcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTAc                                                                                                                  >  1:1581395/1‑149 (MQ=255)
                                     aaaaGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATgcgc                                                                                                                                                 >  4:338717/1‑118 (MQ=255)
                                     aaaaGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATgcgc                                                                                                                                                 <  3:338717/118‑1 (MQ=255)
                                     aaaaGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTAttt                                                                                                                    >  3:941401/1‑145 (MQ=255)
                                     aaaaGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACAt                                                                                                                >  4:1035089/1‑149 (MQ=255)
                                            aTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATa                                                                                                         >  2:470013/1‑149 (MQ=255)
                                            aTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATa                                                                                                         >  3:792425/1‑149 (MQ=255)
                                            aTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAAACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATa                                                                                                         <  2:967532/149‑1 (MQ=255)
                                             tCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAg                                                                                                        <  1:669386/149‑1 (MQ=255)
                                             tCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAg                                                                                                        >  1:1397608/1‑149 (MQ=255)
                                              ccTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGa                                                                                                       >  2:1139140/1‑149 (MQ=255)
                                                ttAGCGAGATGGGATATACCCACGTTTAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTACGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGatt                                                                                                     >  1:956870/1‑149 (MQ=255)
                                                ttAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAAGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGatt                                                                                                     >  4:1218513/1‑149 (MQ=255)
                                                    cGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTTCTTTTACACCTTATTTACATCCCCATAAattat                                                                                                   >  4:945808/1‑149 (MQ=255)
                                                      agaTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGc                                                                                               >  2:1370978/1‑149 (MQ=255)
                                                            gATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGc                                                                                         >  1:832645/1‑149 (MQ=255)
                                                                   ccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTAAACCCTATTTTCCTTCCCCATAGTTTTTTTGCTTCAGCTccaaa                                                                                   >  4:900303/1‑145 (MQ=255)
                                                                   ccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTca                                                                                                                              >  3:1338660/1‑105 (MQ=255)
                                                                   ccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaaa                                                                                  >  1:110313/1‑149 (MQ=255)
                                                                   ccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaaa                                                                                  >  4:79433/1‑149 (MQ=255)
                                                                   ccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaaa                                                                                  >  2:1041577/1‑149 (MQ=255)
                                                                    cACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAAt                                                                                 >  4:380051/1‑149 (MQ=255)
                                                                      cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATAc                                                                               <  4:758437/149‑1 (MQ=255)
                                                                       gTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGc                                                                                               <  3:960032/132‑1 (MQ=255)
                                                                       gTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGc                                                                                               >  4:960032/1‑132 (MQ=255)
                                                                        ttGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGc                                                                             <  1:1121712/149‑1 (MQ=255)
                                                                                 ccGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCt                                                                    <  2:1315245/149‑1 (MQ=255)
                                                                                 ccGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCt                                                                    <  1:746479/149‑1 (MQ=255)
                                                                                   ggTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTgg                                                                  <  2:846246/149‑1 (MQ=255)
                                                                                   ggTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTgg                                                                  <  3:217541/149‑1 (MQ=255)
                                                                                         cggAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGGGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCCCAAATACGCTTTTTCCCTGCTaaaaa                                                            >  2:474197/3‑149 (MQ=255)
                                                                                          tGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaaaaa                                                           >  3:978596/1‑149 (MQ=255)
                                                                                          tGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaaaaa                                                           >  2:185475/1‑149 (MQ=255)
                                                                                          tGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaaaaa                                                           >  4:32278/1‑149 (MQ=255)
                                                                                           gAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑ttttttttt                                                                                                                                    >  1:932139/1‑75 (MQ=255)
                                                                                              aGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAt                                                       <  2:956870/149‑1 (MQ=255)
                                                                                                  aTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTccc                                                                     <  4:921835/131‑1 (MQ=255)
                                                                                                  aTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTcc                                                   >  1:1003063/1‑149 (MQ=255)
                                                                                                   tCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGa                                                        >  3:1202505/1‑143 (MQ=255)
                                                                                                   tCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCt                                                  <  4:1348366/149‑1 (MQ=255)
                                                                                                    cGCAATGCCGAAGGTCAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGc                                                <  4:765761/149‑1 (MQ=255)
                                                                                                     gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGc                                                <  1:594950/149‑1 (MQ=255)
                                                                                                      cAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCg                                               >  3:287584/1‑149 (MQ=255)
                                                                                                            cGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACt                                          >  2:1464132/1‑149 (MQ=255)
                                                                                                            cGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTa                                         >  4:330249/1‑149 (MQ=255)
                                                                                                               aGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTaaaa                                      <  4:169001/149‑1 (MQ=255)
                                                                                                               aGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTaaaa                                      <  2:830562/149‑1 (MQ=255)
                                                                                                                 gTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAAcc                                    >  4:196741/1‑149 (MQ=255)
                                                                                                                  tCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAAccc                                   >  4:940188/1‑149 (MQ=255)
                                                                                                                   cAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTCTTTCCATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCt                                  >  1:1029069/1‑149 (MQ=255)
                                                                                                                   cAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCt                                  <  4:287584/149‑1 (MQ=255)
                                                                                                                   cAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCt                                  >  2:1368888/1‑149 (MQ=255)
                                                                                                                    aaaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGc                                                >  3:637940/1‑134 (MQ=255)
                                                                                                                    aaaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGc                                                <  4:637940/134‑1 (MQ=255)
                                                                                                                    aaaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCtt                                 <  3:919427/149‑1 (MQ=255)
                                                                                                                     aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTtg                                <  1:1370978/149‑1 (MQ=255)
                                                                                                                      aGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTGCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTtgt                               <  4:978596/149‑1 (MQ=255)
                                                                                                                       ggTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTtgtg                              <  3:1218513/149‑1 (MQ=255)
                                                                                                                            aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCaa                         >  3:646877/1‑149 (MQ=255)
                                                                                                                            aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCaa                         <  2:561632/149‑1 (MQ=255)
                                                                                                                             aaGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAt                        <  3:1333529/149‑1 (MQ=255)
                                                                                                                             aaGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCACCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCTCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAt                        >  4:16086/1‑149 (MQ=255)
                                                                                                                              aGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAtt                       <  3:834203/149‑1 (MQ=255)
                                                                                                                               gAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAtt                       <  1:377889/149‑1 (MQ=255)
                                                                                                                                aaTTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGa                     <  2:856214/149‑1 (MQ=255)
                                                                                                                                  ttCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGGCa                   >  4:113596/1‑149 (MQ=255)
                                                                                                                                  ttCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACa                   >  4:203713/1‑149 (MQ=255)
                                                                                                                                     aCCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGtt                >  4:411048/1‑149 (MQ=255)
                                                                                                                                       cATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTa              >  1:727389/1‑149 (MQ=255)
                                                                                                                                       aaaaaacggcgCATAAGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACa                   <  3:773800/138‑1 (MQ=255)
                                                                                                                                          gAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAttt           >  4:902941/1‑149 (MQ=255)
                                                                                                                                          aaacggcgCATAAGCGCCGTTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAt             >  4:1334359/4‑149 (MQ=255)
                                                                                                                                          aaacggcgCATAAGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAttt           >  1:1245788/4‑149 (MQ=255)
                                                                                                                                          aaacggcgCATAAGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGAATCTTTGCGTCCGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCTTGACTAAAAACCCTTCGCTCAATTGACAGTTTAttt           >  1:8696/4‑149 (MQ=255)
                                                                                                                                            acggcgCATAAGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTtgt                               >  2:1484640/2‑127 (MQ=255)
                                                                                                                                            acggcgCATAAGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTtgt                               <  1:1484640/126‑1 (MQ=255)
                                                                                                                                             cggcgCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCt        <  2:252995/149‑1 (MQ=255)
                                                                                                                                              ggcgCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTg       >  3:1055378/1‑149 (MQ=255)
                                                                                                                                                   ttATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGag  >  3:1492105/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
AAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAG  >  NC_000913/1369879‑1370174

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