Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,370,026 T→A 15.6% 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.6% 394.0 / 62.0 160intergenic (+23/‑190)pspE/ycjMthiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein
Reads supporting (aligned to +/- strand):  ref base T (60/75);  new base A (6/19);  total (66/94)
Fisher's exact test for biased strand distribution p-value = 7.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.19e-01

Only 100 of 162 total aligned reads displayed.

GAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCC  >  NC_000913/1369932‑1370125
                                                                                              |                                                                                                    
gagaTGGGATATGCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATgc                                                                                                 >  1:2591228/1‑100 (MQ=255)
 agaTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATgcg                                                                                                >  2:282691/1‑100 (MQ=255)
  gaTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATgcgc                                                                                               >  1:964733/1‑100 (MQ=255)
  gaTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAgcgc                                                                                               <  2:1099645/100‑1 (MQ=255)
     gggATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgc                                                                                            <  1:3248898/100‑2 (MQ=255)
      ggATATATCCACGTTGAGAAAGCCGTTGGTTTGGCAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgct                                                                                           <  2:1407054/100‑3 (MQ=255)
      ggATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgct                                                                                           <  1:2808018/100‑3 (MQ=255)
      ggATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgct                                                                                           <  2:1888277/100‑3 (MQ=255)
       gATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttt                                                                                         <  1:3551888/100‑1 (MQ=255)
        atatACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctca                                                                                         <  1:2621388/100‑5 (MQ=255)
         tataCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttttt                                                                                       <  1:1868062/100‑1 (MQ=255)
         tataCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcat                                                                                        <  1:559032/100‑6 (MQ=255)
         tataCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcat                                                                                        <  1:387934/100‑6 (MQ=255)
          attcccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatg                                                                                       <  1:1894401/97‑7 (MQ=255)
          ataCCCACGTTTAGAACGCCGGTGGCTTGAACGACATCCCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAACGGCTCATAAGCGCCgctcat                                                                                        <  2:2765835/99‑6 (MQ=255)
           taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttttttt                                                                                     >  1:198239/1‑100 (MQ=255)
           taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttttttt                                                                                     <  1:1589739/100‑1 (MQ=255)
           taCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatgg                                                                                      <  2:1123595/100‑8 (MQ=255)
            aCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑tttttttt                                                                                    >  2:452067/1‑100 (MQ=255)
            aCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatggt                                                                                     <  1:3576634/100‑9 (MQ=255)
            aCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatggt                                                                                     <  1:2477799/100‑9 (MQ=255)
            aCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatggt                                                                                     <  1:1689405/100‑9 (MQ=255)
            aCCCACGTTGAGAACGCCGGTGGCCTGAAACACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatgg                                                                                      <  2:411984/99‑8 (MQ=255)
             cccACGTTGAGAACGCCGTTGGCCTGAAAGATATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCATAAGCGCCgctcatgg                                                                                      <  2:905070/98‑8 (MQ=255)
             cccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttttttttt                                                                                   >  2:3248696/1‑100 (MQ=255)
             cccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTAACCATGAGCGGCGATTATGCGCCG‑ttttttttt                                                                                   >  1:1507570/1‑100 (MQ=255)
              ccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTc                                                                                  <  2:723162/100‑1 (MQ=255)
               cACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCt                                                                                 <  1:2943357/100‑1 (MQ=255)
               cACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCt                                                                                 <  1:2056223/100‑1 (MQ=255)
                aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑ttttttttttct                                                                                >  2:1469525/1‑97 (MQ=255)
                   ttGAGAACGGCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTc                                                                              <  2:2331642/100‑1 (MQ=255)
                   ttGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTc                                                                              <  1:1446348/100‑1 (MQ=255)
                   ttGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTca                                                                             >  2:1937863/1‑100 (MQ=255)
                    tGAGAACGCGGGTGGCATGAAAGACATCGCAAGGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTGATGCGCCG‑TTTTTTTTTCTGTca                                                                             <  2:2485996/99‑1 (MQ=255)
                    tGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTcac                                                                            >  2:87742/1‑100 (MQ=255)
                    tGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTca                                                                             >  2:1231930/1‑99 (MQ=255)
                     gagaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTcaca                                                                           <  2:3702170/100‑1 (MQ=255)
                       gaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACAcc                                                                         <  1:667357/100‑1 (MQ=255)
                         aCGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCtt                                                                       >  1:538318/1‑100 (MQ=255)
                         aCGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCtt                                                                       >  2:3154355/1‑100 (MQ=255)
                          cGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTa                                                                      >  1:3426864/1‑100 (MQ=255)
                             cGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTAttt                                                                   <  1:200230/100‑1 (MQ=255)
                                tGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACa                                                                <  1:2259185/100‑1 (MQ=255)
                                tGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACa                                                                <  2:1741008/100‑1 (MQ=255)
                                  gacTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATc                                                              <  2:59565/98‑1 (MQ=255)
                                    cTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATcc                                                             <  2:1630182/100‑1 (MQ=255)
                                      gAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCa                                                          >  1:2717158/1‑100 (MQ=255)
                                       aaaGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCAt                                                         <  2:2127522/100‑1 (MQ=255)
                                         aGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTTTCACACCTTATTTACATCCCCATAg                                                       >  2:2680263/1‑100 (MQ=255)
                                           aCATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGat                                                     >  1:1023273/1‑100 (MQ=255)
                                           aCATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGat                                                     <  1:3154376/100‑1 (MQ=255)
                                            cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTTTTTTTTTTCTTATTTACATCCCTCTAGat                                                     >  2:466473/1‑99 (MQ=255)
                                            cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGatt                                                    >  1:289050/1‑100 (MQ=255)
                                            cATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGatt                                                    >  1:691012/1‑100 (MQ=255)
                                             aTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCTTAGatt                                                    >  2:567093/1‑99 (MQ=255)
                                              tCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTTTTTTTCAACTTATTTACATCCGCATAGa                                                      >  2:2346192/1‑96 (MQ=255)
                                                gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTAttt                                                >  1:1228869/1‑100 (MQ=255)
                                                gCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTAttt                                                <  2:1809552/100‑1 (MQ=255)
                                                 cAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTTTTTTTTTCCTTATTTACATCCCCATAGATTAttt                                                >  2:2401136/1‑99 (MQ=255)
                                                   aTGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCg                                             <  1:3726893/100‑1 (MQ=255)
                                                    tGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGt                                            <  2:1885855/100‑1 (MQ=255)
                                                      ccGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCa                                          <  2:523580/100‑1 (MQ=255)
                                                          aGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTc                                      <  2:213815/100‑1 (MQ=255)
                                                           ggTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCt                                       >  2:822170/1‑98 (MQ=255)
                                                            gTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCACATAGATTATTTGCGTCAGCTcac                                    >  2:334089/1‑100 (MQ=255)
                                                             tCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTcaca                                   <  1:3790728/100‑1 (MQ=255)
                                                             tCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTcaca                                   >  1:100972/1‑100 (MQ=255)
                                                                aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAAt                                >  1:2413625/1‑100 (MQ=255)
                                                                 aGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCT‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGcaca                                     <  2:2316170/94‑1 (MQ=255)
                                                                  ggTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTTTCACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAAt                                >  1:1790100/1‑99 (MQ=255)
                                                                      aaaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCtt                          >  2:3568077/1‑100 (MQ=255)
                                                                       aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCttt                         <  2:2211670/100‑1 (MQ=255)
                                                                       aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCttt                         <  2:1988243/100‑1 (MQ=255)
                                                                       aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCttt                         >  1:260072/1‑100 (MQ=255)
                                                                         aGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCttttt                       <  2:1275413/100‑1 (MQ=255)
                                                                          gAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCGGCTCACAAATACGCTTTTTc                      <  2:2413613/100‑1 (MQ=255)
                                                                           aaTTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTcc                     <  1:1261420/100‑1 (MQ=255)
                                                                           aaTTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTcc                     <  1:2227532/100‑1 (MQ=255)
                                                                            agtCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTccc                    <  2:2669676/98‑1 (MQ=255)
                                                                            aTTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTTTCCCACCATCTTTACATCCCCCTAGAATATTTGCGTCAGCTCACAAATACGCCTTTTccc                    >  2:2276268/1‑100 (MQ=255)
                                                                             ttCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCt                   <  2:926495/100‑1 (MQ=255)
                                                                             ttCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCt                   <  1:2953188/100‑1 (MQ=255)
                                                                              tcgccATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTg                  <  1:1748031/97‑1 (MQ=255)
                                                                              tCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCGGTAACACCTTATTTACATCCTCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTg                  <  2:1231860/100‑1 (MQ=255)
                                                                               cACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTg                  <  1:2530817/100‑1 (MQ=255)
                                                                               cACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTgg                 <  1:87742/100‑1 (MQ=255)
                                                                                aaaaaaaacggcgCATAAGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGt                >  1:451593/9‑100 (MQ=255)
                                                                                 ccATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTa               <  2:677019/100‑1 (MQ=255)
                                                                                 aaaaaaacggcgCATAAGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTa               >  1:686984/8‑100 (MQ=255)
                                                                                   aaaaacggcgCATAAGCGCCGTTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaa              >  1:667297/6‑100 (MQ=255)
                                                                                    tGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaaa            >  1:2609623/1‑100 (MQ=255)
                                                                                       gcggcgCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt         <  1:872049/100‑1 (MQ=255)
                                                                                          gcgcTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt         <  2:3409164/97‑1 (MQ=255)
                                                                                           cgcTTATGCGGCG‑TTTTTTTTTCGTTCGCACCGTATTTACATCCCCATAGATTATTTGCGGCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGa       <  2:3554918/98‑1 (MQ=255)
                                                                                           cgcTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAtt     >  2:3746789/1‑100 (MQ=255)
                                                                                            gcTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAttt    >  1:773425/1‑100 (MQ=255)
                                                                                            gcTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAttt    <  1:245665/100‑1 (MQ=255)
                                                                                            gcTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTACCAGGGaaaaagcgtatt    <  2:2268410/100‑8 (MQ=255)
                                                                                            gcTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACAGCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAttt    <  1:3290281/100‑1 (MQ=255)
                                                                                              ttATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTcc  <  2:198239/100‑1 (MQ=255)
                                                                                              |                                                                                                    
GAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCC  >  NC_000913/1369932‑1370125

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