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.5% 705.0 / ‑4.8 218S101P (TCG→CCG) pdxKpyridoxal/pyridoxine/pyridoxamine kinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (113/104);  minor base T (1/0);  total (114/104)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 223 total aligned reads displayed.

TATCGCCAATCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATT  >  NC_000913/2534867‑2535046
                                                                                          |                                                                                         
tatCGCCAATCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTg                                                                                         >  1:6703860/1‑93 (MQ=255)
   cGCCAATCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCg                                                                                      >  2:5638020/1‑93 (MQ=255)
      cAATCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTc                                                                                   >  2:1576669/1‑93 (MQ=255)
      cAATCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTc                                                                                   >  1:3219384/1‑93 (MQ=255)
      cAATCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTc                                                                                   <  1:646468/93‑1 (MQ=255)
      cAATCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGCTc                                                                                   >  1:4554358/1‑93 (MQ=255)
      cAATCACCGCATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTc                                                                                   <  1:6602554/93‑1 (MQ=255)
       aaTCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTcc                                                                                  >  1:380918/1‑93 (MQ=255)
       aaTCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTcc                                                                                  >  1:5600073/1‑93 (MQ=255)
       aaTCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTcc                                                                                  >  1:3382821/1‑93 (MQ=255)
       aaTCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTcc                                                                                  <  1:2294077/93‑1 (MQ=255)
       aaTCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTcc                                                                                  >  1:23136/1‑93 (MQ=255)
        aTCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTccc                                                                                 <  2:1527510/93‑1 (MQ=255)
        aTCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTccc                                                                                 >  2:391549/1‑93 (MQ=255)
            ccGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCatat                                                                             >  2:5475204/1‑93 (MQ=255)
             cGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCAtata                                                                            >  1:1370631/1‑93 (MQ=255)
             cGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCAtata                                                                            <  1:111466/93‑1 (MQ=255)
               gATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGc                                                                          <  2:3777467/93‑1 (MQ=255)
                aTCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGcc                                                                         >  2:1115331/1‑93 (MQ=255)
                   gACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGt                                                                      <  1:6175966/93‑1 (MQ=255)
                    aCCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTg                                                                     >  1:2690303/1‑93 (MQ=255)
                     ccATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCTGTgg                                                                    >  2:4227676/1‑93 (MQ=255)
                     ccATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTgg                                                                    <  2:675589/93‑1 (MQ=255)
                      cATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGt                                                                   >  1:2182671/1‑93 (MQ=255)
                       aTGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGtt                                                                  >  2:5900173/1‑93 (MQ=255)
                       aTGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGtt                                                                  >  2:4525255/1‑93 (MQ=255)
                         gATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGCCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTAc                                                                <  1:4708380/93‑1 (MQ=255)
                          aTCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGTCCGTGGTTACa                                                               >  1:2856323/1‑93 (MQ=255)
                          aTCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACa                                                               >  2:755340/1‑93 (MQ=255)
                           tCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAg                                                              <  2:4727329/93‑1 (MQ=255)
                             aaTAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCa                                                            >  2:1168694/1‑93 (MQ=255)
                              aTAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCAc                                                           <  1:4845530/93‑1 (MQ=255)
                               tAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACg                                                          <  1:121403/93‑1 (MQ=255)
                                aGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGa                                                         >  1:2777161/1‑93 (MQ=255)
                                 ggTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGaa                                                        <  1:87903/93‑1 (MQ=255)
                                 ggTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGaa                                                        >  2:6365589/1‑93 (MQ=255)
                                  gTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAg                                                       <  2:2815555/93‑1 (MQ=255)
                                   tCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGt                                                      <  2:4124231/93‑1 (MQ=255)
                                   tCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGt                                                      <  2:1325376/93‑1 (MQ=255)
                                   gcAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGt                                                      <  2:4083012/92‑1 (MQ=255)
                                    cAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGtt                                                     <  2:322441/93‑1 (MQ=255)
                                    cAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGtt                                                     <  2:2455350/93‑1 (MQ=255)
                                    cAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGtt                                                     <  1:187139/93‑1 (MQ=255)
                                    cAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGtt                                                     <  1:641181/93‑1 (MQ=255)
                                     aGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGTTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTg                                                    >  2:5091933/1‑93 (MQ=255)
                                      ggATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTgg                                                   <  2:828346/93‑1 (MQ=255)
                                      ggATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTgg                                                   >  2:3866306/1‑93 (MQ=255)
                                      ggATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTgg                                                   >  2:1400328/1‑93 (MQ=255)
                                      ggATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTgg                                                   <  2:144093/93‑1 (MQ=255)
                                      ggATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTgg                                                   >  2:4759189/1‑93 (MQ=255)
                                      ggATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTgg                                                   <  2:4446384/93‑1 (MQ=255)
                                      ggATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTgg                                                   >  2:362015/1‑93 (MQ=255)
                                      ggATGGTCGTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTgg                                                   >  1:6518447/1‑93 (MQ=255)
                                       gATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTAGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgc                                                  >  2:4695958/1‑93 (MQ=255)
                                         tGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgc                                                >  1:1023794/1‑93 (MQ=255)
                                         tGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgc                                                >  1:143742/1‑93 (MQ=255)
                                         tGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgc                                                <  2:5545229/93‑1 (MQ=255)
                                           gTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgcag                                              <  1:5996078/93‑1 (MQ=255)
                                           gTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgcag                                              >  1:6107367/1‑93 (MQ=255)
                                            tCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgcagc                                             >  1:4479277/1‑93 (MQ=255)
                                             cTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgcagcg                                            <  2:2221059/93‑1 (MQ=255)
                                              tttGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgcagcgc                                           >  1:5545076/1‑93 (MQ=255)
                                                tGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCAt                                         <  2:3184226/93‑1 (MQ=255)
                                                 gCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATc                                        >  1:5688559/1‑93 (MQ=255)
                                                  cGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCa                                       <  1:127868/93‑1 (MQ=255)
                                                  cGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCa                                       <  1:1752197/93‑1 (MQ=255)
                                                   gTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCAc                                      >  1:5210704/1‑93 (MQ=255)
                                                   gTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCAc                                      >  2:6192245/1‑93 (MQ=255)
                                                      gcgcAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCt                                   <  1:5454889/93‑1 (MQ=255)
                                                        gcAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTcc                                 <  1:3868037/93‑1 (MQ=255)
                                                        gcAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTcc                                 <  1:314897/93‑1 (MQ=255)
                                                        gcAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTcc                                 <  2:2919090/93‑1 (MQ=255)
                                                         cAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCt                                >  1:5646203/1‑93 (MQ=255)
                                                          cgTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTg                               <  1:2817341/92‑1 (MQ=255)
                                                           gTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGa                              >  2:2101419/1‑93 (MQ=255)
                                                           gTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGa                              >  2:1899142/1‑93 (MQ=255)
                                                            tCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGaa                             <  2:2332311/93‑1 (MQ=255)
                                                             cAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAg                            >  1:378626/1‑93 (MQ=255)
                                                             cAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAg                            <  2:1159563/93‑1 (MQ=255)
                                                               gCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAgcg                          <  1:4349411/93‑1 (MQ=255)
                                                                  aCTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGcgcac                       >  1:5185954/1‑93 (MQ=255)
                                                                  aCTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGcgcac                       >  2:5810520/1‑93 (MQ=255)
                                                                    tCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGcgcacgc                     <  2:6266905/93‑1 (MQ=255)
                                                                        cAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAAt                 <  1:1070451/93‑1 (MQ=255)
                                                                         aaGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATa                <  1:4130527/93‑1 (MQ=255)
                                                                         aaGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATa                >  2:134483/1‑93 (MQ=255)
                                                                         aaGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATa                <  2:334398/93‑1 (MQ=255)
                                                                          aGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAg               >  2:1756631/1‑93 (MQ=255)
                                                                          aGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAg               >  1:1724728/1‑93 (MQ=255)
                                                                          aGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAg               >  1:6300100/1‑93 (MQ=255)
                                                                           ggATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGc              <  1:3356143/93‑1 (MQ=255)
                                                                            gATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACTAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGcc             >  2:2314184/1‑93 (MQ=255)
                                                                             aTTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCg            >  2:2862065/1‑93 (MQ=255)
                                                                              ttttGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGc           <  2:4640196/93‑1 (MQ=255)
                                                                                ttGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTa         >  2:3452126/1‑93 (MQ=255)
                                                                                 tGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTaa        <  1:4315477/93‑1 (MQ=255)
                                                                                    tttGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAAcc     >  1:2048660/1‑93 (MQ=255)
                                                                                     ttGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCa    <  2:6691497/93‑1 (MQ=255)
                                                                                       gCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCAtt  <  2:6331065/93‑1 (MQ=255)
                                                                                       gCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCAtt  <  1:3637805/93‑1 (MQ=255)
                                                                                          |                                                                                         
TATCGCCAATCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATT  >  NC_000913/2534867‑2535046

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