Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,551,614 T→C L382P (CTG→CCG)  uxuB → D‑mannonate oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,551,6140TC99.7% 935.1 / ‑5.7 302L382P (CTG→CCG) uxuBD‑mannonate oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (159/141);  total (160/142)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 329 total aligned reads displayed.

CGCTGCTGATCGAGCGTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCTGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCGTGGCTGGCTGGATGCGTTACACCCAGG  >  NC_000913/4551524‑4551703
                                                                                          |                                                                                         
 gctgctGATCGAGCGTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCCGAAGTTACCGCAGCGTccgct                                                                                        <  1:1368718/93‑1 (MQ=255)
 gctgctGATCGAGCGTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTccgct                                                                                        >  1:281425/1‑93 (MQ=255)
  ctgctgATCGAGCGTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTccgctg                                                                                       <  2:5562116/93‑1 (MQ=255)
   tgctgATCGAGCGTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTgg                                                                                      >  1:4004557/1‑93 (MQ=255)
    gctgATCGAGCGTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGa                                                                                     <  2:5900837/93‑1 (MQ=255)
      tgATCGAGCGTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGAcc                                                                                   >  2:1279853/1‑93 (MQ=255)
       gATCGAGCGTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGAccc                                                                                  <  2:2633696/93‑1 (MQ=255)
       gATCGAGCGTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGAccc                                                                                  <  1:4477704/93‑1 (MQ=255)
       gATCGAGCGTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGAccc                                                                                  >  2:2307760/1‑93 (MQ=255)
       gATCGAGCGTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGAccc                                                                                  <  2:6451070/93‑1 (MQ=255)
        aTCGAGCGTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCg                                                                                 >  2:1795475/1‑93 (MQ=255)
         tCGAGCGTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCgg                                                                                <  1:2293881/93‑1 (MQ=255)
         tCGAGCGTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCgg                                                                                <  2:2219422/93‑1 (MQ=255)
         tCGAGCGTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCgg                                                                                >  1:7403116/1‑93 (MQ=255)
         tCGAGCATTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCgg                                                                                <  2:4728473/93‑1 (MQ=255)
          cGAGCGTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCggg                                                                               >  1:3058795/1‑92 (MQ=255)
          cGAGCGTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGt                                                                               <  2:2137121/93‑1 (MQ=255)
          cGAGCGTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGt                                                                               <  2:499394/93‑1 (MQ=255)
           gAGCGTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTTCGATGGACGGCAGACAGAAGTTACCGCAGCGTCCGCTGGACCCGGTg                                                                              >  2:3027778/1‑93 (MQ=255)
           gAGCGTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTg                                                                              <  1:2147252/93‑1 (MQ=255)
           gAGCGTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTg                                                                              >  2:3718839/1‑93 (MQ=255)
           gAGCGTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTg                                                                              <  1:6861835/93‑1 (MQ=255)
             gCGTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCg                                                                            >  2:93710/1‑93 (MQ=255)
             gCGTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCg                                                                            <  1:3189285/93‑1 (MQ=255)
               gTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTc                                                                          <  2:910043/93‑1 (MQ=255)
               gTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTc                                                                          >  1:6118825/1‑93 (MQ=255)
                tttCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGTGTCCGCTGGACCCGGTGCGTCt                                                                         >  1:211866/1‑93 (MQ=255)
                 ttCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGGGCGTCtc                                                                        >  1:55028/1‑92 (MQ=255)
                    aGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCAc                                                                     <  2:6804841/93‑1 (MQ=255)
                     gCAAACCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCAcc                                                                    >  1:7836594/1‑93 (MQ=255)
                       aaCCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTg                                                                  >  2:7113312/1‑93 (MQ=255)
                       aaCCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTg                                                                  >  2:1828330/1‑93 (MQ=255)
                          ccGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCaa                                                               >  2:7814028/1‑93 (MQ=255)
                          ccGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCaa                                                               >  2:2050842/1‑93 (MQ=255)
                            gTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCaaaa                                                             <  2:733799/93‑1 (MQ=255)
                             tctcTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAAc                                                            <  1:3553445/93‑1 (MQ=255)
                             tctcTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAAc                                                            <  2:4274717/93‑1 (MQ=255)
                             tctcTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAAc                                                            >  1:7902358/1‑93 (MQ=255)
                             gctctGCGTCGCCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAAc                                                            <  2:4415164/92‑1 (MQ=255)
                               tctGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAAcgg                                                          >  2:6127321/1‑93 (MQ=255)
                               tctGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAAcgg                                                          >  1:1047261/1‑93 (MQ=255)
                               tctGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAAcgg                                                          >  1:4657523/1‑93 (MQ=255)
                               tctGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAAcgg                                                          >  2:4875/1‑93 (MQ=255)
                                ctGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAAcggc                                                         <  2:3850385/93‑1 (MQ=255)
                                ctGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAAcggc                                                         <  1:3220182/93‑1 (MQ=255)
                                ctGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAAcggc                                                         >  1:2445243/1‑93 (MQ=255)
                                 tGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAAcggcg                                                        >  1:3818438/1‑93 (MQ=255)
                                 tGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAAcggcg                                                        >  2:2760075/1‑93 (MQ=255)
                                   cGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACggcggc                                                      <  1:3701618/93‑1 (MQ=255)
                                    gTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCa                                                     >  2:6152477/1‑93 (MQ=255)
                                      cACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGc                                                   <  2:2696266/93‑1 (MQ=255)
                                      cACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGc                                                   >  2:3509320/1‑93 (MQ=255)
                                      cACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGc                                                   >  1:1657102/1‑93 (MQ=255)
                                      cACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGc                                                   >  1:7060381/1‑93 (MQ=255)
                                       cccGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCt                                                  <  1:848151/92‑1 (MQ=255)
                                        ccGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTg                                                 >  1:6245606/1‑93 (MQ=255)
                                         cGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTgg                                                >  2:7453073/1‑93 (MQ=255)
                                           tACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCg                                              >  2:4214902/1‑93 (MQ=255)
                                            aCCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGt                                             <  1:664264/93‑1 (MQ=255)
                                            aCCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGt                                             <  1:4921022/93‑1 (MQ=255)
                                              cTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTcc                                           >  1:3910324/1‑92 (MQ=255)
                                                ggCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCAcc                                         <  2:2447081/93‑1 (MQ=255)
                                                 gCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGGGTCACCt                                        <  2:6667868/93‑1 (MQ=255)
                                                  cAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGTGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTg                                       <  2:6373095/93‑1 (MQ=255)
                                                  cAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTg                                       >  1:1013653/1‑93 (MQ=255)
                                                  cAGATGGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTg                                       <  1:4083187/93‑1 (MQ=255)
                                                   aGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTgg                                      <  2:5596664/93‑1 (MQ=255)
                                                    gATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGgc                                     <  2:6847252/93‑1 (MQ=255)
                                                       tGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCt                                  <  2:7087043/93‑1 (MQ=255)
                                                       tGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCt                                  <  1:5955118/93‑1 (MQ=255)
                                                         cGATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTgg                                <  2:3117194/93‑1 (MQ=255)
                                                          gATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTggg                               >  1:5792489/1‑93 (MQ=255)
                                                          gATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTggg                               >  2:1027977/1‑93 (MQ=255)
                                                          gATGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTggg                               >  1:7308071/1‑93 (MQ=255)
                                                           aTGGACGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGc                              >  1:5330745/1‑93 (MQ=255)
                                                            tGGATGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCg                             >  1:7741024/1‑93 (MQ=255)
                                                               aCGGCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCgggg                          >  1:1887266/1‑93 (MQ=255)
                                                                 ggCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCGtggct                        <  1:5677227/93‑1 (MQ=255)
                                                                  gCAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCGtggctg                       <  1:6744215/93‑1 (MQ=255)
                                                                   cAGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCGtggctgg                      <  1:5010632/93‑1 (MQ=255)
                                                                    aGCCAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCGtggctggc                     <  1:432770/93‑1 (MQ=255)
                                                                      ccAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCGTGgctggctg                   <  1:333300/93‑1 (MQ=255)
                                                                       cAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCGTGGctggctgg                  >  1:1336867/1‑93 (MQ=255)
                                                                       cAGAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCGTGGctggctgg                  >  1:5090067/1‑93 (MQ=255)
                                                                         gAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCGTGGCTGGCTGGAt                >  2:5741822/1‑93 (MQ=255)
                                                                         gAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCGTGGCTGGCTGGAt                >  1:2752687/1‑93 (MQ=255)
                                                                         gAAGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCGTGGCTGGCTGGAt                >  2:1393149/1‑93 (MQ=255)
                                                                          aaGTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCGTGGCTGGCTGGATg               <  2:3908535/93‑1 (MQ=255)
                                                                            gTTACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCGTGGCTGGCTGGATGCg             >  2:3582936/1‑93 (MQ=255)
                                                                             ttACCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCGTGGCTGGCTGGGTGCGt            >  2:327979/1‑93 (MQ=255)
                                                                               tccGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCGTGGCTGGCTGGATGCGTTa          <  1:3431181/92‑1 (MQ=255)
                                                                               aCCGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCGTGGCTGGCTGGATGCGTTa          <  2:5082378/93‑1 (MQ=255)
                                                                                ccGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCGTGGCTGGCTGGATGCGTTac         <  1:6468702/93‑1 (MQ=255)
                                                                                ccGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCGTGGCTGGCTGGATGCGTTac         >  1:3892468/1‑93 (MQ=255)
                                                                                 cGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCGTGGCTGGCTGGATGCGTTaca        <  1:1469341/93‑1 (MQ=255)
                                                                                 cGCAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCGTGGCTGGCTGGATGCGTTaca        >  2:1710600/1‑93 (MQ=255)
                                                                                   cAGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCGTGGCTGGCTGGATGCGTTACAcc      <  2:5983957/93‑1 (MQ=255)
                                                                                    aGCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCGTGGCTGGCTGGATGGGTTACAccc     <  2:5241003/93‑1 (MQ=255)
                                                                                     gCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCGTGGCTGGCTGGATGCGTTACACCCa    >  2:4593415/1‑93 (MQ=255)
                                                                                     gCGTCCGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCGTGGCTGGCTGGATGCGTTACACCCa    >  2:7695715/1‑93 (MQ=255)
                                                                                          |                                                                                         
CGCTGCTGATCGAGCGTTTCAGCAACCCGTCTCTGCGTCACCGTACCTGGCAGATTGCGATGGACGGCAGCCAGAAGTTACCGCAGCGTCTGCTGGACCCGGTGCGTCTGCACCTGCAAAACGGCGGCAGCTGGCGTCACCTGGCGCTGGGCGTGGCTGGCTGGATGCGTTACACCCAGG  >  NC_000913/4551524‑4551703

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