Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,070,600 C→T 6.0% P215P (CCG→CCA fbaA ← fructose‑bisphosphate aldolase, class II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,070,6000CT 6.0% 214.3 / 2.7 66P215P (CCG→CCAfbaAfructose‑bisphosphate aldolase, class II
Reads supporting (aligned to +/- strand):  ref base C (32/30);  new base T (2/2);  total (34/32)
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.38e-01

ATACGAAGTTCAGGCTGTTGTGCGGCAGGTTGTGTTTCTTGGAAACATATTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTT  >  NC_000913/3070455‑3070745
                                                                                                                                                 |                                                                                                                                                 
aTACGAAGTTCAGGCTGTTGTGCGGCAGGTTGTGTTTCTTGGAAACATATTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCgg                                                                                                                                                 >  7:22923/1‑148 (MQ=255)
aTACGAAGTTCAGGCTGTTGTGCGGCAGGTTGTGTTTCTTGGAAACATATTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCgg                                                                                                                                                 <  7:293192/148‑1 (MQ=255)
  aCGAAGTTCAGGCTGTTGTGCGGCAGGTTGTGTTTCTTGGAAACATATTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGc                                                                                                                                               >  4:236885/1‑148 (MQ=255)
       gTTCAGGCTGTTGTGCGGCAGGTTGTGTTTCTTGGAAACATATTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGAtt                                                                                                                                          >  3:158323/1‑148 (MQ=255)
              cTGTTGTGCGGCAGGTTGTGTTTCTTGGAAACATATTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTAACCGGCATGAAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCt                                                                                                                                     >  1:82990/1‑146 (MQ=255)
                    tgCGGCAGGTTGTGTTTCTTGGAAACATATTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCgg                                                                                                                             <  3:236885/148‑1 (MQ=255)
                    tgCGGCAGGTTGTGTTTCTTGGAAACATATTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCgg                                                                                                                             <  5:2114/148‑1 (MQ=255)
                    tgCGGCAGGTTGTGTTTCTTGGAAACATATTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCgg                                                                                                                             <  8:22923/148‑1 (MQ=255)
                       ggCAGGTTGTGTTTCTTGGAAACATATTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGgtgt                                                                                                                          <  3:83964/148‑1 (MQ=255)
                           ggTTGTGTTTCTTGGAAACATATTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGc                                                                                                                      <  7:30158/148‑1 (MQ=255)
                             ttGTGTTTCTTGGAAACATATTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGTGGGCTGATTTTGCTCAGTTCgg                                                                                                                             <  5:216655/139‑1 (MQ=255)
                             ttGTGTTTCTTGGAAACATATTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGTGGGCTGATTTTGCTCAGTTCgg                                                                                                                             >  6:216655/1‑139 (MQ=255)
                              tgtgTTTCTTGGAAACATATTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGt                                                                                                                    >  1:136915/1‑147 (MQ=255)
                                tgtTTCTTGGAAACATATTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTa                                                                                                                         <  2:208912/140‑1 (MQ=255)
                                tgtTTCTTGGAAACATATTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTa                                                                                                                         >  1:208912/1‑140 (MQ=255)
                                tgtTTCTTGGAAACATATTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAAt                                                                                                                 >  7:139091/1‑148 (MQ=255)
                                    tCTTGGAAACATATTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGGTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCaa                                                                                                              <  2:193811/147‑1 (MQ=255)
                                       tGGAAACATATTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTc                                                                                                          >  5:198496/1‑148 (MQ=255)
                                          aaaCATATTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTTTACGTTACCGAAGGACGCTGCGATTGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTc                                                                                                       <  2:66669/148‑1 (MQ=255)
                                           aaCATATTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGgtgt                                                                                                                          <  8:298187/128‑1 (MQ=255)
                                           aaCATATTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGgtgt                                                                                                                          >  7:298187/1‑128 (MQ=255)
                                              atatTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGc                                                                                                   <  4:158323/148‑1 (MQ=255)
                                               tatTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGTGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCt                                                                                                  >  7:327941/1‑148 (MQ=255)
                                               tatTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCt                                                                                                  >  1:163581/1‑148 (MQ=255)
                                               tatTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGCGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCt                                                                                                  >  3:246974/1‑148 (MQ=255)
                                                 ttCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGTGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTg                                                                                                 <  8:327941/147‑1 (MQ=255)
                                                    cTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGgtg                                                                                             >  1:104625/1‑148 (MQ=255)
                                                    cTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGgt                                                                                              <  7:132190/147‑1 (MQ=255)
                                                     tGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTAACCGGCTTGTAAACAACTTGTACGTTACCGAAGGCCGCTGCGAAGGTGAAACGCGGTCTGATTTTGCTCAGTTCGGTGTATGCGTAATAAACGTCTTCCGGCTGGgtg                                                                                             >  1:110985/1‑147 (MQ=255)
                                                      gagaATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTa                                                                                           >  8:193140/1‑148 (MQ=255)
                                                      gagaATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTa                                                                                           >  6:13314/1‑148 (MQ=255)
                                                       agaATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGgtgt                                                                                            <  8:139091/146‑1 (MQ=255)
                                                          aTCACGCAGGATGGTCGTAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGt                                                                                       >  2:76411/1‑148 (MQ=255)
                                                          aTCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTAc                                                                                          >  6:94077/1‑145 (MQ=255)
                                                          aTCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTAc                                                                                          <  5:94077/145‑1 (MQ=255)
                                                           tCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTg                                                                                      <  2:163581/148‑1 (MQ=255)
                                                            cACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTg                                                                                      >  4:101751/1‑147 (MQ=255)
                                                              cGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTg                                                                                      >  6:159292/1‑145 (MQ=255)
                                                                       gTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGCGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTcc                                                                           <  4:246974/147‑1 (MQ=255)
                                                                        tCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCAt                                                                         >  1:881/1‑148 (MQ=255)
                                                                        tCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCAt                                                                         <  5:65033/148‑1 (MQ=255)
                                                                           gAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTa                                                                                                                         <  7:179942/97‑1 (MQ=255)
                                                                           gAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTa                                                                                                                         >  8:179942/1‑97 (MQ=255)
                                                                               cAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCtg                                                                  >  5:271983/1‑148 (MQ=255)
                                                                               cAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCtg                                                                  >  4:37735/1‑148 (MQ=255)
                                                                               cAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCtg                                                                  >  4:183099/1‑148 (MQ=255)
                                                                                 gAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCtgtt                                                                <  6:271983/148‑1 (MQ=255)
                                                                                    ccACGTTACCCGGCTTGTAAACTCCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAAAGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTc                                                             <  2:82990/148‑1 (MQ=255)
                                                                                       cGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCAc                                                          >  8:207747/1‑148 (MQ=255)
                                                                                               ggCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGtctt                                                  >  1:280749/1‑148 (MQ=255)
                                                                                                 cTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGtcttct                                                <  6:198496/148‑1 (MQ=255)
                                                                                                   tGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGAGTCCATGTGGCTGTTGTCCACGCCGTcttctt                                               >  1:82406/1‑147 (MQ=255)
                                                                                                    gTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTcttctt                                               <  7:193140/146‑1 (MQ=255)
                                                                                                    gTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCttcttc                                              <  5:159292/147‑1 (MQ=255)
                                                                                                     tAAACACCGTGTACGTTACCGAAGGACGCAGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGACTTCCGGCTGGGTGTACAGAGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTcac                                            >  5:290789/1‑148 (MQ=255)
                                                                                                           ccGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCGGTTGTCCACGCCGTCTTCTTCACCACCgg                                      <  5:13314/148‑1 (MQ=255)
                                                                                                             gtgtACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTg                                    >  3:51297/1‑148 (MQ=255)
                                                                                                                   gTTACCGAAGGAGGCTGAGTTGGTGAATCGCGGGCTGATTTTGCTCAGTTCTGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAg                                                                                <  8:102273/98‑1 (MQ=255)
                                                                                                                    ttACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGGTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGc                                                                               >  7:102273/1‑98 (MQ=255)
                                                                                                                          aaGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGa                       >  6:108295/1‑148 (MQ=255)
                                                                                                                           aGGACGCTGCGATGGTGAAACGCGTGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCAt                                                                         <  8:163893/97‑1 (MQ=255)
                                                                                                                           aGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCAt                                                                         >  7:163893/1‑97 (MQ=255)
                                                                                                                             gACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGAttt                    <  2:881/148‑1 (MQ=255)
                                                                                                                                gCTGGGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCa                 <  1:76411/148‑1 (MQ=255)
                                                                                                                                    cGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGt             <  3:37735/148‑1 (MQ=255)
                                                                                                                                     gATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTc            <  3:183099/148‑1 (MQ=255)
                                                                                                                                               cgcgGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGc                                                                    >  4:14481/1‑82 (MQ=255)
                                                                                                                                               cgcgGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGc                                                                    <  3:14481/82‑1 (MQ=255)
                                                                                                                                               cgcgGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGAttt  <  4:51297/148‑1 (MQ=255)
                                                                                                                                               cgcgGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGAttt  <  2:136915/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
ATACGAAGTTCAGGCTGTTGTGCGGCAGGTTGTGTTTCTTGGAAACATATTCCTGAGAATCACGCAGGATGGTCGGAGTCAGAACCACGTTACCCGGCTTGTAAACACCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTT  >  NC_000913/3070455‑3070745

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