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.4% 496.0 / ‑5.0 159S101P (TCG→CCG) pdxKpyridoxal/pyridoxine/pyridoxamine kinase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (84/74);  total (85/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.70e-01

Only 100 of 162 total aligned reads displayed.

TCGCCAATCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATT  >  NC_000913/2534869‑2535046
                                                                                        |                                                                                         
tCGCCAATCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGcc                                                                                       >  2:3757617/1‑93 (MQ=255)
tCGCCAATCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGcc                                                                                       <  2:1578391/93‑1 (MQ=255)
  gCCAATCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGt                                                                                     <  2:3742717/93‑1 (MQ=255)
   ccAATCACCGGATCGACCATGATCAATAGGTCAGGATTGTCTTTGCGTAGCGCAGCCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGtt                                                                                    >  2:2672410/1‑93 (MQ=255)
     aaTCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTcc                                                                                  >  2:3201246/1‑93 (MQ=255)
     aaTCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTcc                                                                                  >  2:2379315/1‑93 (MQ=255)
      aTCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTccc                                                                                 >  1:3598407/1‑93 (MQ=255)
      aTCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTccc                                                                                 >  2:26117/1‑93 (MQ=255)
       tCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCa                                                                                <  2:860063/93‑1 (MQ=255)
       tCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCa                                                                                >  1:674525/1‑93 (MQ=255)
         aCCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCTGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCata                                                                              <  2:4549116/93‑1 (MQ=255)
          ccGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCatat                                                                             >  1:1425442/1‑93 (MQ=255)
           cGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCAtata                                                                            >  1:1505854/1‑93 (MQ=255)
           cGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCAtata                                                                            <  1:3337208/93‑1 (MQ=255)
             gATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGc                                                                          >  1:4131541/1‑93 (MQ=255)
              aTCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGcc                                                                         >  1:3496456/1‑93 (MQ=255)
              aTCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGcc                                                                         <  2:3831726/93‑1 (MQ=255)
              aTCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGcc                                                                         <  1:1354064/93‑1 (MQ=255)
                cGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCg                                                                       <  1:4240132/93‑1 (MQ=255)
                   ccATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTgg                                                                    >  2:4139844/1‑93 (MQ=255)
                   ccATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTgg                                                                    <  1:629034/93‑1 (MQ=255)
                    cATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGt                                                                   >  2:1779140/1‑93 (MQ=255)
                    cATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGt                                                                   >  1:137223/1‑93 (MQ=255)
                     aTGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGtt                                                                  >  1:2495467/1‑93 (MQ=255)
                     aTGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGtt                                                                  >  2:3893572/1‑93 (MQ=255)
                      tGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATCGCCCGTGGTTa                                                                 >  2:1262359/1‑93 (MQ=255)
                      tGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTa                                                                 >  1:951541/1‑93 (MQ=255)
                       gATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTAc                                                                <  2:3204081/93‑1 (MQ=255)
                       gATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTAc                                                                >  2:3318664/1‑93 (MQ=255)
                        aTCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACa                                                               <  2:221807/93‑1 (MQ=255)
                        aTCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACa                                                               >  1:2005829/1‑93 (MQ=255)
                         tCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAg                                                              <  2:2695211/93‑1 (MQ=255)
                         tCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAg                                                              >  2:3890020/1‑93 (MQ=255)
                         tCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAg                                                              <  1:3887656/93‑1 (MQ=255)
                          cAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGc                                                             >  1:1131523/1‑93 (MQ=255)
                          cAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGc                                                             <  1:1987860/93‑1 (MQ=255)
                           aaTAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCa                                                            >  1:2504959/1‑93 (MQ=255)
                           aaTAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCa                                                            <  1:4458165/93‑1 (MQ=255)
                            aTAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCAc                                                           <  2:2334672/93‑1 (MQ=255)
                                gTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAg                                                       >  2:3376196/1‑93 (MQ=255)
                                gTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAg                                                       <  2:4066625/93‑1 (MQ=255)
                                gTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAg                                                       <  1:1998212/93‑1 (MQ=255)
                                   aGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTg                                                    >  1:3786183/1‑93 (MQ=255)
                                    ggATGGTCTTTGCGCAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTgg                                                   <  1:3479833/93‑1 (MQ=255)
                                       tGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgc                                                >  2:4223566/1‑93 (MQ=255)
                                       tGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcga                                                >  2:3266562/1‑92 (MQ=255)
                                         gTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgcag                                              >  1:4508532/1‑93 (MQ=255)
                                         gTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgcag                                              <  1:1670702/93‑1 (MQ=255)
                                          tCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgcagc                                             >  2:4354364/1‑93 (MQ=255)
                                          tCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgcagc                                             <  1:1390572/93‑1 (MQ=255)
                                          tCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgcagc                                             >  2:2228318/1‑93 (MQ=255)
                                          tCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgcagc                                             <  2:3123795/93‑1 (MQ=255)
                                             ttGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGgcgcagcgca                                          <  2:2087167/93‑1 (MQ=255)
                                              tGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCAt                                         <  2:1530435/93‑1 (MQ=255)
                                                 gTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCAc                                      >  2:2071229/1‑93 (MQ=255)
                                                  tAGCGCAGTCACCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACg                                     <  1:1147678/93‑1 (MQ=255)
                                                   aGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGc                                    <  1:1871106/93‑1 (MQ=255)
                                                    gcgcAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCt                                   >  1:4002425/1‑93 (MQ=255)
                                                    gcgcAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCt                                   <  2:3618110/93‑1 (MQ=255)
                                                      gcAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTcc                                 <  2:4287027/93‑1 (MQ=255)
                                                      gcAGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTcc                                 <  2:3305820/93‑1 (MQ=255)
                                                        aGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTg                               >  1:3120681/1‑93 (MQ=255)
                                                        aGTCAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTg                               >  2:4255030/1‑93 (MQ=255)
                                                           cAGCCACTCGGCGAGAATTTTGATTTGCGGGGCCGTTCCCATATAGCCCGTAGTTACAGCACGAAGTTGGCGCAGCGCAACACGATCCTGGAg                            <  2:2272937/93‑1 (MQ=255)
                                                           cAGCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAg                            <  1:3893572/93‑1 (MQ=255)
                                                             gCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAgcg                          >  1:103880/1‑93 (MQ=255)
                                                             gCCACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGCTGGCGCAGCGCATCACGCTCCTGAAgcg                          >  2:652483/1‑93 (MQ=255)
                                                               cACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGcgca                        <  1:4065338/93‑1 (MQ=255)
                                                               cACTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGcgca                        <  1:807525/93‑1 (MQ=255)
                                                                aCTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGcgcac                       >  2:1907308/1‑93 (MQ=255)
                                                                aCTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGcgcac                       <  1:3318841/93‑1 (MQ=255)
                                                                aCTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGcgcac                       >  1:3456953/1‑93 (MQ=255)
                                                                aCTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGcgcac                       >  2:102796/1‑93 (MQ=255)
                                                                aCTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGcgcac                       <  1:4440168/93‑1 (MQ=255)
                                                                aCTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGcgcac                       >  1:644910/1‑93 (MQ=255)
                                                                aCTCGGCAAGGATTTTGATTTACGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGcgcac                       >  2:2985068/1‑93 (MQ=255)
                                                                 cTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGcgcacg                      >  1:1023452/1‑93 (MQ=255)
                                                                 cTCGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGcgcacg                      >  1:42026/1‑93 (MQ=255)
                                                                   cGGCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGcgcacgca                    >  2:309195/1‑93 (MQ=255)
                                                                     gCAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCaaa                  <  1:3610955/93‑1 (MQ=255)
                                                                      cAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAAt                 >  2:2913795/1‑93 (MQ=255)
                                                                      cAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAAt                 <  2:1748729/93‑1 (MQ=255)
                                                                      cAAGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAAt                 >  1:3979650/1‑93 (MQ=255)
                                                                       aaGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATa                >  1:818851/1‑93 (MQ=255)
                                                                       aaGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATa                >  1:1585660/1‑93 (MQ=255)
                                                                       aaGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATa                <  1:2902342/93‑1 (MQ=255)
                                                                       aaGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATa                <  1:1287719/93‑1 (MQ=255)
                                                                        aGGATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAg               >  1:2621715/1‑93 (MQ=255)
                                                                         ggATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGc              >  2:1142992/1‑93 (MQ=255)
                                                                         ggATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGc              >  2:3595971/1‑93 (MQ=255)
                                                                         ggATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGc              >  1:4309079/1‑93 (MQ=255)
                                                                         ggATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGAACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGc              >  1:3041011/1‑93 (MQ=255)
                                                                          gATTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGcc             <  1:1929124/93‑1 (MQ=255)
                                                                           aTTTTGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCg            >  1:270144/1‑93 (MQ=255)
                                                                           aTTTTGATTCGCGGTGCCGTTCCCATATAGCCCCTGGTTACAGCACGACGTTGGCGCAGCGCATCACTCTCCTGAAGCGCACGCAAATacgcg            >  1:2151512/1‑89 (MQ=255)
                                                                            ttttGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGc           >  1:4540695/1‑93 (MQ=255)
                                                                            ttttGATTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGc           <  1:2724520/93‑1 (MQ=255)
                                                                                 aTTTGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAAc      >  2:180295/1‑93 (MQ=255)
                                                                                    tGCGGTGCCGTTCCTATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCAt   <  2:278783/93‑1 (MQ=255)
                                                                                    tGCGGTGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCAt   <  2:1270980/93‑1 (MQ=255)
                                                                                        |                                                                                         
TCGCCAATCACCGGATCGACCATGATCAATAGGTCAGGATGGTCTTTGCGTAGCGCAGTCAGCCACTCGGCAAGGATTTTGATTTGCGATGCCGTTCCCATATAGCCCGTGGTTACAGCACGAAGTTGGCGCAGCGCATCACGCTCCTGAAGCGCACGCAAATAGCCGCTAAACCATT  >  NC_000913/2534869‑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: