Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 882,722 G→A R5C (CGC→TGC)  deoR ← deoxyribose‑5‑phosphate‑inducible deoxyribose operon transcriptional repressor; repressor of nupG and tsx

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913882,7220GA96.8% 316.1 / 0.2 94R5C (CGC→TGC) deoRdeoxyribose‑5‑phosphate‑inducible deoxyribose operon transcriptional repressor; repressor of nupG and tsx
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base A (45/46);  total (45/49)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

CATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGAT  >  NC_000913/882627‑882818
                                                                                               |                                                                                                
cATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTACGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACg                                                                                               >  2:462954/1‑99 (MQ=255)
 aTCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACgtg                                                                                             >  1:1109385/1‑100 (MQ=255)
  tctcCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACgtgt                                                                                            >  1:1033695/1‑100 (MQ=255)
   ctcCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGtt                                                                                           >  2:692148/1‑100 (MQ=255)
     ccGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTc                                                                                         <  1:1045189/100‑1 (MQ=255)
      caaaaCCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTcc                                                                                        >  1:1574287/3‑100 (MQ=255)
        aaaCCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCAt                                                                                      <  1:868307/100‑1 (MQ=255)
        aaaCCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCAt                                                                                      <  2:458708/100‑1 (MQ=255)
         aaCCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGTTTCCATa                                                                                     <  1:367799/100‑1 (MQ=255)
         aaCCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATa                                                                                     <  2:941819/100‑1 (MQ=255)
         aaCCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATa                                                                                     <  1:609838/100‑1 (MQ=255)
         aaCCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATa                                                                                     <  2:1288455/100‑1 (MQ=255)
          aCCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAGCGTGTTTCCATaa                                                                                    >  1:1362286/1‑100 (MQ=255)
          aCCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATaa                                                                                    >  1:1753566/1‑100 (MQ=255)
          aCCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATaa                                                                                    >  1:1471333/1‑100 (MQ=255)
           ccccAAGCAGGGCGGCGGCGTCTTTTAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAAt                                                                                   >  2:827826/1‑100 (MQ=255)
           ccccAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAAt                                                                                   >  1:226484/1‑100 (MQ=255)
           ccccAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAAt                                                                                   >  2:1192825/1‑100 (MQ=255)
                aGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGTTTCCATAATCCctc                                                                              <  1:1942912/100‑1 (MQ=255)
                 gCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCctct                                                                             <  2:1574282/100‑1 (MQ=255)
                  ctgggCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTg                                                                            <  2:843685/98‑1 (MQ=255)
                  cAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTg                                                                            <  2:1181001/100‑1 (MQ=255)
                    gggCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGag                                                                          >  1:843687/1‑99 (MQ=255)
                    gggCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGaa                                                                          <  1:1591112/100‑1 (MQ=255)
                     ggcggcggcGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAAt                                                                         <  2:430693/100‑1 (MQ=255)
                     ggcggcggcGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAAt                                                                         <  2:530899/100‑1 (MQ=255)
                       cggcggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAg                                                                       >  2:2004110/1‑100 (MQ=255)
                       cggcggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAg                                                                       >  2:829514/1‑100 (MQ=255)
                        ggcggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGt                                                                      <  2:1392936/100‑1 (MQ=255)
                        ggcggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGt                                                                      <  1:818353/100‑1 (MQ=255)
                           ggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTAt                                                                   <  2:126841/100‑1 (MQ=255)
                            gcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTAtt                                                                  >  1:944374/1‑100 (MQ=255)
                            gcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTAtt                                                                  <  1:1834446/100‑1 (MQ=255)
                            gcgTCTTTAAGATGTAACTTATCACTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTAtt                                                                  <  2:2084363/100‑1 (MQ=255)
                             cgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTg                                                                 <  2:1471329/100‑1 (MQ=255)
                              gTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGa                                                                <  2:1753561/100‑1 (MQ=255)
                               tCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGaa                                                               <  1:592689/100‑1 (MQ=255)
                                cTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAg                                                              >  2:1110259/1‑100 (MQ=255)
                                cTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAg                                                              <  1:913659/100‑1 (MQ=255)
                                cTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAg                                                              >  2:2178954/1‑100 (MQ=255)
                                 tttAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGc                                                             >  2:189829/1‑100 (MQ=255)
                                 tttAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGc                                                             >  1:894443/1‑100 (MQ=255)
                                 tttAAGATGTAACTTATCGCTGCGTTTTAATTCTTGAAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTGCCATAATCCCTCTGAATAGTTATTGAAGc                                                             >  2:774186/1‑100 (MQ=255)
                                   tAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGa                                                           >  2:1518255/1‑100 (MQ=255)
                                   tAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGa                                                           >  1:1273868/1‑100 (MQ=255)
                                      gATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGcc                                                        >  1:2041492/1‑100 (MQ=255)
                                       aTGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGcc                                                        <  1:2070371/99‑1 (MQ=255)
                                          tAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTc                                                    <  2:467373/100‑1 (MQ=255)
                                           aaCTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCa                                                   >  2:1680067/1‑100 (MQ=255)
                                            aCTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCaa                                                  >  2:2053184/1‑100 (MQ=255)
                                            aCTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCaa                                                  >  2:1255042/1‑100 (MQ=255)
                                              ttATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAAta                                                >  2:439831/1‑100 (MQ=255)
                                              ttATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAAta                                                <  1:1110261/100‑1 (MQ=255)
                                               tATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAAtac                                               >  1:2461/1‑100 (MQ=255)
                                                  cGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAAtactac                                            <  1:2160910/100‑1 (MQ=255)
                                                    cTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTacac                                          >  2:1876763/1‑100 (MQ=255)
                                                    cTGCGTTTTAATTCTTGCAGCAGCTGCCCGAGACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTaca                                           <  2:1849556/99‑1 (MQ=255)
                                                      gCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACtt                                        >  2:2287046/1‑100 (MQ=255)
                                                         ttttAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACttttt                                     <  2:1646072/100‑1 (MQ=255)
                                                          tttAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTa                                    <  2:944372/100‑1 (MQ=255)
                                                          tttAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTa                                    <  1:2147156/100‑1 (MQ=255)
                                                            tAATTCTTGCAGTAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGc                                  >  1:2025764/1‑100 (MQ=255)
                                                            tAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGc                                  >  2:157680/1‑100 (MQ=255)
                                                            tAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGc                                  <  2:1721797/100‑1 (MQ=255)
                                                            tAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGc                                  >  2:2113165/1‑100 (MQ=255)
                                                             aaTTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCa                                 <  1:2028151/100‑1 (MQ=255)
                                                              aTTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCag                                <  1:2053191/100‑1 (MQ=255)
                                                              aTTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCag                                >  1:1020521/1‑100 (MQ=255)
                                                                 cTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAgaga                             <  1:1980942/100‑1 (MQ=255)
                                                                   tGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGTCGCTCAATACTACACTTTTTAGCAGAGATc                           >  1:1914113/1‑100 (MQ=255)
                                                                   tGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATc                           >  1:2250374/1‑100 (MQ=255)
                                                                   tGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATc                           <  1:1066708/100‑1 (MQ=255)
                                                                   tGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATc                           <  1:692149/100‑1 (MQ=255)
                                                                   tGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATc                           <  2:2250366/100‑1 (MQ=255)
                                                                   tGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATc                           <  2:1362283/100‑1 (MQ=255)
                                                                   tGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATc                           <  1:189831/100‑1 (MQ=255)
                                                                     cagcagCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAg                         >  2:280491/1‑100 (MQ=255)
                                                                      agcagcTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGt                        <  2:2461/100‑1 (MQ=255)
                                                                      agcagcTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGt                        >  1:1673136/1‑100 (MQ=255)
                                                                       gcagcTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTc                       <  2:783497/100‑1 (MQ=255)
                                                                        cagcTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCa                      <  1:462956/100‑1 (MQ=255)
                                                                          gcTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACg                    >  1:424890/1‑100 (MQ=255)
                                                                          gcTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACg                    >  1:167541/1‑100 (MQ=255)
                                                                           cTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGc                   >  1:2113752/1‑100 (MQ=255)
                                                                             gCCCGATACGCTCTTCGCGACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCAc                 <  1:888804/100‑1 (MQ=255)
                                                                             gCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCAc                 <  2:925762/100‑1 (MQ=255)
                                                                               ccGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCAccc               <  2:1020519/100‑1 (MQ=255)
                                                                                  aTACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGc            <  1:114148/100‑1 (MQ=255)
                                                                                   tACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGcc           >  1:1549805/1‑100 (MQ=255)
                                                                                    aCGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCt          >  1:213994/1‑100 (MQ=255)
                                                                                      gCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCttt        >  1:1626549/1‑100 (MQ=255)
                                                                                        tctTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGc      <  2:1776984/100‑1 (MQ=255)
                                                                                         ctTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCg     >  1:963666/1‑100 (MQ=255)
                                                                                          ttCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCgg    <  2:2203310/100‑1 (MQ=255)
                                                                                            cgcaACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGAt  >  1:2339555/1‑100 (MQ=255)
                                                                                               |                                                                                                
CATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGAT  >  NC_000913/882627‑882818

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