Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 882,722 G→A 100% 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,7220GA99.3% 471.3 / ‑4.6 141R5C (CGC→TGC) deoRdeoxyribose‑5‑phosphate‑inducible deoxyribose operon transcriptional repressor; repressor of nupG and tsx
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (61/79);  total (61/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

Only 100 of 142 total aligned reads displayed.

GTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGAT  >  NC_000913/882625‑882818
                                                                                                 |                                                                                                
  cATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACgt                                                                                              >  1:1076476/1‑100 (MQ=255)
   aTCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGCTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACgt                                                                                              >  2:2021586/1‑99 (MQ=255)
    tctcCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACgtgt                                                                                            >  1:1501857/1‑100 (MQ=255)
       ccGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTc                                                                                         <  1:2713937/100‑1 (MQ=255)
       ccGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTc                                                                                         <  1:1925487/100‑1 (MQ=255)
         gAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCa                                                                                       <  1:3183995/100‑1 (MQ=255)
         gAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCa                                                                                       >  2:3600903/1‑100 (MQ=255)
            aCCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATaa                                                                                    <  1:2326501/100‑1 (MQ=255)
             ccccAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATagt                                                                                   >  1:3670009/1‑98 (MQ=255)
             ccccAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAAt                                                                                   >  1:4039749/1‑100 (MQ=255)
             ccccAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAAt                                                                                   >  1:1494695/1‑100 (MQ=255)
             ccccAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAAt                                                                                   <  2:137843/100‑1 (MQ=255)
             ccccAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAAt                                                                                   >  2:221293/1‑100 (MQ=255)
              cccAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATc                                                                                  >  2:3989687/1‑100 (MQ=255)
              cccAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATc                                                                                  <  2:3887079/100‑1 (MQ=255)
                cAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTAAATAATcc                                                                                 <  2:746992/99‑1 (MQ=255)
                  aGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCctc                                                                              <  1:4037833/100‑1 (MQ=255)
                  aGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCctc                                                                              <  1:698497/100‑1 (MQ=255)
                    cAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTg                                                                            <  1:3600930/100‑1 (MQ=255)
                    cAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTg                                                                            <  1:3894734/100‑1 (MQ=255)
                    cAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTg                                                                            <  2:2341247/100‑1 (MQ=255)
                    cAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTg                                                                            <  1:2595873/100‑1 (MQ=255)
                    cAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTg                                                                            <  2:4127996/100‑1 (MQ=255)
                      gggCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGaa                                                                          >  1:2276127/1‑100 (MQ=255)
                      gggCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGaa                                                                          >  1:1512643/1‑100 (MQ=255)
                       ggcggcggcGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAAt                                                                         <  2:3084884/100‑1 (MQ=255)
                         cggcggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAg                                                                       <  1:3645790/100‑1 (MQ=255)
                          ggcggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGt                                                                      >  2:2690755/1‑100 (MQ=255)
                           gcggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGtt                                                                     <  1:1339916/100‑1 (MQ=255)
                            cggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTa                                                                    <  1:1613354/100‑1 (MQ=255)
                             ggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTAt                                                                   <  2:1369082/100‑1 (MQ=255)
                             ggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTAt                                                                   <  1:3310536/100‑1 (MQ=255)
                              gcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTAtt                                                                  <  2:600878/100‑1 (MQ=255)
                               cgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTg                                                                 <  2:42935/100‑1 (MQ=255)
                               cgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTg                                                                 <  2:3569614/100‑1 (MQ=255)
                                gTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGa                                                                <  2:4039719/100‑1 (MQ=255)
                                  cTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAg                                                              >  2:1708466/1‑100 (MQ=255)
                                  cTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAg                                                              >  2:1107483/1‑100 (MQ=255)
                                   tttAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTTCCCGATACGCTAATAAAAAAAAGATTTAAAAATCCCTCTGAATAGTTATTGAAGc                                                             <  2:920538/100‑1 (MQ=255)
                                   tttAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGc                                                             >  2:3849973/1‑100 (MQ=255)
                                      aaGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTTGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGaa                                                               >  2:2041780/1‑95 (MQ=255)
                                      aaGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAg                                                          <  2:3082824/100‑1 (MQ=255)
                                       aGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGc                                                         <  2:2902950/100‑1 (MQ=255)
                                       aGATGTAACTTATCGCTGCGTTTTAATTCTTCCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGc                                                         >  2:1508871/1‑100 (MQ=255)
                                        gATGTAAGTTATCGTTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGcc                                                        >  1:2187628/1‑100 (MQ=255)
                                        gATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGcc                                                        >  2:599093/1‑100 (MQ=255)
                                         aTGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATATGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCg                                                       >  2:2260585/1‑100 (MQ=255)
                                            tAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTc                                                    >  1:100337/1‑100 (MQ=255)
                                             aaCTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCa                                                   >  1:1961317/1‑100 (MQ=255)
                                             aaCTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGATCa                                                   >  2:1514891/1‑100 (MQ=255)
                                                 tATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAAtac                                               <  2:1679174/100‑1 (MQ=255)
                                                 tATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAAtac                                               <  2:4096666/100‑1 (MQ=255)
                                                     gCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTaca                                           >  1:2720871/1‑100 (MQ=255)
                                                      cTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTacac                                          <  2:3966023/100‑1 (MQ=255)
                                                      cTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTacac                                          >  1:2936634/1‑100 (MQ=255)
                                                      cTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTacac                                          >  1:465732/1‑100 (MQ=255)
                                                      cTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCACTCAATACTacac                                          >  1:665502/1‑100 (MQ=255)
                                                       tGCGTTTCACTTCTTGCAGCATCTGCCTGTTACGCTCTTCGTAACGTGTGGCCATAAGCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAAtact                                              <  2:1974962/95‑1 (MQ=255)
                                                         cGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACttt                                       <  2:3867837/100‑1 (MQ=255)
                                                         cGTTTTAATTCTTGCAGCAGCAGCCCGATAAGCTCTTCGCAACGTGTTTCCATAATCCCTCCGAATAGTTAATGAAGAGAGTCCCCAAATCCTACACttt                                       >  2:1104880/1‑100 (MQ=255)
                                                          gTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACtttt                                      <  1:2684043/100‑1 (MQ=255)
                                                           ttttAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACttttt                                     >  1:290510/1‑100 (MQ=255)
                                                           ttttAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACttttt                                     <  1:1036744/100‑1 (MQ=255)
                                                           ttttAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACttttt                                     <  2:3099514/100‑1 (MQ=255)
                                                           ttttAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACttttt                                     <  1:2792731/100‑1 (MQ=255)
                                                           ttttAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACttttt                                     >  2:1647141/1‑100 (MQ=255)
                                                            tttAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACtttttn                                    <  1:3053731/100‑2 (MQ=255)
                                                             ttAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACCACACTTTTTAg                                   <  2:2174786/100‑1 (MQ=255)
                                                              tAAGTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGc                                  <  2:3232136/100‑1 (MQ=255)
                                                               acttcgttcagcagCTGCCCGATAGGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCa                                 <  2:3277914/92‑1 (MQ=255)
                                                               aaTTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCa                                 >  2:2643124/1‑100 (MQ=255)
                                                               aaTTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCa                                 >  1:551651/1‑100 (MQ=255)
                                                                aTTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCag                                >  2:2879658/1‑100 (MQ=255)
                                                                  tCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCagag                              <  1:2099281/100‑1 (MQ=255)
                                                                   cTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAgaga                             <  1:1809791/100‑1 (MQ=255)
                                                                   cTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAgaga                             >  1:175589/1‑100 (MQ=255)
                                                                   cTTGCAGCAGCTGCCCGATACGCGCTTCGCAGCGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCag                                <  2:3933744/97‑1 (MQ=255)
                                                                     tGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTCGCAGAGATc                           >  1:2877519/1‑100 (MQ=255)
                                                                     tGCAGCAGCTGCCCGATACGATCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATc                           >  2:2553537/1‑100 (MQ=255)
                                                                      gcaggaGCTGCCCGTTACGCGCGTCGCAGCGTGGTGCCATAATCCGTGTGAATATTGATTGAAGCGACCCGCTCAATACTACACTTGTTAGCAGAGATc                           <  2:4118830/99‑1 (MQ=255)
                                                                      gcagcaGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCa                          >  2:3756372/1‑100 (MQ=255)
                                                                      gcagcaGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCa                          >  1:2331066/1‑100 (MQ=255)
                                                                       cagcagCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAg                         <  2:3155719/100‑1 (MQ=255)
                                                                        agcagcTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGt                        >  1:1654682/1‑100 (MQ=255)
                                                                         tctgcTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTc                       <  1:2371447/97‑1 (MQ=255)
                                                                          cagcTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCa                      >  2:2790294/1‑100 (MQ=255)
                                                                          cagcTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCa                      <  1:1415026/100‑1 (MQ=255)
                                                                           agcTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCAc                     <  2:2509731/100‑1 (MQ=255)
                                                                            gcTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCaga                    >  2:2371432/1‑98 (MQ=255)
                                                                             cTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGc                   >  1:2376306/1‑100 (MQ=255)
                                                                              tGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCa                  >  1:2842339/1‑100 (MQ=255)
                                                                                 ccGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATTGTTATTGAAGCGCGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCAccc               <  2:1452091/100‑1 (MQ=255)
                                                                                 ccGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGc                   >  2:3371073/1‑96 (MQ=255)
                                                                                 ccGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACCACACTTTTTAGCAGAGATCAGTCACGCAccc               <  2:2091737/100‑1 (MQ=255)
                                                                                    aTACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGc            <  1:3297181/100‑1 (MQ=255)
                                                                                      aCGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATACTTATTGAAGCCAGCCGCTCAATACTACTCTTTTTAGCAGAGATCAGTCCCGCACCCAGcc           <  1:3962999/99‑1 (MQ=255)
                                                                                       cGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCtt         >  2:3611710/1‑100 (MQ=255)
                                                                                          tctTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGc      <  1:2203960/100‑1 (MQ=255)
                                                                                          tctTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGc      <  1:3171154/100‑1 (MQ=255)
                                                                                             tCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGa   >  1:4117189/1‑100 (MQ=255)
                                                                                                 |                                                                                                
GTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGGAT  >  NC_000913/882625‑882818

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: