Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,534,957 A→G S101P (TCG→CCG)  pdxK ← pyridoxal/pyridoxine/pyridoxamine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,534,9570AG99.6% 769.7 / NA 248S101P (TCG→CCG) pdxKpyridoxal/pyridoxine/pyridoxamine kinase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (120/127);  total (121/127)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 252 total aligned reads displayed.

ATCGCCAATCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATT  >  NC_000913/2534868‑2535046
                                                                                         |                                                                                         
atCGCCAATCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGc                                                                                        >  2:2195600/1‑93 (MQ=255)
    ccAATCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGtt                                                                                    <  2:6940408/93‑1 (MQ=255)
     cAATCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTc                                                                                   >  1:4022581/1‑93 (MQ=255)
      aaTCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTcc                                                                                  >  2:4164434/1‑93 (MQ=255)
      aaTCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTcc                                                                                  >  1:2024825/1‑93 (MQ=255)
       aTCACCGGTTCGACCTTGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAGGGATTTTGCTTTGCGGTGCCGTTccc                                                                                 <  2:7048065/93‑1 (MQ=255)
       aTCACCGGATCTGCCATGATCAATAGGTCATGATGGTCTTTTGGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTTCCGTTccc                                                                                 <  2:5632307/93‑1 (MQ=255)
       aTCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTccc                                                                                 >  2:6321094/1‑93 (MQ=255)
       aTCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTccc                                                                                 >  1:5483233/1‑93 (MQ=255)
        tCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCa                                                                                <  2:2860680/93‑1 (MQ=255)
        tCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCa                                                                                <  1:3883830/93‑1 (MQ=255)
        tCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCa                                                                                >  1:6087628/1‑93 (MQ=255)
        tCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCa                                                                                >  2:1401590/1‑93 (MQ=255)
          aCCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCata                                                                              <  2:1714680/93‑1 (MQ=255)
            cGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCAtata                                                                            <  1:5817778/93‑1 (MQ=255)
             ggATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAg                                                                           >  2:1090524/1‑93 (MQ=255)
             ggATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAg                                                                           <  1:7887973/93‑1 (MQ=255)
               aTCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGcc                                                                         <  1:5293670/93‑1 (MQ=255)
                tCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGccc                                                                        <  2:4398597/93‑1 (MQ=255)
                 cGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCg                                                                       <  2:3213497/93‑1 (MQ=255)
                  gACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGt                                                                      >  1:1699463/1‑93 (MQ=255)
                   aCCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGGGGTGCCGTTCCCATATAGCGCGTg                                                                     >  1:6665333/1‑93 (MQ=255)
                   aCCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTg                                                                     >  1:2338002/1‑93 (MQ=255)
                   aCCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTg                                                                     >  2:7313651/1‑93 (MQ=255)
                      aTGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGtt                                                                  >  1:833351/1‑93 (MQ=255)
                       tGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTa                                                                 >  2:1516156/1‑93 (MQ=255)
                         aTCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCGCGTGGTTACa                                                               >  1:6319526/1‑93 (MQ=255)
                         aTCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACa                                                               <  2:4106410/93‑1 (MQ=255)
                          tCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAg                                                              >  2:7518745/1‑93 (MQ=255)
                            aaTAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCa                                                            >  1:7107031/1‑93 (MQ=255)
                            aaTAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCa                                                            >  2:3002011/1‑93 (MQ=255)
                             aTAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCAc                                                           <  2:1770411/93‑1 (MQ=255)
                             aTAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCAc                                                           <  2:5831620/93‑1 (MQ=255)
                                 gTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAg                                                       >  2:6122524/1‑93 (MQ=255)
                                  tCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGt                                                      <  1:4779874/93‑1 (MQ=255)
                                   cAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGtt                                                     <  2:6051225/93‑1 (MQ=255)
                                    aGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTg                                                    <  1:1233570/93‑1 (MQ=255)
                                     ggATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTgg                                                   >  1:4982300/1‑93 (MQ=255)
                                      gATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgc                                                  <  2:2318726/93‑1 (MQ=255)
                                      gATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgc                                                  >  2:379144/1‑93 (MQ=255)
                                        tGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgc                                                <  2:5964162/93‑1 (MQ=255)
                                           tCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgcagc                                             <  1:3392324/93‑1 (MQ=255)
                                           tCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgcagc                                             <  1:3694765/93‑1 (MQ=255)
                                           tCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgcagc                                             <  1:7798020/93‑1 (MQ=255)
                                              ttGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgcagcgca                                          <  1:5369682/93‑1 (MQ=255)
                                              ttGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgcagcgca                                          >  2:4597987/1‑93 (MQ=255)
                                               tGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCAt                                         <  1:2961691/93‑1 (MQ=255)
                                                gCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATc                                        <  1:1463322/93‑1 (MQ=255)
                                                 cGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCa                                       <  2:5151314/93‑1 (MQ=255)
                                                    aGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGc                                    >  1:1117287/1‑93 (MQ=255)
                                                    aGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGc                                    <  1:3178724/93‑1 (MQ=255)
                                                        cAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCt                                <  2:7390560/93‑1 (MQ=255)
                                                         aGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTg                               >  1:2567962/1‑93 (MQ=255)
                                                         aGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTg                               >  1:6926286/1‑93 (MQ=255)
                                                         aGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTg                               >  1:6137272/1‑93 (MQ=255)
                                                          gTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGa                              <  2:4762724/93‑1 (MQ=255)
                                                          gTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGa                              <  1:3490917/93‑1 (MQ=255)
                                                          gTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAACCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGa                              >  2:3310559/1‑93 (MQ=255)
                                                           tCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGaa                             <  1:4790547/93‑1 (MQ=255)
                                                            cAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAg                            <  1:4687120/93‑1 (MQ=255)
                                                             aGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAgc                           >  2:1475438/1‑93 (MQ=255)
                                                              gCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAgcg                          >  2:3728480/1‑93 (MQ=255)
                                                              gCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAgcg                          >  2:6384075/1‑93 (MQ=255)
                                                               ccACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAgcgc                         >  2:4906547/1‑93 (MQ=255)
                                                                cACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGcgca                        <  2:7777745/93‑1 (MQ=255)
                                                                cACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGcgca                        >  1:4412771/1‑93 (MQ=255)
                                                                cACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGcgca                        >  2:44152/1‑93 (MQ=255)
                                                                cACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGcgca                        <  1:5130271/93‑1 (MQ=255)
                                                                 ccTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGcgcac                       <  2:6683556/92‑1 (MQ=255)
                                                                 aCTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGcgcac                       >  2:724940/1‑93 (MQ=255)
                                                                 aCTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGcgcac                       <  2:5793413/93‑1 (MQ=255)
                                                                 aCTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGcgcac                       >  1:3750759/1‑93 (MQ=255)
                                                                 aCTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGcgcac                       <  2:1653683/93‑1 (MQ=255)
                                                                 aCTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGcgcac                       <  1:7539080/93‑1 (MQ=255)
                                                                  cTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGcgcacg                      >  2:6869999/1‑93 (MQ=255)
                                                                     ggCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCaa                   >  2:1830829/1‑93 (MQ=255)
                                                                      gCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCaaa                  <  1:979110/93‑1 (MQ=255)
                                                                      gCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCaaa                  <  1:1788190/93‑1 (MQ=255)
                                                                      gCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCaaa                  <  1:4550356/93‑1 (MQ=255)
                                                                       cAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAAt                 >  1:7244270/1‑93 (MQ=255)
                                                                        aaGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATa                >  1:7412150/1‑93 (MQ=255)
                                                                         aGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAg               >  2:5299343/1‑93 (MQ=255)
                                                                         aGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAg               >  2:5555190/1‑93 (MQ=255)
                                                                         aGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAg               >  2:25512/1‑93 (MQ=255)
                                                                          ggATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGc              <  2:839154/93‑1 (MQ=255)
                                                                           ggttttGATTTGCGGTGCCGTTCCCATATAGCACGTGGTTACAGCACGAAGTTGGCGCAGCGCATCGCGCTCCTAAAACGCACGCTAATAGcc             >  2:1350868/3‑93 (MQ=255)
                                                                           gATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGcc             <  1:4722330/93‑1 (MQ=255)
                                                                            aTTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCg            >  2:3371854/1‑93 (MQ=255)
                                                                              tttGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCt          >  1:5234912/1‑93 (MQ=255)
                                                                               ttGATTTGCGGTGCCGATCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTa         <  1:7043363/93‑1 (MQ=255)
                                                                                 gATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTaaa       <  2:3486802/93‑1 (MQ=255)
                                                                                 gATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTaaa       <  1:6365006/93‑1 (MQ=255)
                                                                                   tttGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAAcc     >  2:2665680/1‑93 (MQ=255)
                                                                                   tttGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAAcc     >  1:2578620/1‑93 (MQ=255)
                                                                                    ttGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCa    <  1:6583290/93‑1 (MQ=255)
                                                                                     tGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCAt   <  2:7694039/93‑1 (MQ=255)
                                                                                     tGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCAt   <  1:1011353/93‑1 (MQ=255)
                                                                                      gCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCAtt  <  2:6087627/93‑1 (MQ=255)
                                                                                      gCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCAtt  <  1:5702067/93‑1 (MQ=255)
                                                                                      gCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCAtt  <  1:561823/93‑1 (MQ=255)
                                                                                         |                                                                                         
ATCGCCAATCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATT  >  NC_000913/2534868‑2535046

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