Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,891,179 G→C 19.2% S240R (AGC→AGG cysJ ← sulfite reductase, alpha subunit, flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,891,1790GC19.2% 332.1 / 69.8 149S240R (AGC→AGGcysJsulfite reductase, alpha subunit, flavoprotein
Reads supporting (aligned to +/- strand):  ref base G (49/69);  new base C (14/14);  total (66/83)
Fisher's exact test for biased strand distribution p-value = 5.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

Only 100 of 155 total aligned reads displayed.

ACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGC  >  NC_000913/2891080‑2891278
                                                                                                   |                                                                                                   
aCCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAgg                                                                                                     <  1:2857130/100‑1 (MQ=255)
aCCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAgg                                                                                                     <  1:3749642/100‑1 (MQ=255)
 cccGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGc                                                                                                    <  1:1114366/100‑1 (MQ=255)
   cGGCTGGTGACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGTGTTATGCCCGGTAATTTTCTGGTTAACAGAGAGGCTa                                                                                                  <  2:752640/100‑1 (MQ=255)
    ggCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGTGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAg                                                                                                 >  2:3765831/1‑100 (MQ=255)
     gCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGc                                                                                                >  1:755460/1‑100 (MQ=255)
     gCTGGTAACGCATGCCCGAGTCACCTAAGTCAATGTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGc                                                                                                >  2:2461304/1‑100 (MQ=255)
      cTGGTAACGCATGCCCGAGTCACCTAATTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCTGTAATTTTCTGGTTAACAGAGAGGCTAGcc                                                                                               <  2:2068914/100‑1 (MQ=255)
        ggTAACGCATGTCCTATTCACCTAAGTCAATTTCGATATGGTGAACGTCTTTGTCATAGTTACGCTCGGTAATTTTCTGGTTAACAGAGAGGCTAGcc                                                                                               <  2:1345853/98‑1 (MQ=255)
          tAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATTGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGccacca                                                                                           <  2:470763/100‑1 (MQ=255)
           aaCGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAg                                                                                          <  1:453322/100‑1 (MQ=255)
            aCGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAgc                                                                                         <  2:2464526/100‑1 (MQ=255)
                aTGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAgcggcg                                                                                     <  1:1764169/100‑1 (MQ=255)
                  gCCCGCGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTTATTAACAGAGAGGCTATCCACCAGCGCcgcg                                                                                   >  2:896106/1‑100 (MQ=255)
                    ccGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTc                                                                                 >  1:519522/1‑100 (MQ=255)
                     cGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCt                                                                                >  1:3846897/1‑100 (MQ=255)
                      gAGTCACCTAAGTCAATTTCGATATGGTGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCt                                                                                <  1:1321578/99‑1 (MQ=255)
                       aGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCttt                                                                              <  2:3213901/100‑1 (MQ=255)
                        gTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTg                                                                             <  1:2217700/100‑1 (MQ=255)
                        gTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTg                                                                             <  2:2146488/100‑1 (MQ=255)
                         tCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTg                                                                             <  1:764648/99‑1 (MQ=255)
                         tCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGc                                                                            <  1:1115312/100‑1 (MQ=255)
                         tCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGc                                                                            <  1:3889026/100‑1 (MQ=255)
                         tCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGc                                                                            <  1:2936039/100‑1 (MQ=255)
                         tCACCTAAGTAAATTTCGATATGGCGAACGTCTTTTTCAGATTTACGTCCGGTAATTTTCTGGTTAACAGAGACGCTAGTCACCAGAGGcgcgcaggtga                                                                            <  2:452071/100‑8 (MQ=255)
                           aCCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCt                                                                           <  1:2771631/99‑1 (MQ=255)
                            ccTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTg                                                                          <  2:3836725/99‑1 (MQ=255)
                            ccTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGt                                                                         >  2:2005271/1‑100 (MQ=255)
                            ccTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGt                                                                         >  1:3179518/1‑100 (MQ=255)
                             cTAAGTCAATTTCGATATGGAGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTa                                                                        >  2:1883851/1‑100 (MQ=255)
                              tAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTAc                                                                       >  1:3891136/1‑100 (MQ=255)
                               aaGTCAATTTCGATATGGCGAACTTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACg                                                                      <  1:197163/100‑1 (MQ=255)
                               aaGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACg                                                                      >  2:1842012/1‑100 (MQ=255)
                               aaGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACg                                                                      <  1:1243245/100‑1 (MQ=255)
                                 gTCAATTTCGATATGGCGAATGTCTTTTTCAGAGATACGCCCGGTAAATTTCTGGTTAACAGAGAGACTAGCCCGCAGCGGCGCGGCtttaat                                                                           >  2:3447517/1‑90 (MQ=255)
                                 gTCAATTTCGACATGACGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACgg                                                                     <  2:7474/99‑1 (MQ=255)
                                  tCAATTTCGATATGGCGAACTTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACg                                                                      <  2:1058234/97‑1 (MQ=255)
                                  tCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGc                                                                   >  2:2842690/1‑100 (MQ=255)
                                   cAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACggg                                                                    >  2:173786/1‑98 (MQ=255)
                                     aTTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTgg                                                                >  1:3244159/1‑100 (MQ=255)
                                     aTTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTgg                                                                >  2:1293239/1‑100 (MQ=255)
                                      tttCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGgt                                                               >  1:1341211/1‑100 (MQ=255)
                                        tCGATATGGTGAACGTCTTTTTCAGAGTTACGCTTGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGgtgt                                                             <  1:2330160/100‑1 (MQ=255)
                                         cGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGtgtg                                                            <  2:3894366/100‑1 (MQ=255)
                                         cGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGtgtg                                                            <  2:863957/100‑1 (MQ=255)
                                         cGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGtgtg                                                            >  1:2484414/1‑100 (MQ=255)
                                             atGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGAtt                                                        <  1:1360732/100‑1 (MQ=255)
                                             atGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGACTAGCCACCAGCGGAGCGTCTTTGCTGTAAGGGCTGGTGTGGAt                                                         >  2:169235/1‑99 (MQ=255)
                                              tGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGAttt                                                       >  2:2034384/1‑100 (MQ=255)
                                               ggCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTc                                                      >  2:1242492/1‑100 (MQ=255)
                                               ggCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTc                                                      >  2:1573621/1‑100 (MQ=255)
                                               ggCGAACGTCTTTTTCAGAGTTAAGCCCGGTAATTTTATGGTTAACAGAGAGGCTAGCCACCAGAGGCGAGTCTTTGCTGTAAGGGCTGGTGTGGATTTc                                                      >  2:2325708/1‑100 (MQ=255)
                                                gCGAACGTCTTTTTCAGAGTTACTCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCa                                                     <  1:920791/100‑1 (MQ=255)
                                                gCGAACGTCTTTTTCAGAGTTACGCCCGGTATTTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCa                                                     <  1:1514201/100‑1 (MQ=255)
                                                gCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCa                                                     >  1:684790/1‑100 (MQ=255)
                                                 cGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCAt                                                    <  1:1242506/100‑1 (MQ=255)
                                                 cGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCAt                                                    <  2:2891720/100‑1 (MQ=255)
                                                  gAACGTCTTTTTCAGAGTTAAGCAAGGTAATTTTCTGGTGAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTCGATTTCAtt                                                   >  2:634339/1‑100 (MQ=255)
                                                   aaCGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCAttt                                                  <  2:1197473/100‑1 (MQ=255)
                                                       tCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACcg                                              >  1:2021885/1‑100 (MQ=255)
                                                        cTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACcgc                                             >  1:2753826/1‑100 (MQ=255)
                                                        cTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACcgc                                             <  2:2669567/100‑1 (MQ=255)
                                                        cTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGACTTTGGTGTACGGGCTGGTGTGGATTTCATTTACcgc                                             >  2:2633843/1‑100 (MQ=255)
                                                           tttCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGcc                                          >  1:3003719/1‑100 (MQ=255)
                                                            ttCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCa                                         >  2:1070320/1‑100 (MQ=255)
                                                               agagTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTa                                      <  2:684782/100‑1 (MQ=255)
                                                                 agTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGc                                    <  1:1200029/100‑1 (MQ=255)
                                                                  gTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGc                                    <  2:1466136/99‑1 (MQ=255)
                                                                  gTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCg                                   <  2:1172155/100‑1 (MQ=255)
                                                                      cGCCCGGTAATTTTCTGCTTAACAGAGAGCCTAGCCACCAGCTGCGCGTCTTTGCTGTACGGGCTGGTGTGTATTTATTTTAACGCCCTAGCAGCTAc                                 >  1:322039/1‑98 (MQ=255)
                                                                        cccGGTAATTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGAt                             >  2:1405478/1‑100 (MQ=255)
                                                                          cGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTg                           <  1:3455904/100‑1 (MQ=255)
                                                                          cGGTAATTTTCTGGTTAACAGAAAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGCAGCCACGGATTg                           >  2:2265732/1‑100 (MQ=255)
                                                                           ggTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGc                          <  1:2663433/100‑1 (MQ=255)
                                                                           ggTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGc                          <  1:1193406/100‑1 (MQ=255)
                                                                             tAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGa                        <  1:1842027/100‑1 (MQ=255)
                                                                              aaTTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGaa                       <  2:1698638/100‑1 (MQ=255)
                                                                              aaTTTTCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGATGAAAAAGATTTAATATACCGCGCAAGTAGCGACGGATTGCGaa                       >  1:2454781/1‑100 (MQ=255)
                                                                               aTTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGAATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAg                      >  1:3326284/1‑100 (MQ=255)
                                                                                ttttCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAgn                     <  2:1154769/100‑2 (MQ=255)
                                                                                  ttCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTg                   <  2:3711578/100‑1 (MQ=255)
                                                                                  ttCTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTg                   <  2:3269472/100‑1 (MQ=255)
                                                                                    cTGGTTAACAGAGAGCCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGcc                 <  1:3666404/100‑1 (MQ=255)
                                                                                      ggTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCcgc               <  1:2395230/100‑1 (MQ=255)
                                                                                      ggTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCcgc               <  2:2499174/100‑1 (MQ=255)
                                                                                      ggTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCcgc               >  1:809471/1‑100 (MQ=255)
                                                                                       gTTAACAGAGAGTCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCcgcg              >  2:2975846/1‑100 (MQ=255)
                                                                                       gTTAACAGAGAGGCTAGTTATTAGCGGCGCGTTTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCcgcg              <  1:581509/100‑1 (MQ=255)
                                                                                       gTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGGTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCcgcg              >  2:2077003/1‑100 (MQ=255)
                                                                                       gTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGTTTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCcgcg              <  2:2996100/100‑1 (MQ=255)
                                                                                       gTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCcgcg              <  1:3248267/100‑1 (MQ=255)
                                                                                         tAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGAc            <  2:2151629/100‑1 (MQ=255)
                                                                                         tAACAGAGAGCCTAGCCACCAGTGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGAc            <  2:939588/100‑1 (MQ=255)
                                                                                          aaCAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGa             <  2:3525557/98‑1 (MQ=255)
                                                                                           aCAGTGCGGCTAGCAACCAGCGGAGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAg          >  2:3729634/1‑100 (MQ=255)
                                                                                               agagGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGgcgc      <  1:3844525/100‑1 (MQ=255)
                                                                                                 agGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCAc    <  2:3589503/100‑1 (MQ=255)
                                                                                                  ggCTAGCCACCAGCGGTGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCAc    >  1:48471/1‑99 (MQ=255)
                                                                                                  ggCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCAcg   >  1:3705919/1‑100 (MQ=255)
                                                                                                  ggCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCAcg   >  1:2839782/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGC  >  NC_000913/2891080‑2891278

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