Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,370,029 T→A 16.5% 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,0290TA16.5% 261.5 / 36.1 110intergenic (+26/‑187)pspE/ycjMthiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein
Reads supporting (aligned to +/- strand):  ref base T (43/47);  new base A (5/13);  total (48/62)
Fisher's exact test for biased strand distribution p-value = 1.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

GCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGC  >  NC_000913/1369930‑1370128
                                                                                                   |                                                                                                      
gCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTAt                                                                                                        <  1:1774199/100‑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)
       gggATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgc                                                                                                 <  1:1929571/100‑2 (MQ=255)
        ggATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgct                                                                                                <  2:1532570/100‑3 (MQ=255)
           tataCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcat                                                                                             <  1:339785/100‑6 (MQ=255)
            ataCCCACGTTGAGACCGCCGGTGGCTTAAAAGACATCGCAATGCCGAACGTTAAAGGTTATAAGAATTCACCATGATCGGCGCATAAGCGCCgctcatg                                                                                            <  1:432116/100‑7 (MQ=255)
            ataCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatg                                                                                            <  2:610211/100‑7 (MQ=255)
             tagccACGTTGAGAACCCCGGTGGCCTGAAAGCCATCCCACTGCCGAAGGTCAAAGGTTAACAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatgg                                                                                           <  1:325444/97‑8 (MQ=255)
             taCCCACGTTGAGAACGCCGGTGGCTTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatgg                                                                                           <  1:611477/100‑8 (MQ=255)
             taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑ttttttt                                                                                         <  1:1640787/100‑1 (MQ=255)
             taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑ttttttt                                                                                         <  1:2398169/100‑1 (MQ=255)
             taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatgg                                                                                           <  1:128414/100‑8 (MQ=255)
             taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatgg                                                                                           <  1:873061/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                                                                                      <  2:1655176/100‑1 (MQ=255)
                 cACGTTGAGAACGCCGGTGGCCTGAAAGACAGCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatg                                                                                            <  1:959519/95‑7 (MQ=255)
                   cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGt                                                                                  >  1:338984/1‑100 (MQ=255)
                   cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGt                                                                                  >  2:2195859/1‑100 (MQ=255)
                      tGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTAGGCGCCG‑‑TTTTGTTTTCTGTcac                                                                               <  1:1078494/100‑1 (MQ=255)
                      tGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTACGCGCCG‑‑TTTTTTTTTCTGTcac                                                                               <  1:1500457/100‑1 (MQ=255)
                        agaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCacac                                                                             <  2:2238525/100‑1 (MQ=255)
                        agaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCacac                                                                             <  2:1812181/100‑1 (MQ=255)
                         gaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACAcc                                                                            >  2:1881509/1‑100 (MQ=255)
                         gaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACAcc                                                                            >  2:1281157/1‑100 (MQ=255)
                         gaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACAcc                                                                            >  2:381406/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:1937473/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:1247339/1‑100 (MQ=255)
                                     ccTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATcc                                                                >  2:804655/1‑100 (MQ=255)
                                     ccTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATcc                                                                >  2:1375533/1‑100 (MQ=255)
                                      cTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATccc                                                               >  2:1843445/1‑100 (MQ=255)
                                      cTGAAAGACATCGCAACGCCGAAGTTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATccc                                                               <  2:467048/100‑1 (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)
                                               aTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGatt                                                       >  1:2163423/1‑100 (MQ=255)
                                                 cGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGattatt                                                    >  1:463400/1‑100 (MQ=255)
                                                  gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTAttt                                                   <  2:711365/100‑1 (MQ=255)
                                                   cAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTg                                                  >  2:262309/1‑100 (MQ=255)
                                                    aaTGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCCTAGATTATTTGc                                                 >  1:2138829/1‑100 (MQ=255)
                                                     aTGCCGAAGGTAAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCg                                                <  2:1543531/100‑1 (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‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCttt                            >  1:1475792/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:2269117/100‑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)
                                                                                 cACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATATATTATTTGCGTCAGCTCACAAATATGCTTTTTCCCTgg                    <  1:1602446/100‑1 (MQ=255)
                                                                                  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)
                                                                                          cggcgCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATg           <  2:1362817/100‑1 (MQ=255)
                                                                                          cggcgCTTATGCGCCG‑‑GTTTTTTTTCGGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATg           <  2:815881/100‑1 (MQ=255)
                                                                                           ggcgCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATg           <  1:858127/99‑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)
                                                                                                   tGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGc  >  1:1305447/1‑100 (MQ=255)
                                                                                                   |                                                                                                      
GCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGC  >  NC_000913/1369930‑1370128

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: