Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,370,026 T→A 7.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 7.1% 352.0 / 10.8 125intergenic (+23/‑190)pspE/ycjMthiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein
Reads supporting (aligned to +/- strand):  ref base T (59/57);  new base A (4/5);  total (63/62)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

Only 100 of 132 total aligned reads displayed.

GAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGT  >  NC_000913/1369877‑1370175
                                                                                                                                                     |                                                                                                                                                        
gAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCt                                                                                                                                                          <  1:907378/150‑1 (MQ=255)
 aaaGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCtt                                                                                                                                                         <  1:1031645/150‑1 (MQ=255)
  aaGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTa                                                                                                                                                        >  1:1049536/1‑150 (MQ=255)
  aaGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTa                                                                                                                                                        <  1:868384/150‑1 (MQ=255)
    gtgtATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATg                                                                                                                                                      >  1:1049215/1‑150 (MQ=255)
        tttGCAATGCAGGACGCCAGTCAGGGAAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGcc                                                                                                                                                  <  2:1064450/149‑1 (MQ=255)
           gCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCAGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑tt                                                                                                                                             >  2:1137342/1‑150 (MQ=255)
              aTGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACTTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcat                                                                                                                                            <  1:1642481/150‑6 (MQ=255)
                gCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATTGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatgg                                                                                                                                          <  2:276921/150‑8 (MQ=255)
                gCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑ttttttt                                                                                                                                        <  2:1419051/150‑1 (MQ=255)
                 ccGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑tttttttt                                                                                                                                       >  2:1036883/1‑150 (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)
                     aCGCCAGTCAGGGCAAGCAAAAGAGATACTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTg                                                                                                                                   <  2:1424719/150‑1 (MQ=255)
                        ccAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTca                                                                                                                                <  1:725056/150‑1 (MQ=255)
                         cAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTcac                                                                                                                               <  2:1049536/150‑1 (MQ=255)
                             cAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCt                                                                                                                           <  1:777065/150‑1 (MQ=255)
                                ggCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACTCCGGTGGCCTGAAAGACATCGCAAAGCCGAACCTCAAAGCTTAAAAGAATTCACCATGAGCGCCGCTTATGCGCCG‑‑TTTATTTTTCAGTCACACCTTAt                                                                                                                        >  1:1046729/1‑150 (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)
                                       aaaaGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATc                                                                                                                 >  2:611735/1‑150 (MQ=255)
                                         aaGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATccc                                                                                                               <  1:776432/150‑1 (MQ=255)
                                          agagATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACAATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATccc                                                                                                               <  1:252460/150‑1 (MQ=255)
                                           gagaTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCa                                                                                                             >  2:152056/1‑150 (MQ=255)
                                             gaTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATa                                                                                                           <  2:861931/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)
                                                     gCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTATGTCACACCTTATTTACATCCCCATAGATTAttt                                                                                                   >  1:358769/1‑150 (MQ=255)
                                                        agaTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCg                                                                                                <  2:357053/150‑1 (MQ=255)
                                                          aTGGTATATACCCACGTGGAGAACGCCGGTGGCCTGAATGACATCGCAATGCCGAAGGCCAAAGGTTAAAAGATTTCACCCTGAGCGGCGGTTATGCGCCG‑‑TTTTTTTTTCTGTCAAACCTTATTTACATCCCCATAGATTATTTGAGTc                                                                                              <  2:251884/150‑1 (MQ=255)
                                                          aTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTc                                                                                              <  1:113395/150‑1 (MQ=255)
                                                          aTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTc                                                                                              <  1:617559/150‑1 (MQ=255)
                                                           tGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAATAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCa                                                                                             <  1:505458/150‑1 (MQ=255)
                                                           tGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCa                                                                                             >  2:577204/1‑150 (MQ=255)
                                                            gggATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAg                                                                                            <  1:1557816/150‑1 (MQ=255)
                                                            gggATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAg                                                                                            >  2:1358744/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: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)
                                                                              gaACGGCGGTGGCCTGAATTACATCGCAATGCCGATGGTCAAAGTTTAAAAGATTTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGAATATTTGCGTCAGCTCACAAATACGCttttt                                                                          <  1:701921/150‑1 (MQ=255)
                                                                               aaCGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTc                                                                         <  1:152056/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)
                                                                                    cGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTgg                                                                    <  1:1573084/150‑1 (MQ=255)
                                                                                         gCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaaa                                                               >  1:643167/1‑150 (MQ=255)
                                                                                          ccTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTATGTCACACATTATTTACATCCCCATAGATTTTTTGCGTCAGCTCCCAACTACGATTTTTCCCTGGTtaaa                                                               >  2:962362/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)
                                                                                                       gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACCCCTTATTTACATCCCCATAGATTATTTGCGTCAGCTAACAAATAAGCTTTTTCCCTGGTATAAAATGAGTTCCTGCg                                                 >  2:1266646/1‑150 (MQ=255)
                                                                                                       gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTATCCCTGGTAAAAAATGATTTCCTGCg                                                 >  1:424845/1‑150 (MQ=255)
                                                                                                       gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACAACTTATTTACATCCCCATAGATTATTTGCGTAAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCg                                                 >  2:66563/1‑150 (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)
                                                                                                                     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)
                                                                                                                                 gAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTg                       <  2:1685185/150‑1 (MQ=255)
                                                                                                                                  aaTTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGAc                     >  2:1218306/1‑150 (MQ=255)
                                                                                                                                   aTTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATATATTATTTGCGTCAGCTAACAAATAAGCTTTTTACCTGGTAAAAAATTATTTGCTGCGTGACTAAAACCCTTGTGCTCAATTGAc                     >  2:1512714/1‑150 (MQ=255)
                                                                                                                                   aTTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTAAGTCACACCTTATTTTCATCCCCATAGATTATTTGCGTCAGCTCAAAAAATCGCATTTTGCCATGTAAAAAAAGAATTCCTCCGTGAAAAAAACCCTTGTGCTAAATTAACa                    >  2:175577/1‑150 (MQ=255)
                                                                                                                                      cACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCAGTCCCGCCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTTAAAAATGATTTCCTGCGTGACAAAAACCCATGTGCTCACTTGACAGtt                 >  2:349178/1‑150 (MQ=255)
                                                                                                                                       aCCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGGTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACGAAAACCCTTGTGGGCAATTGACAGttt                >  2:313925/1‑150 (MQ=255)
                                                                                                                                        ccATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTGACACCTTATTTACATCCCAATAGATTATTTTCGTCAGCTCTCAAATACGCTTTTTACCTGGTAAAAAATGATTTCCTGGGTGACTAAAACCCTAGTGCGCAATTGACAGTTTa               >  1:1103591/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)
                                                                                                                                          aaaaacggcgCATAAGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAtt             >  1:1247361/6‑150 (MQ=255)
                                                                                                                                           aaaacggcgCATAAGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTACCTGGTAAAAAATGATTTACTGCGTGAATAAAACCCGTGTGCGCAATTGACAGGTTAttt            >  1:1573918/5‑150 (MQ=255)
                                                                                                                                             aGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTACTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTAttttg          >  2:516856/1‑149 (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)
                                                                                                                                                gtagattATGAGCCT‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGATCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGc       <  2:158272/145‑1 (MQ=255)
                                                                                                                                                  cgcTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTTTTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCgg     <  1:1200802/150‑1 (MQ=255)
                                                                                                                                                     ttATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGagt  <  1:519341/150‑1 (MQ=255)
                                                                                                                                                     ttATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGagt  <  1:611735/150‑1 (MQ=255)
                                                                                                                                                     ttATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTAGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGagt  <  1:1020008/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                        
GAAAGTGTATTGCAATGCCGGACGCCAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTGCGGAGT  >  NC_000913/1369877‑1370175

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: