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.1% 364.4 / ‑4.6 109R5C (CGC→TGC) deoRdeoxyribose‑5‑phosphate‑inducible deoxyribose operon transcriptional repressor; repressor of nupG and tsx
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (56/52);  total (57/52)
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.12e-01

Only 100 of 109 total aligned reads displayed.

GTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGG  >  NC_000913/882625‑882816
                                                                                                 |                                                                                              
gTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAAc                                                                                              <  1:655930/100‑1 (MQ=255)
 tCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACg                                                                                             <  2:45125/100‑1 (MQ=255)
  cATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACgt                                                                                            >  2:1405187/1‑100 (MQ=255)
    tctcCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACgtgt                                                                                          >  1:3655497/1‑100 (MQ=255)
       ccGAAACCCCAAGCAGGGCGGCTGCGTCTTTGAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTc                                                                                       <  2:2787432/100‑1 (MQ=255)
       ccGAAACCCCAAGCAGGGCGGCGGCGTCTTTAGGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTc                                                                                       <  2:3276668/100‑1 (MQ=255)
       ccGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGTGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTc                                                                                       >  2:3971125/1‑100 (MQ=255)
        cGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTcc                                                                                      >  2:1731108/1‑100 (MQ=255)
        cGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTcc                                                                                      >  2:2663226/1‑100 (MQ=255)
        cGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTcc                                                                                      >  2:3021620/1‑100 (MQ=255)
        cGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTc                                                                                       <  2:3606330/99‑1 (MQ=255)
         gAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCa                                                                                     >  1:3479198/1‑100 (MQ=255)
           aaCCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCAt                                                                                    <  2:1263708/99‑1 (MQ=255)
             ccccAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAAt                                                                                 <  1:1518052/100‑1 (MQ=255)
             ccccAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAAt                                                                                 >  1:1002416/1‑100 (MQ=255)
              cccAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATc                                                                                <  2:11896/100‑1 (MQ=255)
              cccAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATc                                                                                >  2:1476379/1‑100 (MQ=255)
              cccAAGCAGGGCGGCGGCGTCTTTAAGAAGTAACTTATCGCTGCGTTCTAATTCTTGCAGCAGCTGCCCGATACGCACTTCTCAACGTGTTTCCATAAt                                                                                 >  1:72132/1‑99 (MQ=255)
                cAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATccc                                                                              >  1:618413/1‑100 (MQ=255)
                  aGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCctc                                                                            <  2:1154194/100‑1 (MQ=255)
                    cAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTg                                                                          <  1:626770/100‑1 (MQ=255)
                    cAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTg                                                                          <  1:3021639/100‑1 (MQ=255)
                     aGGGCTGCTGCGTCGTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCTATACGCTCGTCGCAACGTGTTTCCATAATCCctct                                                                           <  2:2374144/98‑1 (MQ=255)
                     aGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGa                                                                         <  1:2039594/100‑1 (MQ=255)
                       ggcggcggcGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGGTACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAAt                                                                       >  2:2246312/1‑100 (MQ=255)
                       ggcggcggcGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAAt                                                                       >  1:2601268/1‑100 (MQ=255)
                         cggcggcgTCTTTAAGATGTAACTTAGCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAg                                                                     <  1:615050/100‑1 (MQ=255)
                          ggcggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTACTTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGt                                                                    >  2:2548935/1‑100 (MQ=255)
                          ggcggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGt                                                                    <  1:1794072/100‑1 (MQ=255)
                           gcggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGtt                                                                   >  2:683152/1‑100 (MQ=255)
                            cggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTa                                                                  <  2:3479685/100‑1 (MQ=255)
                              gcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTAtt                                                                <  1:3336667/100‑1 (MQ=255)
                              gcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTAtt                                                                >  1:246572/1‑100 (MQ=255)
                                gTCTTTAAGATGTACCTTATCGCTGCGTTTTAATACTTGCAGCAGCGGCCCGATACGCTCTTCGCAACGTGTGTCCATAAACCCGCTGAATAGTTATTg                                                               <  2:3906804/99‑1 (MQ=255)
                                gTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGa                                                              >  1:2697380/1‑100 (MQ=255)
                                 tCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGGTATTGaa                                                             <  1:2163877/100‑1 (MQ=255)
                                  cTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAg                                                            <  2:1387998/100‑1 (MQ=255)
                                  cTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAg                                                            <  1:1647544/100‑1 (MQ=255)
                                       aGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGc                                                       <  2:4018898/100‑1 (MQ=255)
                                        gATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGcc                                                      <  2:1002410/100‑1 (MQ=255)
                                        gATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGcc                                                      >  2:1348956/1‑100 (MQ=255)
                                         aTGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCg                                                     <  2:3051432/100‑1 (MQ=255)
                                          tggaaCTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGc                                                    <  2:636378/97‑1 (MQ=255)
                                          tGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGc                                                    <  2:1076560/100‑1 (MQ=255)
                                          tGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGACCAAGACTCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGc                                                    <  2:202329/100‑1 (MQ=255)
                                           gTAACTTATCGCTGCGTTTTAATTCTTGCAGCACCTGCTCGATACGCTCTTCGCAACTTGATTCCATAATCCCTCTCAATAGTTTTTGaaccaaacagc                                                    >  2:835972/1‑90 (MQ=255)
                                            tAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTc                                                  >  1:4056936/1‑100 (MQ=255)
                                            tAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTc                                                  <  1:1978769/100‑1 (MQ=255)
                                            tAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATACCTCTGACTAGTTATTGAAGCGAGCCGCTc                                                  >  2:3528655/1‑100 (MQ=255)
                                             aaCTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCa                                                 >  2:176549/1‑100 (MQ=255)
                                               cTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAAt                                               <  1:2377130/100‑1 (MQ=255)
                                                   tCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAAtacta                                           >  1:719161/1‑100 (MQ=255)
                                                     gctgctTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCACCGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAAtactac                                          <  2:2293479/99‑1 (MQ=255)
                                                     gCTGCGTTTTAATTCTTGCAGTAGCTGCCCGAGACGCTCTTCGCAACGTGTTTCCATACTCCCTCTGAATAGCTATTGAAGCGACCCGCGCAATACTaca                                         >  2:2383094/1‑100 (MQ=255)
                                                     gCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTaca                                         >  2:3605099/1‑100 (MQ=255)
                                                      cTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAAtactac                                          <  2:1450752/98‑1 (MQ=255)
                                                      cTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTacac                                        >  2:1609840/1‑100 (MQ=255)
                                                       tGCGTTTTAATTCTTGCAGCAGCTGCCCTATTCGCTCTTCGCAACGTGTTTCCATATTCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAAtactac                                          <  2:893524/97‑1 (MQ=255)
                                                        gCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACtt                                      >  1:3731986/1‑100 (MQ=255)
                                                         cGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACttt                                     <  1:3855357/100‑1 (MQ=255)
                                                         cGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACttt                                     <  2:1878556/100‑1 (MQ=255)
                                                           ttttAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACttttt                                   >  1:1049635/1‑100 (MQ=255)
                                                           tgttAATTCTTGCAGCAACTGCCCGATACGCTCGTCGCAACGTGTTTCCAGATTCCCTATTAATAATTACTGAAGCTAGCCGCTCAATCCTACAAttt                                     >  2:308582/3‑98 (MQ=255)
                                                            tttAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTa                                  >  1:3590950/1‑100 (MQ=255)
                                                              tAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGc                                <  1:828186/100‑1 (MQ=255)
                                                              tAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGc                                >  1:2039259/1‑100 (MQ=255)
                                                               aaTTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCa                               >  1:327184/1‑100 (MQ=255)
                                                                aTTCTTGCAGCAGCTGCCCGTTACGCTCTTCGCAACGTGTTTCCATAATCCTTTTTTTTTTTTAGAGACGGCGGCCACCCAATACTACACTTTTTAGCag                              <  2:77658/100‑1 (MQ=255)
                                                                aTTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCag                              >  2:4076133/1‑100 (MQ=255)
                                                                 ttCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCaga                             >  2:1241914/1‑100 (MQ=255)
                                                                 ttCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCaga                             <  2:72132/100‑1 (MQ=255)
                                                                  tCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCagag                            >  1:2229847/1‑100 (MQ=255)
                                                                   cTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAgaga                           >  1:1173069/1‑100 (MQ=255)
                                                                    ttGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGAt                          >  2:2102962/1‑100 (MQ=255)
                                                                     tGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATc                         <  1:1728807/100‑1 (MQ=255)
                                                                     tGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATc                         <  1:3470705/100‑1 (MQ=255)
                                                                     tGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATc                         >  2:1918566/1‑100 (MQ=255)
                                                                      gcagcaGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCa                        <  1:1248214/100‑1 (MQ=255)
                                                                      gcagcaGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCa                        >  1:3492201/1‑100 (MQ=255)
                                                                       cagcagCTGCCCGATACGCTCTTCGCAACGTGTTTCCTTCATCCCTCTGAATAGTTAGCGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAg                       <  2:1891093/100‑1 (MQ=255)
                                                                          cagcTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCa                    <  2:1455124/100‑1 (MQ=255)
                                                                           agcTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCAc                   <  1:2284584/100‑1 (MQ=255)
                                                                            gcTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACg                  <  2:2708500/100‑1 (MQ=255)
                                                                            gcTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTCAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACg                  <  1:2384068/100‑1 (MQ=255)
                                                                               gCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTCTTGAAGCGAGCCGCTCAATACTACACTTTTTAACAGAGATCAGTCACGCAc               >  2:2615814/1‑100 (MQ=255)
                                                                               gCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCAc               <  2:273844/100‑1 (MQ=255)
                                                                               gCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCAc               >  2:1155487/1‑100 (MQ=255)
                                                                               gCCCGATACGCTCTTCGCAACGTGTATCCATAATAACTCTGAATAGTTATTGAAGCAAGCCGCTCAATATTAAAATTTTTAGCAGAGATCAGTAAc                   >  2:1395893/1‑96 (MQ=255)
                                                                                 ccGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCAccc             >  1:754311/1‑100 (MQ=255)
                                                                                  cGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCa            >  1:4032422/1‑100 (MQ=255)
                                                                                  cGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCa            <  1:67314/100‑1 (MQ=255)
                                                                                       cGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCtt       >  1:653706/1‑100 (MQ=255)
                                                                                       cGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCtt       >  2:3727193/1‑100 (MQ=255)
                                                                                       cGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCtt       >  1:451662/1‑100 (MQ=255)
                                                                                       cGCTCTTCGCAACGTGTTTCCATAATCCCTCTAAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCATCCtt       >  1:3924797/1‑100 (MQ=255)
                                                                                       cGCTCTTCGCAACGTGTTTCCATAATACGTGTGATGAGTTATTGATGCGAGCAGTTCAACACTACACTTTGTAGCAGAGATCAGTCACGCACCCAGCCtt       >  2:2586466/1‑100 (MQ=255)
                                                                                        gCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTGTTAGCAGAGAGCAGTCACGCACCCAGCCttt      >  2:3102384/1‑100 (MQ=255)
                                                                                           ctTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCg   <  1:1807823/100‑1 (MQ=255)
                                                                                            ttCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCgg  <  1:2201099/100‑1 (MQ=255)
                                                                                            ttCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCgg  >  2:3683704/1‑100 (MQ=255)
                                                                                                 |                                                                                              
GTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGG  >  NC_000913/882625‑882816

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