Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 175,150:1 +G 15.8% coding (44/1422 nt) clcA → H(+)/Cl(‑) exchange transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913175,1491.G15.8% 374.2 / 76.8 120coding (43/1422 nt)clcAH(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base . (48/53);  new base G (11/8);  total (59/61)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.32e-01

Only 100 of 123 total aligned reads displayed.

ACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGC  >  NC_000913/175053‑175248
                                                                                                 |                                                                                                   
aCGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGcc                                                                                                  <  2:1431878/100‑1 (MQ=255)
 cGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCt                                                                                                 <  2:303879/100‑1 (MQ=255)
  gTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTg                                                                                                <  1:2045021/100‑1 (MQ=255)
  gTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTg                                                                                                >  2:2139155/1‑100 (MQ=255)
   tGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGc                                                                                               <  1:2555246/100‑1 (MQ=255)
     cGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCGCGCGGCCTGCgg                                                                                             <  1:1850121/100‑2 (MQ=255)
      gAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCGCGCGGCCTGCggt                                                                                            <  1:2385501/100‑3 (MQ=255)
      gAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCGCGCGGCCTGCggt                                                                                            <  1:80100/100‑3 (MQ=255)
      gAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCGCGCGGCCTGCggt                                                                                            <  2:1771309/100‑3 (MQ=255)
       atgaaaaTGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCGCGCGGCCTGCggtg                                                                                           <  1:2533183/97‑4 (MQ=255)
       aGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACg                                                                                           <  1:378522/100‑1 (MQ=255)
       aGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCGCGCGGCCTGCggtg                                                                                           <  2:1693259/100‑4 (MQ=255)
        gTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCGCGCGGCCTGCggtgt                                                                                          <  1:2139028/100‑5 (MQ=255)
         tAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCa                                                                                         <  2:321013/100‑1 (MQ=255)
           aaaTGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGa                                                                                       >  1:805405/1‑100 (MQ=255)
            aaTGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCGCGCGGCCTGCggtgtttct                                                                                      <  2:1437282/100‑9 (MQ=38)
               gCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAAc                                                                                   <  2:709637/100‑1 (MQ=255)
                ccGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACt                                                                                  >  2:2287270/1‑100 (MQ=255)
                ccGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACt                                                                                  <  1:2583093/100‑1 (MQ=255)
                 cGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAAc                                                                                   >  1:2215408/1‑98 (MQ=255)
                 cGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTg                                                                                 <  1:785756/100‑1 (MQ=255)
                      tCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCg                                                                            <  2:878016/100‑1 (MQ=255)
                      tCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCg                                                                            >  2:115832/1‑100 (MQ=255)
                      tCTTATTTGTATGACCAATAAGTGATCATTGGATCAAAACTGATACTACCTCTTTAGAAAAACAGCAGGCCGCGCAACTGCGACGCAGACAAAAGATTc                                                                             >  1:1155693/1‑99 (MQ=255)
                       cTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGc                                                                           >  2:178873/1‑100 (MQ=255)
                       cTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGc                                                                           <  1:74356/100‑1 (MQ=255)
                       cTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGc                                                                           >  1:1460390/1‑100 (MQ=255)
                        ttATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGcc                                                                          <  1:960568/100‑1 (MQ=255)
                         tATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCa                                                                         >  1:953165/1‑100 (MQ=255)
                         tATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGACGCGCGCCTGCGACGCAGACAACTGATTCGCCa                                                                         >  1:598214/1‑100 (MQ=255)
                             tGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAActt                                                                     <  2:2434592/100‑1 (MQ=255)
                               tATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAActtct                                                                   <  1:1978654/100‑1 (MQ=255)
                               tATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAActtct                                                                   <  1:2545013/100‑1 (MQ=255)
                                aTGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATCCGCCAActtc                                                                    >  1:849806/1‑98 (MQ=255)
                                 tGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTg                                                                 >  1:2022361/1‑100 (MQ=255)
                                 tGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTg                                                                 >  1:1001043/1‑100 (MQ=255)
                                   aCCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGACAACTTCTTGAg                                                               <  2:1703552/100‑1 (MQ=255)
                                    ccAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAgc                                                              >  2:2223082/1‑100 (MQ=255)
                                    ccAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAgc                                                              <  2:2308784/100‑1 (MQ=255)
                                    ccAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAg                                                               <  1:1032425/99‑1 (MQ=255)
                                     cAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAgcg                                                             <  2:214973/100‑1 (MQ=255)
                                      aaTAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAgcgc                                                            >  1:1523812/1‑100 (MQ=255)
                                      aaTAAGTGATCATTGGATGAAAACTGATACTCCCTCTTAAAAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTgcgcgc                                                            >  1:1594233/1‑100 (MQ=255)
                                       aTAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGcgcg                                                           <  2:1549260/100‑1 (MQ=255)
                                       aTAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGcgcg                                                           >  1:1549251/1‑100 (MQ=255)
                                       aTAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGcgcg                                                           <  1:2369672/100‑1 (MQ=255)
                                         aaGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGAt                                                         <  2:638890/100‑1 (MQ=255)
                                         aaGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCNTGAGCGCGAt                                                         <  2:874762/100‑1 (MQ=255)
                                          aGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATa                                                        >  1:2009027/1‑100 (MQ=255)
                                              aTCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAAc                                                    <  2:677029/100‑1 (MQ=255)
                                              aTCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAAc                                                    <  2:2011789/100‑1 (MQ=255)
                                                 aTTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTTGCCAACTTCTCGAGAGCGATAAAAcccc                                                 >  2:2202931/1‑100 (MQ=255)
                                                   tGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGt                                               <  1:1530351/100‑1 (MQ=255)
                                                    ggATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGtt                                              >  1:87679/1‑100 (MQ=255)
                                                    ggATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGtt                                              <  1:2214579/100‑1 (MQ=255)
                                                     gATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTa                                             >  2:1316845/1‑100 (MQ=255)
                                                      aTGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAg                                            <  1:1286315/100‑1 (MQ=255)
                                                      aTGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAg                                            >  1:1147797/1‑100 (MQ=255)
                                                       tGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGc                                           <  1:184207/100‑1 (MQ=255)
                                                            aCTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCAttt                                      >  2:2224167/1‑100 (MQ=255)
                                                            aCTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCAttt                                      <  2:1135339/100‑1 (MQ=255)
                                                              tGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTg                                    <  2:1869101/100‑1 (MQ=255)
                                                              tGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTg                                    <  2:1966840/100‑1 (MQ=255)
                                                                aTACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGtt                                  <  2:810614/100‑1 (MQ=255)
                                                                 tACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGttt                                 <  2:2158679/100‑1 (MQ=255)
                                                                 tACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGttt                                 >  2:364486/1‑100 (MQ=255)
                                                                   cTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTAt                               >  2:1569848/1‑100 (MQ=255)
                                                                    tCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATg                              <  2:726048/100‑1 (MQ=255)
                                                                     cccTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATgg                             <  1:1178763/100‑1 (MQ=255)
                                                                      ccTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATggc                            >  2:525735/1‑100 (MQ=255)
                                                                       ctctTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATggcg                           >  2:1488927/1‑100 (MQ=255)
                                                                       ctctTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATggcg                           <  2:1386334/100‑1 (MQ=255)
                                                                       ctctTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATggcg                           >  2:1290241/1‑100 (MQ=255)
                                                                          tttAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCa                        >  1:2549219/1‑100 (MQ=255)
                                                                           ttAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAg                       <  1:2311283/100‑1 (MQ=255)
                                                                            tgtctgcgtcGCAGGCGCGCGGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAgt                      >  2:2318438/10‑100 (MQ=38)
                                                                              gAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAgtcg                    >  2:552989/1‑100 (MQ=255)
                                                                              gAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAgtcg                    >  2:386276/1‑100 (MQ=255)
                                                                              gAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAgtcg                    >  2:1414807/1‑100 (MQ=255)
                                                                              gAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAgtcg                    <  2:1143713/100‑1 (MQ=255)
                                                                              gAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAgtcg                    >  1:479392/1‑100 (MQ=255)
                                                                               ctgcgtcGCAGGCGCGCGGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAgtcgt                   >  2:920330/7‑100 (MQ=255)
                                                                                tgcgtcGCAGGCGCGCGGCCTGCGACGCAGACAACTGCTTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGGAgtcgtc                  >  1:1388968/6‑100 (MQ=255)
                                                                                tgcgtcGCAGGCGCGCGGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAgtcgtc                  >  1:1343506/6‑100 (MQ=255)
                                                                                aaCACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAgtcgtc                  >  2:1991259/1‑100 (MQ=255)
                                                                                 gcgtcGCAGGCGCGCGGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGtcgtcg                 >  2:324392/5‑100 (MQ=255)
                                                                                 acacCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGtcgtcg                 >  2:2319122/1‑100 (MQ=255)
                                                                                  cgtcGCAGGCGCGCGGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCgg                >  1:1721053/4‑100 (MQ=255)
                                                                                   gtcGCAGGCGCGCGGCCTGCGACGCAGACAACTGCTTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGc               >  1:1778305/3‑100 (MQ=255)
                                                                                   gtcGCAGGCGCGCGGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGc               >  2:2202157/3‑100 (MQ=255)
                                                                                   gtcGCAGGCGCGCGGCCTGAGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGc               >  1:1566703/3‑100 (MQ=255)
                                                                                   acCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGc               <  2:2353946/100‑1 (MQ=255)
                                                                                    tcGCAGGCGCGCGGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCa              >  2:1895951/2‑100 (MQ=255)
                                                                                    ccGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGc               <  1:407530/99‑1 (MQ=255)
                                                                                    ccGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCa              >  2:147074/1‑100 (MQ=255)
                                                                                    ccGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCa              <  1:2223067/100‑1 (MQ=255)
                                                                                     cGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCAc             <  1:223202/100‑1 (MQ=255)
                                                                                            cgcgcgCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCttgtt      >  2:1661665/1‑100 (MQ=255)
                                                                                             gcgcgcCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCttgttg     <  1:2117432/100‑1 (MQ=255)
                                                                                                cgcCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGc  >  2:574355/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
ACGTGCGAGTAAAATGCCGTCATCTTATTTGTATGACCAATAAGTGATCATTGGATGAAAACTGATACTCCCTCTTTAGAAACACCGCAGGCCGCGCGCCTGCGACGCAGACAACTGATTCGCCAACTTCTTGAGCGCGATAAAACCCCGTTAGCCATTTTGTTTATGGCGGCAGTCGTCGGCACGCTTGTTGGGC  >  NC_000913/175053‑175248

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