Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 53,376 G→A 16.6% P14L (CCC→CTC)  pdxA ← 4‑hydroxy‑L‑threonine phosphate dehydrogenase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091353,3760GA16.6% 456.2 / 65.3 192P14L (CCC→CTC) pdxA4‑hydroxy‑L‑threonine phosphate dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base G (83/77);  new base A (18/14);  total (101/91)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 207 total aligned reads displayed.

GCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGC  >  NC_000913/53277‑53475
                                                                                                   |                                                                                                   
 cccGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCggg                                                                                                    <  1:4124430/100‑1 (MQ=255)
     gTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCg                                                                                                >  2:1161115/1‑100 (MQ=255)
        ggTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCg                                                                                             >  1:4060076/1‑100 (MQ=255)
        ggTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCg                                                                                             <  1:2165964/100‑1 (MQ=255)
         gTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCgg                                                                                            >  2:2983782/1‑100 (MQ=255)
          tAAGGAGAGTGGCATCGGCACAAACAACCAGTGCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGAGCCCGCCggg                                                                                           <  1:1189140/100‑1 (MQ=255)
            aGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAg                                                                                         <  1:627139/100‑1 (MQ=255)
             ggAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGt                                                                                        >  1:908711/1‑100 (MQ=255)
             ggAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACGCACGCTGTGCAAGCTGGACTACTAAGTCCGGGCCAATCCCGCCGGGCTCGCCGGGAg                                                                                         >  2:339229/1‑99 (MQ=255)
              gagaGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTg                                                                                       >  1:2976718/1‑100 (MQ=255)
              gagaGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTg                                                                                       >  1:3875048/1‑100 (MQ=255)
              gagaGTGGCAGCGGCACAAACAACAAGTTCGAGCGGCCGCTCACGCTGTGCAAGCTGCACATCGAAGTCCGGGCCAATCCGGGCGAGCCCGCCGGGatt                                                                                        <  2:1573901/99‑3 (MQ=255)
               agagTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGAGCCCGCCGGGATTgg                                                                                      <  1:452206/100‑2 (MQ=255)
               agagTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGATGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTg                                                                                       >  2:3026865/1‑99 (MQ=255)
                 agTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATc                                                                                    <  2:3841228/100‑1 (MQ=255)
                  gTGGCATCGGCACAACCAACCCGTTCGACCGGCCACTCACGCTGTGCAAGCGGGACAACTAAGTCCGGGCCAATCCCGGCGGGCCCGCCGGGAGTGATc                                                                                    <  2:2307644/99‑1 (MQ=255)
                  gTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATca                                                                                   >  1:1263718/1‑100 (MQ=255)
                   tGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATcac                                                                                  >  2:1405406/1‑100 (MQ=255)
                     gCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCacaa                                                                                >  2:3084357/1‑100 (MQ=255)
                       aTCGGCACAAACACCCAGTTCGACCGGCCACGCACGCTGTGCAAGCTGGGCAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCacaac                                                                               <  2:3646680/99‑1 (MQ=255)
                         cGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACg                                                                            <  2:693420/100‑1 (MQ=255)
                         cGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACg                                                                            <  1:2146027/100‑1 (MQ=255)
                          ggCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGt                                                                           >  1:3574556/1‑100 (MQ=255)
                          ggCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGt                                                                           >  2:3145357/1‑100 (MQ=255)
                          ggCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGt                                                                           >  2:65110/1‑100 (MQ=255)
                          ggCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGt                                                                           >  1:1960344/1‑100 (MQ=255)
                             acaaacaaCCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTgg                                                                        <  2:2195208/100‑1 (MQ=255)
                              caaccaacCAGTTCGACCGGCCCCTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCAAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACTTTggg                                                                       <  2:3497887/100‑1 (MQ=255)
                              caaacaaCCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTggg                                                                       >  1:545456/1‑100 (MQ=255)
                              caaacaaCCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTggg                                                                       >  2:586355/1‑100 (MQ=255)
                              caaacaaCCAGTCCGACGGGCCACTCAGGCTGTGCAGGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACCACACGTTgg                                                                        <  2:1522885/99‑1 (MQ=255)
                               aaacaaCCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGt                                                                      <  1:3119787/100‑1 (MQ=255)
                                   aacCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTa                                                                  >  2:2814278/1‑100 (MQ=255)
                                    acCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTaa                                                                 >  1:1359753/1‑100 (MQ=255)
                                    acCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTaa                                                                 <  1:1607618/100‑1 (MQ=255)
                                     ccAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAAc                                                                <  1:1688778/100‑1 (MQ=255)
                                      cAGTTCGACCGGTCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAAcc                                                               >  1:3660733/1‑100 (MQ=255)
                                       aGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCa                                                              >  1:2528982/1‑100 (MQ=255)
                                        gTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCAt                                                             >  2:119978/1‑100 (MQ=255)
                                        gTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGGTTTAACCAt                                                             >  2:2167538/1‑100 (MQ=255)
                                         ttCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCAtt                                                            <  2:3198229/100‑1 (MQ=255)
                                          tCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTa                                                           >  2:3060391/1‑100 (MQ=255)
                                           cGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAg                                                          >  1:3974455/1‑100 (MQ=255)
                                           cGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAg                                                          <  1:3386291/100‑1 (MQ=255)
                                           cGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAg                                                          >  1:1008196/1‑100 (MQ=255)
                                                 gCCACTCACGCTGTGCAAGCTGGACCACTAAGTCCGGGCCAATATCGGAGAGCTCGCCTGGAGTGATCACAACACGTGGGGTTATACCCATTAGTGGCt                                                     <  2:521988/99‑1 (MQ=255)
                                                      tCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGGGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGa                                                >  2:2618262/1‑99 (MQ=255)
                                                       cACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGAtt                                              <  2:2048201/100‑1 (MQ=255)
                                                       cACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCGCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGAtt                                              >  2:2099861/1‑100 (MQ=255)
                                                        aCGCTGTGCAAGCTGGCCAACTAAGTCCGGGCCCATCCCGGCGGGCTCGCCGGGAGTGATCACACCACGTTGGTTTTTAACCATGAGCTGCTCAGGAttt                                             >  2:1712534/1‑100 (MQ=255)
                                                        aCGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGAttt                                             >  1:3442835/1‑100 (MQ=255)
                                                         cGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACGCGCTGGGCTTTAACCATCAGGTACTCAGGAtttt                                            >  2:612049/1‑100 (MQ=255)
                                                         cGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGAtttt                                            >  2:2864397/1‑100 (MQ=255)
                                                         cGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGAtttt                                            >  1:3908081/1‑100 (MQ=255)
                                                          gCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTa                                           <  1:2564612/100‑1 (MQ=255)
                                                            tgtgCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAAc                                         <  1:738215/100‑1 (MQ=255)
                                                            tgtgCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAAc                                         <  1:3239354/100‑1 (MQ=255)
                                                            tgtgCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAAc                                         >  1:4157228/1‑100 (MQ=255)
                                                            tgtgCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGAACACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAAc                                         >  1:1751810/1‑100 (MQ=255)
                                                              tgCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGt                                       <  2:1252504/100‑1 (MQ=255)
                                                               gCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTa                                      <  2:3357698/100‑1 (MQ=255)
                                                               gCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTa                                      >  1:3802472/1‑100 (MQ=255)
                                                                 aaGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAgg                                    <  2:918551/100‑1 (MQ=255)
                                                                 aaGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAgg                                    <  1:2558920/100‑1 (MQ=255)
                                                                  aGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGTGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAga                                    <  2:2316393/99‑2 (MQ=255)
                                                                  aGCTGGACAACTAAGTCCGGGCCAAACCCGGCGGGTTAGCCCGGAGTGATCACACCACGTTGCGTTTTAACCATAAGTTGCTCAGGATTTTAACg                                        >  2:2769926/1‑95 (MQ=255)
                                                                      ggACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTg                               <  2:848496/100‑1 (MQ=255)
                                                                        aCAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGc                             <  2:3974428/100‑1 (MQ=255)
                                                                          aaCTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCa                            <  2:253402/99‑1 (MQ=255)
                                                                           aCTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACg                          <  2:2253384/100‑1 (MQ=255)
                                                                              aaGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACgttg                       <  1:833395/100‑1 (MQ=255)
                                                                                gTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACgttgtt                     >  1:2679840/1‑100 (MQ=255)
                                                                                gTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACgttgtt                     >  2:439425/1‑100 (MQ=255)
                                                                                  ccGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTcc                   <  1:972436/100‑1 (MQ=255)
                                                                                  ccGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTcc                   <  2:3802448/100‑1 (MQ=255)
                                                                                  ccGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTcc                   >  1:3831275/1‑100 (MQ=255)
                                                                                   gtgatcactCCCGGCGAGCCCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCt                  >  1:530976/9‑100 (MQ=38)
                                                                                    tgatcactCCCGGCGAGCCCGCCGGGATTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTg                 <  1:2693234/93‑1 (MQ=37)
                                                                                    gggCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTg                 <  1:1928500/100‑1 (MQ=255)
                                                                                    gggCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTg                 <  1:2302787/100‑1 (MQ=255)
                                                                                      atcactCCCGGCGAGCCCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCa               >  2:3859657/6‑100 (MQ=255)
                                                                                       tccgatctggcggGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCAt              <  2:2225912/92‑1 (MQ=255)
                                                                                       tcactCCCGGCGAGCCCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCAt              >  2:4168343/5‑100 (MQ=255)
                                                                                       ccAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCAt              <  1:3145371/100‑1 (MQ=255)
                                                                                         actCCCGGCGAGCCCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATcc            >  2:2382875/3‑100 (MQ=255)
                                                                                         aaTCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATcc            <  2:454551/100‑1 (MQ=255)
                                                                                         aaTCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATcc            >  1:327709/1‑100 (MQ=255)
                                                                                         aaTCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTAGTTCCTGCATcc            <  1:2757963/100‑1 (MQ=255)
                                                                                          ctCCCGGCGAGCCCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCa           >  2:1885828/2‑100 (MQ=255)
                                                                                          ctCCCGGCGAGCCCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCa           >  1:2101520/2‑100 (MQ=255)
                                                                                          ctCCCGGCGAGCCCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCa           >  1:1384823/2‑100 (MQ=255)
                                                                                          ctCCCGGCGAGCCCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCa           >  2:1364304/2‑100 (MQ=255)
                                                                                          ctCCCGGCGAGCCCGCCGGGAGTGATCACAACACGTTGGGTTTTAACAATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCa           >  1:2017230/2‑100 (MQ=255)
                                                                                            cccGGCGAGCCCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGc         >  1:3864795/1‑100 (MQ=255)
                                                                                            cccGCCGAGCCCGCCAGGAGCGATCACAACACGTTGGGTTTTAACCATTGGTTGCTCAGGAGATTAACGTAAGAGCTGGCACGTTGTTCCTGCATCCAg          >  2:208060/1‑99 (MQ=255)
                                                                                               ggcggGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCAtagatcgga      >  1:2225917/1‑92 (MQ=255)
                                                                                               ggcggGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTtg      <  1:305453/100‑1 (MQ=255)
                                                                                                 cggGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTtgct    <  1:782324/100‑1 (MQ=255)
                                                                                                  gggCTCGCCGGGAGTGATCACAACACGTGGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTtgctg   <  2:2918287/100‑1 (MQ=255)
                                                                                                   ggCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTtgctgc  >  1:156474/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GCCCGGTTGGTAAGGAGAGTGGCATCGGCACAAACAACCAGTTCGACCGGCCACTCACGCTGTGCAAGCTGGACAACTAAGTCCGGGCCAATCCCGGCGGGCTCGCCGGGAGTGATCACAACACGTTGGGTTTTAACCATTAGTTGCTCAGGATTTTAACGTAGGCGCTGGCACGTTGTTCCTGCATCCAGCTTGCTGC  >  NC_000913/53277‑53475

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: