Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,299,266 T→C N503S (AAC→AGC)  napA ← putative nitrate reductase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,299,2660TC99.6% 723.2 / NA 223N503S (AAC→AGC) napAputative nitrate reductase 3
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (87/135);  minor base G (1/0);  total (88/135)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 234 total aligned reads displayed.

CCCTCTTTCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTC  >  NC_000913/2299176‑2299355
                                                                                          |                                                                                         
cccTCTTTCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTg                                                                                         >  2:1768773/1‑93 (MQ=255)
 ccTCTTTCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTgc                                                                                        >  2:3836682/1‑93 (MQ=255)
  ctctTTCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTgcg                                                                                       <  1:1156506/93‑1 (MQ=255)
   tctTTCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTgcgc                                                                                      <  2:1166686/93‑1 (MQ=255)
    ctTTCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGcgcg                                                                                     >  2:3029317/1‑93 (MQ=255)
    ctTTCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGcgcg                                                                                     >  2:1545363/1‑93 (MQ=255)
     tgtctcTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCgg                                                                                    <  2:4924618/91‑1 (MQ=255)
      ttCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGa                                                                                   <  1:2571099/93‑1 (MQ=255)
       tctcTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGAt                                                                                  <  1:518023/93‑1 (MQ=255)
        ctctACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATc                                                                                 <  2:2391872/93‑1 (MQ=255)
         tctACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATcg                                                                                >  2:5986078/1‑93 (MQ=255)
          ctACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATcgc                                                                               >  2:6606112/1‑93 (MQ=255)
             cccACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGcc                                                                            <  1:2106671/93‑1 (MQ=255)
             cccACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGcc                                                                            <  1:3862690/93‑1 (MQ=255)
             cccACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGcc                                                                            >  2:1444643/1‑93 (MQ=255)
                acaTTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGc                                                                         >  2:2450558/1‑93 (MQ=255)
                acaTTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGc                                                                         >  2:3473988/1‑93 (MQ=255)
                acaTTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGc                                                                         >  1:2595447/1‑93 (MQ=255)
                 caTTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGcc                                                                        >  2:5394437/1‑93 (MQ=255)
                 caTTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGcc                                                                        <  2:5183641/93‑1 (MQ=255)
                  aTTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGccc                                                                       <  2:1381262/93‑1 (MQ=255)
                   ttGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCg                                                                      >  1:1253718/1‑93 (MQ=255)
                   ttGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCg                                                                      >  1:4543314/1‑93 (MQ=255)
                    tGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCgg                                                                     <  1:3683096/93‑1 (MQ=255)
                      cgctCGGCAAGATCAATTGTGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCa                                                                   <  1:4154769/90‑1 (MQ=255)
                        gTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATa                                                                 <  1:1349740/93‑1 (MQ=255)
                             cAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGctc                                                            <  1:285552/93‑1 (MQ=255)
                             cAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGctc                                                            >  2:2165079/1‑93 (MQ=255)
                              aGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGctct                                                           <  1:72281/93‑1 (MQ=255)
                              aGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGctct                                                           >  1:4441686/1‑93 (MQ=255)
                               ggATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCtt                                                          >  2:4331963/1‑93 (MQ=255)
                               ggATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCtt                                                          <  1:4323879/93‑1 (MQ=255)
                                gATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTc                                                         <  2:2219930/93‑1 (MQ=255)
                                gATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTc                                                         <  1:5833819/93‑1 (MQ=255)
                                  tCAAGTCGGCGGCCAGCGCTCTGACTGTCGGATCCGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGTCATACGCTCTTCAt                                                       <  1:2641941/93‑1 (MQ=255)
                                  tCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCAt                                                       >  1:1388479/1‑93 (MQ=255)
                                   cAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCAtt                                                      <  2:3924269/93‑1 (MQ=255)
                                    aaGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTg                                                     >  1:479411/1‑93 (MQ=255)
                                    aaGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTg                                                     >  1:1695818/1‑93 (MQ=255)
                                     aGTCGGCGGCCAGCGCACTGACTGTCGGATCCGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGa                                                    <  1:6376440/93‑1 (MQ=255)
                                     aGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGa                                                    <  1:5812415/93‑1 (MQ=255)
                                     aGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGa                                                    <  1:2365294/93‑1 (MQ=255)
                                      gTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATGGCGCCCGCCCCGCGTAGCGTCTTCCTTGAt                                                   >  1:69596/1‑93 (MQ=255)
                                      gTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGAt                                                   <  1:4812464/93‑1 (MQ=255)
                                       tCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATg                                                  >  1:5044582/1‑93 (MQ=255)
                                        cggcggCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGt                                                 <  2:3132034/93‑1 (MQ=255)
                                        cggcggCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGt                                                 <  2:6275706/93‑1 (MQ=255)
                                         ggcggcCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGtt                                                <  1:1813900/93‑1 (MQ=255)
                                          gcggcCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTc                                               <  2:5826095/93‑1 (MQ=255)
                                          gcggcCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTc                                               <  1:6407744/93‑1 (MQ=255)
                                          gcggcCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTc                                               >  2:436507/1‑93 (MQ=255)
                                          gcggcCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTc                                               <  1:267050/93‑1 (MQ=255)
                                          gcggcCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTc                                               <  2:3599768/93‑1 (MQ=255)
                                           cggcCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTcg                                              <  2:6108013/93‑1 (MQ=255)
                                              ccAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTcggcc                                           <  1:5074954/93‑1 (MQ=255)
                                              ccAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTcggcc                                           >  2:6397599/1‑93 (MQ=255)
                                               cAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTcggccc                                          <  1:4051386/93‑1 (MQ=255)
                                               cAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTcggccc                                          <  1:5008055/93‑1 (MQ=255)
                                                aGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTcggcccg                                         >  1:5637833/1‑93 (MQ=255)
                                                aGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGACAGCCCGGCATACGCTCTTCATTGATGTTcggcccg                                         >  2:5540569/1‑93 (MQ=255)
                                                 gcgcACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCTCCAGCCCGGCATACGCTCTTCATTGATGTTcggcccgg                                        <  1:3409784/93‑1 (MQ=255)
                                                  ggcACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTcggcccggc                                       <  1:1479871/92‑1 (MQ=255)
                                                   gcACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTcggcccgggc                                      >  1:601260/1‑91 (MQ=255)
                                                   gcACTGACTGTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTcggcccggcc                                      >  1:5414557/1‑93 (MQ=255)
                                                   gcACTGACTGTCGGATACGGATCGGAGACGATGAAGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTcggcccggcc                                      <  1:2065728/93‑1 (MQ=255)
                                                     actgactgTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTg                                    >  1:896507/1‑93 (MQ=255)
                                                      ctgactgTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGc                                   <  1:5717936/93‑1 (MQ=255)
                                                      ctgactgTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGc                                   >  2:3362789/1‑93 (MQ=255)
                                                         actgTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCAtg                                >  1:3855316/1‑93 (MQ=255)
                                                         actgTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCAtg                                >  2:6682593/1‑93 (MQ=255)
                                                          ctgTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCAtgt                               <  2:1061018/93‑1 (MQ=255)
                                                           tgTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCAtgtt                              <  1:1053656/93‑1 (MQ=255)
                                                           tgTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCAtgtt                              <  1:5978198/93‑1 (MQ=255)
                                                            gTCGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCAtgttg                             >  1:3036027/1‑93 (MQ=255)
                                                              cGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCAtgttgtc                           >  1:792620/1‑92 (MQ=255)
                                                              cGGATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATgttgtt                           >  2:2658506/1‑93 (MQ=255)
                                                               ggATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGttgttg                          >  2:2028851/1‑93 (MQ=255)
                                                               ggATACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGttgttg                          <  2:6170652/93‑1 (MQ=255)
                                                                 aTACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTa                       >  1:2193130/1‑93 (MQ=255)
                                                                  tACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTa                       <  2:6492517/93‑1 (MQ=255)
                                                                  tACGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTa                       >  1:3319904/1‑93 (MQ=255)
                                                                    cGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTaca                     <  2:2678817/93‑1 (MQ=255)
                                                                    cGGATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTaca                     <  2:6431557/93‑1 (MQ=255)
                                                                      gATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTAcaca                   <  2:78860/93‑1 (MQ=255)
                                                                      gATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTAcaca                   >  2:68709/1‑93 (MQ=255)
                                                                      gATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTAcaca                   <  2:110284/93‑1 (MQ=255)
                                                                      gATCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTAcaca                   <  1:4881896/93‑1 (MQ=255)
                                                                       aTCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACAt                  <  2:3246533/93‑1 (MQ=255)
                                                                       aTCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACAt                  <  1:2749519/93‑1 (MQ=255)
                                                                        tCGGAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATg                 >  2:1398977/1‑93 (MQ=255)
                                                                          ggAGACGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGt               <  2:4229895/93‑1 (MQ=255)
                                                                           gagaCGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTc              >  1:4231031/1‑93 (MQ=255)
                                                                           gagaCGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTc              >  2:5634804/1‑93 (MQ=255)
                                                                           gagaCGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTc              <  1:354580/93‑1 (MQ=255)
                                                                           gagaCGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTc              >  1:6217050/1‑93 (MQ=255)
                                                                           gagaCGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTc              <  1:5703346/93‑1 (MQ=255)
                                                                            agaCGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTcc             >  1:5240465/1‑93 (MQ=255)
                                                                             gaCGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCa            >  1:1828017/1‑93 (MQ=255)
                                                                             gaCGATGATGAAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCa            <  1:2938508/93‑1 (MQ=255)
                                                                                     tgaAGCTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACAt    >  2:1750236/1‑93 (MQ=255)
                                                                                          |                                                                                         
CCCTCTTTCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTC  >  NC_000913/2299176‑2299355

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