Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,653,372 T→C 15.2% intergenic (+33/+785) sseA → / ← sseB 3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,653,3720TC15.2% 286.8 / 42.4 125intergenic (+33/+785)sseA/sseB3‑mercaptopyruvate sulfurtransferase/rhodanase‑like enzyme, sulfur transfer from thiosulfate
Reads supporting (aligned to +/- strand):  ref base T (45/61);  new base C (9/10);  total (54/71)
Fisher's exact test for biased strand distribution p-value = 8.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.88e-01

Only 100 of 126 total aligned reads displayed.

CTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCAC  >  NC_000913/2653274‑2653469
                                                                                                  |                                                                                                 
 tGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTgg                                                                                                 <  1:2454805/100‑1 (MQ=255)
 tGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTg                                                                                                  <  1:1441806/99‑1 (MQ=255)
  ccgatctggCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGc                                                                                                <  1:1136501/93‑1 (MQ=255)
   tACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCg                                                                                               <  2:2505541/100‑1 (MQ=255)
   tACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCg                                                                                               <  2:1597077/100‑1 (MQ=255)
    acgacgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCgg                                                                                              <  1:1997969/100‑1 (MQ=255)
     cgacgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGt                                                                                             >  2:1060160/1‑100 (MQ=255)
     caacgGCGCATGGGGTGAATGGGGCGCACGAGCAGATTTACCAGTTGAGCCCGTAAAATAAGCTTTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGt                                                                                             >  2:1748618/3‑100 (MQ=255)
      gacgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGtt                                                                                            <  2:1337373/100‑1 (MQ=255)
        cgGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGtttt                                                                                          <  2:1011893/100‑1 (MQ=255)
         ggCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGttttt                                                                                         >  1:58773/1‑100 (MQ=255)
          gcgcATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGtttttt                                                                                        <  2:2527173/100‑1 (MQ=255)
          gcgcATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGtttttt                                                                                        >  2:2385620/1‑100 (MQ=255)
           cgcATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTa                                                                                       <  1:1232110/100‑1 (MQ=255)
           cgcATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTa                                                                                       <  2:588992/100‑1 (MQ=255)
            gcATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTAt                                                                                      <  2:2248422/100‑1 (MQ=255)
             cATGGAGTGAATGGGGCGCGCGGGCAGATTTGCCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTAtt                                                                                     <  1:1745038/100‑1 (MQ=255)
              aTGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTg                                                                                    <  2:508955/100‑1 (MQ=255)
              aTGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTg                                                                                    <  2:1439320/100‑1 (MQ=255)
              aTGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTATTg                                                                                    <  2:330531/100‑1 (MQ=255)
              aTGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGTTTTTTATTg                                                                                    <  2:1319824/100‑1 (MQ=255)
               tGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGc                                                                                   <  2:714330/100‑1 (MQ=255)
               tGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGc                                                                                   <  2:1233072/100‑1 (MQ=255)
               tGGAGTGAATGGGGCGCGCGGGAAGATTTACCGGTTGAGCGAGTTAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGc                                                                                   <  2:537191/100‑1 (MQ=255)
                ggAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCt                                                                                  <  1:2050821/100‑1 (MQ=255)
                 gAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTa                                                                                 >  2:1133774/1‑100 (MQ=255)
                  aGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAg                                                                                <  1:1839989/100‑1 (MQ=255)
                          gggCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTc                                                                        >  2:2017640/1‑100 (MQ=255)
                             cgcgcgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTcgcg                                                                     <  2:1979893/100‑1 (MQ=255)
                             cgcgcgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTcgcg                                                                     <  1:1935975/100‑1 (MQ=255)
                                gcgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGcc                                                                  <  1:294824/100‑1 (MQ=255)
                                gcgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGcc                                                                  <  1:2477505/100‑1 (MQ=255)
                                 cgGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCt                                                                 >  1:1177765/1‑100 (MQ=255)
                                   ggCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCttt                                                               >  1:1349836/1‑100 (MQ=255)
                                     cAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTcc                                                             >  1:2018446/1‑100 (MQ=255)
                                      aGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCca                                                            <  2:1472176/100‑1 (MQ=255)
                                       gATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcag                                                           >  2:767370/1‑100 (MQ=255)
                                       gATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcag                                                           >  1:579085/1‑100 (MQ=255)
                                       gATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcag                                                           >  1:1291953/1‑100 (MQ=255)
                                        aTTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcag                                                           <  1:530066/99‑1 (MQ=255)
                                         tttACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcagca                                                         >  2:246062/1‑100 (MQ=255)
                                            aCCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGt                                                      >  2:1268023/1‑100 (MQ=255)
                                             ccGGTTGAGCCCGTAAAATAAGCTTTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTGGTCTGGTTCGCGGCCTTTCCAGCAGGtt                                                     <  2:2001193/100‑1 (MQ=25)
                                               ggTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGa                                                   >  2:114196/1‑100 (MQ=255)
                                                 ttGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACt                                                 <  2:2085336/100‑1 (MQ=255)
                                                   gAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTtg                                               <  1:1160503/100‑1 (MQ=255)
                                                     gCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTtgtg                                             >  2:2397177/1‑100 (MQ=255)
                                                      ccAGGGAAATAGGTATGTGCCAGGCCATAAAGAACTGCCGAATTTGGCGGTGTTTTATTGCGAGTTTGGTTCGCGGCCTTTCCAGCAGGTTGACGtgtg                                             <  1:2192038/99‑1 (MQ=21)
                                                       cAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGtt                                           <  2:1397131/100‑1 (MQ=255)
                                                       cAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGtt                                           <  2:1929003/100‑1 (MQ=255)
                                                       cAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGtt                                           >  1:1237193/1‑100 (MQ=255)
                                                           gAAATAATTATTTTACAGGCAATAATCTATCGCCGGATTTGGCGGTTTTTTATTGCTAGTCTTGTTCGCGGCCTTTCCAGCAGGTTGACGTGTGTTa                                          <  1:1355319/97‑1 (MQ=25)
                                                           gAAATAAGTATTTTACAGGCAATAAAAAACCGCCGGATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACa                                        <  1:1143607/99‑1 (MQ=255)
                                                           gAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACAt                                       >  1:2332811/1‑100 (MQ=255)
                                                           gAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACAt                                       >  2:2570163/1‑100 (MQ=255)
                                                            aaaTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATg                                      <  2:579087/100‑1 (MQ=255)
                                                            aaaTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATg                                      <  1:2522234/100‑1 (MQ=255)
                                                            aaaTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATg                                      <  2:1177769/100‑1 (MQ=255)
                                                            aaaTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATg                                      <  2:124729/100‑1 (MQ=255)
                                                              aTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAg                                    >  2:897440/1‑100 (MQ=255)
                                                              aTAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAg                                    <  2:507580/100‑1 (MQ=255)
                                                               tAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGc                                   <  2:2018458/100‑1 (MQ=255)
                                                               tAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGc                                   >  1:466769/1‑100 (MQ=255)
                                                                aaGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTATCGTTCGCGGCCTTTCGAGCATGTTTCCTTGTGTTACATGAGc                                   >  1:1017441/1‑99 (MQ=37)
                                                                    aTTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGc                              >  1:684315/1‑100 (MQ=255)
                                                                     ccagactagCAATAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCa                             >  1:2181734/9‑100 (MQ=25)
                                                                     ccagactagCAATAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCa                             >  1:224191/9‑100 (MQ=25)
                                                                       ttACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAgg                           <  2:1163896/100‑1 (MQ=255)
                                                                       ttACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAgg                           <  1:1961559/100‑1 (MQ=255)
                                                                        tACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTGCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAgg                           <  1:430949/99‑1 (MQ=255)
                                                                         aCAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTg                         <  2:1877286/100‑1 (MQ=255)
                                                                          ctagCAATAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGc                        >  2:1117394/4‑100 (MQ=37)
                                                                          cAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGc                        <  1:1177375/100‑1 (MQ=255)
                                                                          cAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTAACATGAGCAACGCGGGTGc                        >  2:2199993/1‑100 (MQ=255)
                                                                           ccgatctAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCt                       <  1:54043/94‑1 (MQ=255)
                                                                           aGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCt                       <  1:2490977/100‑1 (MQ=255)
                                                                           aGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCt                       <  2:192805/100‑1 (MQ=255)
                                                                           aGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCt                       <  2:780650/100‑1 (MQ=255)
                                                                            ggCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCtt                      <  1:720212/100‑1 (MQ=255)
                                                                            ggCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCtt                      >  2:663912/1‑100 (MQ=255)
                                                                             gCAATAAAAAACCGCCGCATTTGGCGGTTTTTTATGGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTc                     <  1:321186/100‑1 (MQ=255)
                                                                             gCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACGTGTGTTACATGAGCAACGCAGGTGCTTc                     >  1:169468/1‑100 (MQ=255)
                                                                               aaTAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCAGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTcc                    >  1:2040430/1‑98 (MQ=37)
                                                                                 tAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAg                 <  1:675942/100‑1 (MQ=255)
                                                                                 tAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAg                 >  2:2588442/1‑100 (MQ=255)
                                                                                 tAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAg                 <  1:1916948/100‑1 (MQ=255)
                                                                                  aaaaaaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGc                >  2:891843/1‑100 (MQ=255)
                                                                                  aaaaaaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCCGGATCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCtagatcgg                >  2:54044/1‑93 (MQ=37)
                                                                                  aaaaaaCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAAGATGACATGTGCGACATGAGCGGCGCAGGTGCGGCACAg                 >  1:1838200/1‑99 (MQ=18)
                                                                                   aaaaaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGTCTTATGTTAGTTGAGCAACGCAGGTGCTTCGCAg                 >  1:2288290/1‑98 (MQ=37)
                                                                                   aaaaaCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCa               >  2:237077/1‑100 (MQ=255)
                                                                                      aaCCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCaaaa            <  1:2470386/100‑1 (MQ=255)
                                                                                        ccgccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAAc           <  1:1906007/99‑1 (MQ=255)
                                                                                        ccgccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACa          >  1:232206/1‑100 (MQ=255)
                                                                                         cgccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACaa         <  2:2445178/100‑1 (MQ=255)
                                                                                          gccgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACCTGAGCAACGCAGGTGCTTCACAGCAAAACAAt        >  2:900047/1‑100 (MQ=255)
                                                                                          gccgAAGTTGGCGGATTTTGATGGCTAGTAAGGTCCGCGGCCTTACCAGCAGGTTGGCTTGTGTTACATGAGCAACGAAGGTACTTCACAGCAAAACa          <  1:1781598/98‑1 (MQ=255)
                                                                                            cgAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTGCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGAAAAACAATAc      <  1:121564/100‑1 (MQ=255)
                                                                                               aTTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCa   <  1:1675099/100‑1 (MQ=255)
                                                                                                tttGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCAc  <  1:2385605/100‑1 (MQ=255)
                                                                                                  |                                                                                                 
CTGTACGACGGCGCATGGAGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCAC  >  NC_000913/2653274‑2653469

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: