Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,370,026 T→A 15.5% 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,0260TA15.5% 287.0 / 35.1 115intergenic (+23/‑190)pspE/ycjMthiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein
Reads supporting (aligned to +/- strand):  ref base T (44/53);  new base A (5/13);  total (49/66)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 117 total aligned reads displayed.

TAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTT  >  NC_000913/1369928‑1370123
                                                                                                  |                                                                                                    
tAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCtt                                                                                                     >  2:2475821/1‑100 (MQ=255)
tAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATCCACCATGAGCGGCGCtt                                                                                                     >  2:1120877/1‑100 (MQ=255)
 aGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTa                                                                                                    <  2:1878262/100‑1 (MQ=255)
  gCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTa                                                                                                    >  2:2320103/1‑99 (MQ=255)
  gCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTAt                                                                                                   <  1:1774199/100‑1 (MQ=255)
   cGAGATGGGATATACCCACGTTGAGAACGCCGGTTGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTa                                                                                                    <  2:1814831/98‑1 (MQ=255)
     agaTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATgcg                                                                                                >  1:386747/1‑100 (MQ=255)
      gaTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATgcgc                                                                                               <  2:46538/100‑1 (MQ=255)
      gaTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAgcgc                                                                                               <  2:1265069/100‑1 (MQ=255)
       aTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGcc                                                                                              <  2:604172/100‑1 (MQ=255)
          ggATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgct                                                                                           <  2:1532570/100‑3 (MQ=255)
             tataCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcat                                                                                        <  1:339785/100‑6 (MQ=255)
              ataCCCACGTTGAGACCGCCGGTGGCTTAAAAGACATCGCAATGCCGAACGTTAAAGGTTATAAGAATTCACCATGATCGGCGCATAAGCGCCgctcatg                                                                                       <  1:432116/100‑7 (MQ=255)
               taCCCACGTTGAGAACGCCGGTGGCTTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatgg                                                                                      <  1:611477/100‑8 (MQ=255)
               taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑ttttttt                                                                                    <  1:1640787/100‑1 (MQ=255)
               taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑ttttttt                                                                                    >  1:362853/1‑100 (MQ=255)
               taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑ttttttt                                                                                    <  1:2398169/100‑1 (MQ=255)
               taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatgg                                                                                      <  1:873061/100‑8 (MQ=255)
               taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatgg                                                                                      <  1:128414/100‑8 (MQ=255)
               gaCCCACGTTGAGAACGCCAGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatgg                                                                                      <  1:1221330/99‑8 (MQ=255)
                aCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGGCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatg                                                                                       <  1:906252/98‑7 (MQ=255)
                 ccccCGTTGCGCACGCCCGTGGCCTGAAAGAGATCGCACTGCCGAAGGTCAAAGGTTACAAGTATTCTCCATGAGCGGCGCTTATGCGCCG‑‑ttttttt                                                                                    <  1:1543169/98‑1 (MQ=255)
                 cccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑ttttttttt                                                                                  >  2:2268797/1‑100 (MQ=255)
                 cccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑ttttttttt                                                                                  >  2:1987933/1‑100 (MQ=255)
                  ccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTc                                                                                 <  1:921100/100‑1 (MQ=255)
                  ccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTc                                                                                 <  1:1873840/100‑1 (MQ=255)
                  ccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTc                                                                                 <  2:1655176/100‑1 (MQ=255)
                   cACGTTGAGAACGCCGGTGGCCTGAAAGACAGCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatg                                                                                       <  1:959519/95‑7 (MQ=255)
                    aCGTTGAGAACGGCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTg                                                                              <  2:1517198/100‑1 (MQ=255)
                     cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGt                                                                             >  2:2195859/1‑100 (MQ=255)
                     cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGt                                                                             >  1:338984/1‑100 (MQ=255)
                        tGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTAGGCGCCG‑‑TTTTGTTTTCTGTcac                                                                          <  1:1078494/100‑1 (MQ=255)
                        tGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTACGCGCCG‑‑TTTTTTTTTCTGTcac                                                                          <  1:1500457/100‑1 (MQ=255)
                          agaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCacac                                                                        <  2:1812181/100‑1 (MQ=255)
                          agaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCacac                                                                        <  2:2238525/100‑1 (MQ=255)
                           gaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACAcc                                                                       >  2:381406/1‑100 (MQ=255)
                           gaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACAcc                                                                       >  2:1881509/1‑100 (MQ=255)
                              cGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTTCTGTCACACCt                                                                      >  2:2278535/1‑100 (MQ=255)
                                ccGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTAtt                                                                  >  1:421858/1‑100 (MQ=255)
                                  ggTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTa                                                                <  2:2502267/100‑1 (MQ=255)
                                   gTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTAc                                                               <  1:1082615/100‑1 (MQ=255)
                                   gTGGCCTAAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTAttt                                                                 <  1:245147/98‑1 (MQ=255)
                                    tGGCCTGAAAGACATCGCAATGCGGAAGGTCAAAGGTTAAAAGTATTCACCATGAGCGGTGCTTATGCGCCT‑‑TTTTTTTTTCTGTCACACCTTATTTACa                                                              <  1:870693/100‑1 (MQ=255)
                                    tGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACa                                                              <  2:1936747/100‑1 (MQ=255)
                                      gCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACCCTTTATTTAc                                                               >  2:23782/1‑98 (MQ=255)
                                      gCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATc                                                            <  2:1301856/100‑1 (MQ=255)
                                       ccTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTTCTGTCACCCCtattttccc                                                              >  1:1084573/1‑91 (MQ=255)
                                       ccTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATcc                                                           >  2:1375533/1‑100 (MQ=255)
                                       ccTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATcc                                                           >  2:804655/1‑100 (MQ=255)
                                       ccTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATcc                                                           >  2:1247339/1‑100 (MQ=255)
                                        cTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATccc                                                          >  2:1843445/1‑100 (MQ=255)
                                          gAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTTTTTCATTTTATTTACATc                                                            >  1:2392254/1‑96 (MQ=255)
                                             aGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGa                                                    <  2:1635251/100‑1 (MQ=255)
                                                cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGatt                                                  >  2:589435/1‑100 (MQ=255)
                                                cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGatt                                                  >  2:589515/1‑100 (MQ=255)
                                                 aTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGatt                                                  >  1:2163423/1‑100 (MQ=255)
                                                 aTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGatta                                                 <  1:1417129/100‑1 (MQ=255)
                                                    gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTAttt                                              <  2:711365/100‑1 (MQ=255)
                                                      aaTGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCCTAGATTATTTGc                                            >  1:2138829/1‑100 (MQ=255)
                                                          ccGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTTTTTTTTCCAGTATTTACATCCCCATAAATTTTTTTCGtta                                        >  1:1184393/1‑98 (MQ=255)
                                                          ccGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCa                                        <  2:1962877/100‑1 (MQ=255)
                                                               ggTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTca                                   <  2:570153/100‑1 (MQ=255)
                                                                  cAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCCTAGATTATTTGCGTCAGCTcaca                                 >  1:1956800/1‑99 (MQ=255)
                                                                  cAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaa                                >  2:164169/1‑100 (MQ=255)
                                                                   aaaGGTTAAAAGAATTCACCATGAGCTGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaaa                               <  2:1261556/100‑1 (MQ=255)
                                                                     aGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATa                             <  1:311311/100‑1 (MQ=255)
                                                                       gTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACg                           <  1:95249/100‑1 (MQ=255)
                                                                        ttAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGc                          >  2:778303/1‑100 (MQ=255)
                                                                         tAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTTTTTTTTTCCTTATTAACATCCCCATAGATTATTTGCGTTACCTTACAAATACGCt                         >  1:2227880/1‑100 (MQ=255)
                                                                         tAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCt                         <  2:1490302/100‑1 (MQ=255)
                                                                          aaaaaaATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCtt                        <  2:125727/100‑1 (MQ=255)
                                                                           aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTTTGGTCCCACCTTATTTAAATCCCCAAAGATTATTTCCTCAAGCTCACAAATACGCtt                        >  2:417618/1‑100 (MQ=255)
                                                                           aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTTTTTTTTTTTTTATTTACATCCCCATAGATTATTTTCGTCACCTCACAAATACGCttt                       >  1:1054163/1‑100 (MQ=255)
                                                                           aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTTACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCttt                       >  2:2216819/1‑100 (MQ=255)
                                                                           aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCttt                       >  1:1475792/1‑100 (MQ=255)
                                                                           aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCttt                       >  2:360516/1‑100 (MQ=255)
                                                                                aTTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTTTTTTTCCCCTTATTTACATCCCCATAGATTTTTTACGTTACCTCCGAAATACGCTTTTccc                   >  1:195845/1‑99 (MQ=255)
                                                                                aTTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTCTCACACCTTATTTACATCCCCATACATCATTTGCGTCACACCACCAATACGCgtttttac                  >  1:1931270/1‑97 (MQ=255)
                                                                                 ttCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCt                 <  1:2515662/100‑1 (MQ=255)
                                                                                 ttCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCACt                 >  1:1788359/1‑100 (MQ=255)
                                                                                 ttCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTg                <  2:265728/100‑1 (MQ=255)
                                                                                  tCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCt                 <  2:386749/99‑1 (MQ=255)
                                                                                  tCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTg                <  2:2564878/100‑1 (MQ=255)
                                                                                  gaaaaaaaaacggcgCATAAGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACCAATACGCTTTTTCCCTg                >  2:1639324/11‑100 (MQ=38)
                                                                                  gaaaaaaaaacggcgCATAAGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTg                >  2:1652553/11‑100 (MQ=38)
                                                                                    aaaaaaaacggcgCATAAGCGCCGCTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTg                >  1:659361/9‑99 (MQ=37)
                                                                                    aCCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGt              <  1:2319206/100‑1 (MQ=255)
                                                                                    aCCATGAGCGGCCCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTgg               <  1:846700/99‑1 (MQ=255)
                                                                                     ccATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTa             >  2:709337/1‑100 (MQ=255)
                                                                                      cATGAGCGGCGCTTATGTGCCG‑‑TTTTTATTTATGTCATGCCTTGCTTACATGACGATAGATTATTTGCCTCAGTCCACAAAAACGCTTTTTCCCTGCTa             >  1:1335121/1‑99 (MQ=255)
                                                                                       aaaaacggcgCATAAGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaa           >  2:580632/6‑100 (MQ=255)
                                                                                          aGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaaaaa        >  2:1599547/1‑100 (MQ=255)
                                                                                           gcggcgCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGAGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt       <  1:773571/100‑1 (MQ=255)
                                                                                            cggcgCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATg      <  2:1362817/100‑1 (MQ=255)
                                                                                            cggcgCTTATGCGCCG‑‑GTTTTTTTTCGGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATg      <  2:815881/100‑1 (MQ=255)
                                                                                             ggcgCTTATGCGCCG‑‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAt    <  2:606662/100‑1 (MQ=255)
                                                                                              gtgcTTATGCGCCG‑‑TGTTTTTTACTGTCACACCTTATTTACATCCCCATAGATTATTTGCATCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt       <  1:720032/95‑1 (MQ=255)
                                                                                              gcgcTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGa     <  1:1754329/99‑1 (MQ=255)
                                                                                                gcTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAttt  <  1:1892576/100‑1 (MQ=255)
                                                                                                 cTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAAAAATACGCTTTTTCCCTGGTAAAAAATGAttt  >  1:515512/1‑99 (MQ=255)
                                                                                                  |                                                                                                    
TAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTT  >  NC_000913/1369928‑1370123

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: