Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,370,029 T→A 7.0% 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 7.0% 374.3 / 12.0 128intergenic (+26/‑187)pspE/ycjMthiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein
Reads supporting (aligned to +/- strand):  ref base T (59/60);  new base A (4/5);  total (63/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 = 9.58e-01

Only 100 of 129 total aligned reads displayed.

GTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGT  >  NC_000913/1369881‑1370178
                                                                                                                                                    |                                                                                                                                                        
gtgtATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATg                                                                                                                                                         >  1:1049215/1‑150 (MQ=255)
    tttGCAATGCAGGACGCCAGTCAGGGAAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGcc                                                                                                                                                     <  2:1064450/149‑1 (MQ=255)
     ttGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGATCGGCTCTTATGCGACg                                                                                                                                                    >  1:158272/1‑150 (MQ=255)
      tGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgc                                                                                                                                                   <  1:618008/150‑2 (MQ=255)
           tGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑tttttt                                                                                                                                            >  1:1552115/1‑150 (MQ=255)
            gCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATTGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatgg                                                                                                                                             <  2:276921/150‑8 (MQ=255)
            gCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACTATGAGCGGCGCTTATGCGCCG‑‑ttttttt                                                                                                                                           >  2:519341/1‑150 (MQ=255)
            gCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑ttttttt                                                                                                                                           <  2:1419051/150‑1 (MQ=255)
            gCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatgg                                                                                                                                             <  1:734893/150‑8 (MQ=255)
               ggACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑ttttttttta                                                                                                                                        >  1:1391509/1‑149 (MQ=255)
               ggACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTc                                                                                                                                        <  1:1364368/150‑1 (MQ=255)
                gACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCt                                                                                                                                       >  2:10940/1‑150 (MQ=255)
                 aCGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTg                                                                                                                                      <  2:198111/150‑1 (MQ=255)
                    ccAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTca                                                                                                                                   <  1:725056/150‑1 (MQ=255)
                     cAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTcac                                                                                                                                  <  2:1049536/150‑1 (MQ=255)
                              caagcaaATGCGTTCCTTAGCGAGATGGGATATAGCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGTCGAAGGTCAAAGGTTAAAAGAATTAACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTAttt                                                                                                                         <  1:141777/150‑1 (MQ=255)
                              caagcaaAAGAGATCCTTAGCGTGATGGGATATACCCACGTTGAGAACGCCGCTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGACCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTAttt                                                                                                                         <  2:480925/150‑1 (MQ=255)
                                 gcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACa                                                                                                                      >  2:193295/1‑150 (MQ=255)
                                  caaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACAt                                                                                                                     <  1:419485/150‑1 (MQ=255)
                                   aaaaGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATc                                                                                                                    >  2:1200802/1‑150 (MQ=255)
                                   aaaaGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATc                                                                                                                    >  2:611735/1‑150 (MQ=255)
                                   aaaaGAGATCCTTAGCGAGATGGGATATACACACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCATCATGAGCGGCGCTGATGCGCCG‑‑TTTTTTTTTCTGTCACAACTTAGTTAAATc                                                                                                                    >  2:1505531/1‑150 (MQ=255)
                                     aaGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATccc                                                                                                                  <  1:776432/150‑1 (MQ=255)
                                       gagaTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCa                                                                                                                >  2:152056/1‑150 (MQ=255)
                                       gagaTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCa                                                                                                                >  1:507576/1‑150 (MQ=255)
                                         gaTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATa                                                                                                              <  2:861931/150‑1 (MQ=255)
                                          aTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAg                                                                                                             <  2:1258599/150‑1 (MQ=255)
                                           tCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGa                                                                                                            >  1:915116/1‑150 (MQ=255)
                                           tCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACAACTTATTTACATCCCCATAGa                                                                                                            >  1:1141369/1‑150 (MQ=255)
                                                aGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTAAGTAACACCTTATTAACAACCCCATAGattaat                                                                                                       >  2:1509726/1‑148 (MQ=255)
                                                    agaTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCg                                                                                                   <  2:357053/150‑1 (MQ=255)
                                                      aTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTc                                                                                                 <  1:617559/150‑1 (MQ=255)
                                                      aTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTc                                                                                                 <  1:113395/150‑1 (MQ=255)
                                                       tGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAATAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCa                                                                                                <  1:505458/150‑1 (MQ=255)
                                                       tGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCa                                                                                                >  2:577204/1‑150 (MQ=255)
                                                              taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTcaca                                                                                         <  1:1714460/150‑1 (MQ=255)
                                                               aCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaa                                                                                        <  2:1401385/150‑1 (MQ=255)
                                                                cccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATTAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaaa                                                                                       <  1:1034994/150‑1 (MQ=255)
                                                                  cACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAAt                                                                                      >  1:207938/1‑150 (MQ=255)
                                                                  cACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAAt                                                                                      >  1:1685185/1‑150 (MQ=255)
                                                                   aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATa                                                                                     >  2:522485/1‑150 (MQ=255)
                                                                   aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATAc                                                                                    >  1:387592/1‑150 (MQ=255)
                                                                        gagaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCAAAGATTATTTGCGTCAGCTCACAAATAAGCtt                                                                                >  1:907834/1‑150 (MQ=255)
                                                                         agaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCATCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCtttt                                                                              <  1:1505531/150‑1 (MQ=255)
                                                                            aCGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTcc                                                                           <  2:299256/150‑1 (MQ=255)
                                                                             cGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTccc                                                                          <  2:1255492/150‑1 (MQ=255)
                                                                              gCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCt                                                                         <  2:215267/150‑1 (MQ=255)
                                                                               ccGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTg                                                                        <  2:572212/150‑1 (MQ=255)
                                                                               ccGGTGGCCTGAAAGACATAGCATTGCCGAAGGTCAAAGTTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTg                                                                        <  2:1649063/150‑1 (MQ=255)
                                                                                cGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTgg                                                                       <  1:1573084/150‑1 (MQ=255)
                                                                                     gCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaaa                                                                  >  1:643167/1‑150 (MQ=255)
                                                                                      ccTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaaa                                                                  <  1:308284/150‑1 (MQ=255)
                                                                                      ccTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTATGTCACACATTATTTACATCCCCATAGATTTTTTGCGTCAGCTCCCAACTACGATTTTTCCCTGGTtaaa                                                                  >  2:962362/1‑150 (MQ=255)
                                                                                      ccTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTTTTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTACCTGGTaaaaa                                                                 >  1:1089927/1‑150 (MQ=255)
                                                                                      ccTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaaaa                                                                 >  2:1491735/1‑150 (MQ=255)
                                                                                       cTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATAATTTGCGTCAGCTCACAAATAAGCGTTTTCCCTGGTaaaaaa                                                                >  1:1276078/1‑150 (MQ=255)
                                                                                             gACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAttt                                                          <  2:11014/150‑1 (MQ=255)
                                                                                             gACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAttt                                                          >  2:1188441/1‑150 (MQ=255)
                                                                                                 tCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTg                                                      <  2:438822/150‑1 (MQ=255)
                                                                                                 tCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTg                                                      <  2:1552115/150‑1 (MQ=255)
                                                                                                 tCGCAATGCCGAAGGTCAAAGGTTAAAACAATTCACAATGAGCGGCGCTTTTGCGCAG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTg                                                      <  1:1137342/150‑1 (MQ=255)
                                                                                                   gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACCCCTTATTTACATCCCCATAGATTATTTGCGTCAGCTAACAAATAAGCTTTTTCCCTGGTATAAAATGAGTTCCTGCg                                                    >  2:1266646/1‑150 (MQ=255)
                                                                                                   gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACAACTTATTTACATCCCCATAGATTATTTGCGTAAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCg                                                    >  2:66563/1‑150 (MQ=255)
                                                                                                        gCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACt                                              <  2:334292/150‑1 (MQ=255)
                                                                                                           gAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTaaa                                           <  1:1013540/150‑1 (MQ=255)
                                                                                                              gggcAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAAcc                                        <  2:821420/147‑1 (MQ=255)
                                                                                                              ggTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTACCTGGTAAAAAATGATTTCCTGCGTGACTAAAAcc                                        >  1:336906/1‑150 (MQ=255)
                                                                                                               gTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACAACTTATTTACATACCCATAGATTATTTGCGTAAGCGAACATATACGCTTTTTCCCTGGTAAAAAATGATTTAATGCGTGACTAAAAcc                                        >  1:270572/1‑150 (MQ=255)
                                                                                                               gTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAAccc                                       >  1:420505/1‑150 (MQ=255)
                                                                                                                tCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCt                                      <  1:1036883/150‑1 (MQ=255)
                                                                                                                 cAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCtt                                     >  2:1219472/1‑150 (MQ=255)
                                                                                                                  aaaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTtg                                    <  1:10940/150‑1 (MQ=255)
                                                                                                                  aaaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTATGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTGACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTTAAACCCTtg                                    >  2:225731/1‑150 (MQ=255)
                                                                                                                   aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTATGTCACACCTTATTTACATCCCCAAAGATTATTTGCGTAAGCTCCCCAATACGCTTTTTACCTGGTAAAAAATGATTTACTGCGTGACTAAAACCCTtgt                                   >  1:1286924/1‑150 (MQ=255)
                                                                                                                    aGGTTAAAAGAATTCACCATGTGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTtgtg                                  <  1:1521165/150‑1 (MQ=255)
                                                                                                                     ggTTAAAAAATTTCCCAATTAGCGGCGCTTATTCGCGT‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGTTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGc                                 <  2:109275/150‑1 (MQ=255)
                                                                                                                         aaaaGAATTCACCATGAGCGGCGCTTATGCGCCC‑‑TTTTTTTTTCAGTCACACCTTATTTTCATCCCCATAGATAATTTGGCTGAGCTCACAAAAACGCTATATAACTGTTAAAAAATGATTAACTGCGTGACTAAAACCCATCTGCGCaa                             >  2:721986/1‑150 (MQ=255)
                                                                                                                          aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCGGGTAAAAAATGATTTCCTGCGTGACTAAAAACCTTGTGCTCAAt                            >  2:1141348/1‑150 (MQ=255)
                                                                                                                           aattaTTCACCATGAGATGCGCTTATTGGCCG‑‑TTTTTTTTTCTGTCACACCTTATATACATCCACATAGATTATTTGCGTCAGTTCACAAATACGGTTTTTCCCTGGTTAAATATGATTTCCTGCTTGTCTTAAACCCTTGTGCTCAAtt                           <  1:293202/146‑1 (MQ=255)
                                                                                                                           aaGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAtt                           <  2:1057663/150‑1 (MQ=255)
                                                                                                                            aGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTg                          <  2:1018636/150‑1 (MQ=255)
                                                                                                                             gAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTg                          <  2:1685185/150‑1 (MQ=255)
                                                                                                                               aTTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACa                       <  1:874137/150‑1 (MQ=255)
                                                                                                                               aTTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACa                       >  2:140377/1‑150 (MQ=255)
                                                                                                                                 tCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATTGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGt                     >  2:1133758/1‑150 (MQ=255)
                                                                                                                                  cACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCAGTCCCGCCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTTAAAAATGATTTCCTGCGTGACAAAAACCCATGTGCTCACTTGACAGtt                    >  2:349178/1‑150 (MQ=255)
                                                                                                                                    aaaaaaacggcgCATAAGCGCCGCTTTTTTTTTTATGTAACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCACACAAATACGCATTTTCCCTGGTAAAAAATTATTTCATGGGTGACTAAAACCCTTGTGATCAATTGACAGtt                    >  2:344770/8‑150 (MQ=255)
                                                                                                                                     cATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATAAGCTTTTTACCTGGTAAAAAATGATTTACTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAt                 >  2:1553469/1‑150 (MQ=255)
                                                                                                                                     cATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTAAATCCCCATATATTATTTGCGTCAGCTCACAAATACGCTATTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCACAATTGACAGTTTAt                 >  2:1560272/1‑150 (MQ=255)
                                                                                                                                     aaaaaacggcgCATAAGCGCCG‑‑TTTTTTTTTGTGTAACACCTTATTTAAATACCCATAGATTATTTGCGTGAGCTCACAAATAAGCTTTTTACCTGGTAAAAAATGATTTCCTGCGTGACTTAAACCCCGTGTGCGCAATTGACAGtttt                  >  2:1702904/7‑149 (MQ=255)
                                                                                                                                     aaTGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACAACTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAt                 <  2:269182/149‑1 (MQ=255)
                                                                                                                                       aaaacggcgCATAAGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTACCTGGTAAAAAATGATTTACTGCGTGAATAAAACCCGTGTGCGCAATTGACAGGTTAttt               >  1:1573918/5‑150 (MQ=255)
                                                                                                                                         aGCGGCGATTATGCTCCG‑‑TTTTTTTTTATGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAATCCCTTGTGCTCAATTGACAGTTTATTTTc             <  1:127688/150‑1 (MQ=255)
                                                                                                                                          gcggcgCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAttttat            >  1:1607776/1‑148 (MQ=255)
                                                                                                                                           cggcgCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTAATTTACATCCCCATAGATTATTTTCGTCATCTCACAAATACGCTTTTTAACTGGTTAAAAATGAATTACTGCGTGAATAAAACCCGTGTGCTCAATTGACAGTTTATTTTCTg           >  1:1520644/1‑150 (MQ=255)
                                                                                                                                              cgcTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTTTTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCgg        <  1:1200802/150‑1 (MQ=255)
                                                                                                                                                 ttATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGagt     <  1:611735/150‑1 (MQ=255)
                                                                                                                                                 ttATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGagt     <  1:1020008/150‑1 (MQ=255)
                                                                                                                                                  tATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTACCTGGTAAAAAATGATTTACTGCGTGACTAAAACCCTTGTGCTAAATTGACAGTTTATTTTCTGCGGagta    >  1:447650/1‑150 (MQ=255)
                                                                                                                                                    tGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCTGagtagt  >  2:417369/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                        
GTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGTAGT  >  NC_000913/1369881‑1370178

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: