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‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTT  >  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‑‑TTTTTTTTTCT‑g                                                                              <  2:1517198/100‑1 (MQ=255)
                     cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑Gt                                                                             >  2:2195859/1‑100 (MQ=255)
                     cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑Gt                                                                             >  1:338984/1‑100 (MQ=255)
                        tGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTAGGCGCCG‑‑TTTTGTTTTCT‑GTcac                                                                          <  1:1078494/100‑1 (MQ=255)
                        tGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTACGCGCCG‑‑TTTTTTTTTCT‑GTcac                                                                          <  1:1500457/100‑1 (MQ=255)
                          agaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCacac                                                                        <  2:1812181/100‑1 (MQ=255)
                          agaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCacac                                                                        <  2:2238525/100‑1 (MQ=255)
                           gaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACAcc                                                                       >  2:381406/1‑100 (MQ=255)
                           gaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACAcc                                                                       >  2:1881509/1‑100 (MQ=255)
                              cGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTTCT‑GTCACACCt                                                                      >  2:2278535/1‑100 (MQ=255)
                                ccGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTAtt                                                                  >  1:421858/1‑100 (MQ=255)
                                  ggTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTa                                                                <  2:2502267/100‑1 (MQ=255)
                                   gTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTAc                                                               <  1:1082615/100‑1 (MQ=255)
                                   gTGGCCTAAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTAttt                                                                 <  1:245147/98‑1 (MQ=255)
                                    tGGCCTGAAAGACATCGCAATGCGGAAGGTCAAAGGTTAAAAGTATTCACCATGAGCGGTGCTTATGCGCCT‑‑TTTTTTTTTCT‑GTCACACCTTATTTACa                                                              <  1:870693/100‑1 (MQ=255)
                                    tGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACa                                                              <  2:1936747/100‑1 (MQ=255)
                                      gCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGT‑TTTTTTTTTCT‑GTCACCCTTTATTTAc                                                               >  2:23782/1‑98 (MQ=255)
                                      gCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATc                                                            <  2:1301856/100‑1 (MQ=255)
                                       ccTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTTCT‑GTCACCCCtattttccc                                                              >  1:1084573/1‑91 (MQ=255)
                                       ccTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATcc                                                           >  2:1375533/1‑100 (MQ=255)
                                       ccTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATcc                                                           >  2:804655/1‑100 (MQ=255)
                                       ccTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATcc                                                           >  2:1247339/1‑100 (MQ=255)
                                        cTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATccc                                                          >  2:1843445/1‑100 (MQ=255)
                                          gAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑TTTTCATTTTATTTACATc                                                            >  1:2392254/1‑96 (MQ=255)
                                             aGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑‑TTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGa                                                    <  2:1635251/100‑1 (MQ=255)
                                                cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGatt                                                  >  2:589435/1‑100 (MQ=255)
                                                cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGatt                                                  >  2:589515/1‑100 (MQ=255)
                                                 aTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGT‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGatt                                                  >  1:2163423/1‑100 (MQ=255)
                                                 aTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGatta                                                 <  1:1417129/100‑1 (MQ=255)
                                                    gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTAttt                                              <  2:711365/100‑1 (MQ=255)
                                                      aaTGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCCTAGATTATTTGc                                            >  1:2138829/1‑100 (MQ=255)
                                                          ccGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTTT‑TTTTTCCAGTATTTACATCCCCATAAATTTTTTTCGtta                                        >  1:1184393/1‑98 (MQ=255)
                                                          ccGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCa                                        <  2:1962877/100‑1 (MQ=255)
                                                               ggTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTca                                   <  2:570153/100‑1 (MQ=255)
                                                                  cAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCCTAGATTATTTGCGTCAGCTcaca                                 >  1:1956800/1‑99 (MQ=255)
                                                                  cAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaa                                >  2:164169/1‑100 (MQ=255)
                                                                   aaaGGTTAAAAGAATTCACCATGAGCTGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaaa                               <  2:1261556/100‑1 (MQ=255)
                                                                     aGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATa                             <  1:311311/100‑1 (MQ=255)
                                                                       gTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACg                           <  1:95249/100‑1 (MQ=255)
                                                                        ttAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGc                          >  2:778303/1‑100 (MQ=255)
                                                                         tAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTTT‑TTTTTTCCTTATTAACATCCCCATAGATTATTTGCGTTACCTTACAAATACGCt                         >  1:2227880/1‑100 (MQ=255)
                                                                         tAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCt                         <  2:1490302/100‑1 (MQ=255)
                                                                          aaaaaaATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCtt                        <  2:125727/100‑1 (MQ=255)
                                                                           aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTTTGGTCCCACCTTATTTAAATCCCCAAAGATTATTTCCTCAAGCTCACAAATACGCtt                        >  2:417618/1‑100 (MQ=255)
                                                                           aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTTT‑TTTTTTTTTTATTTACATCCCCATAGATTATTTTCGTCACCTCACAAATACGCttt                       >  1:1054163/1‑100 (MQ=255)
                                                                           aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTTACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCttt                       >  2:2216819/1‑100 (MQ=255)
                                                                           aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCttt                       >  1:1475792/1‑100 (MQ=255)
                                                                           aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCttt                       >  2:360516/1‑100 (MQ=255)
                                                                                aTTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTTT‑TTTTCCCCTTATTTACATCCCCATAGATTTTTTACGTTACCTCCGAAATACGCTTTTccc                   >  1:195845/1‑99 (MQ=255)
                                                                                aTTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑CTCACACCTTATTTACATCCCCATACATCATTTGCGTCACACCACCAATACGCgtttttac                  >  1:1931270/1‑97 (MQ=255)
                                                                                 ttCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCt                 <  1:2515662/100‑1 (MQ=255)
                                                                                 ttCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCACt                 >  1:1788359/1‑100 (MQ=255)
                                                                                 ttCACCATGAGCGGCGCTTATGCGCCG‑‑‑TTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTg                <  2:265728/100‑1 (MQ=255)
                                                                                  tCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCt                 <  2:386749/99‑1 (MQ=255)
                                                                                  tCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTg                <  2:2564878/100‑1 (MQ=255)
                                                                                  gaaaaaaaaacggcgCATAAGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACCAATACGCTTTTTCCCTg                >  2:1639324/11‑100 (MQ=38)
                                                                                  gaaaaaaaaacggcgCATAAGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTg                >  2:1652553/11‑100 (MQ=38)
                                                                                    aaaaaaaacggcgCATAAGCGCCGC‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTg                >  1:659361/9‑99 (MQ=37)
                                                                                    aCCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGt              <  1:2319206/100‑1 (MQ=255)
                                                                                    aCCATGAGCGGCCCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTgg               <  1:846700/99‑1 (MQ=255)
                                                                                     ccATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTa             >  2:709337/1‑100 (MQ=255)
                                                                                      cATGAGCGGCGCTTATGTGCCG‑‑TTTTTATTTAT‑GTCATGCCTTGCTTACATGACGATAGATTATTTGCCTCAGTCCACAAAAACGCTTTTTCCCTGCTa             >  1:1335121/1‑99 (MQ=255)
                                                                                       aaaaacggcgCATAAGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaa           >  2:580632/6‑100 (MQ=255)
                                                                                          aGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaaaaa        >  2:1599547/1‑100 (MQ=255)
                                                                                           gcggcgCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGAGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt       <  1:773571/100‑1 (MQ=255)
                                                                                            cggcgCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATg      <  2:1362817/100‑1 (MQ=255)
                                                                                            cggcgCTTATGCGCCG‑‑GTTTTTTTTCG‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATg      <  2:815881/100‑1 (MQ=255)
                                                                                             ggcgCTTATGCGCCG‑‑‑TTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAt    <  2:606662/100‑1 (MQ=255)
                                                                                              gtgcTTATGCGCCG‑‑TGTTTTTTACT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCATCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt       <  1:720032/95‑1 (MQ=255)
                                                                                              gcgcTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGa     <  1:1754329/99‑1 (MQ=255)
                                                                                                gcTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAttt  <  1:1892576/100‑1 (MQ=255)
                                                                                                 cTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCAAAAATACGCTTTTTCCCTGGTAAAAAATGAttt  >  1:515512/1‑99 (MQ=255)
                                                                                                  |                                                                                                    
TAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCT‑GTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTT  >  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: ‑