Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,658,592 C→T 49.5% A175V (GCC→GTC)  ynfE → putative selenate reductase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,658,5920CT49.5% 6.6 / 115.0 89A175V (GCC→GTC) ynfEputative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (21/24);  new base T (18/26);  total (39/50)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

GCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACA  >  NC_000913/1658493‑1658687
                                                                                                   |                                                                                               
gCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGc                                                                                                 <  2:3227956/100‑1 (MQ=255)
  cTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCt                                                                                               >  2:1532917/1‑100 (MQ=255)
  cTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCt                                                                                               >  2:1437958/1‑100 (MQ=255)
   tGAAGAAAACCGTCGCACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTc                                                                                              >  1:1577040/1‑100 (MQ=255)
    gaagacAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGtctc                                                                                              <  1:463612/99‑1 (MQ=255)
    gaagaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTcg                                                                                             >  1:791993/1‑100 (MQ=255)
    gaagaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTcg                                                                                             >  2:183856/1‑100 (MQ=255)
     aagaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTcgg                                                                                            <  2:1494495/100‑1 (MQ=255)
     aagaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTcgg                                                                                            >  2:1994747/1‑100 (MQ=255)
      agaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTcggc                                                                                           <  2:1389502/100‑1 (MQ=255)
       gaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTcggcg                                                                                          <  1:994148/100‑1 (MQ=255)
         aaaCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGt                                                                                        >  1:1642550/1‑100 (MQ=255)
            ccGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCaa                                                                                     >  2:2947649/1‑100 (MQ=255)
             cGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCaaa                                                                                    >  1:705284/1‑100 (MQ=255)
              gTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAAc                                                                                   <  2:2213541/100‑1 (MQ=255)
                cGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATACGACCCGCTCTTCGCCATACGTCTTCGCGGTCAAACGc                                                                                 >  1:407074/1‑100 (MQ=255)
                 gAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGcc                                                                                <  2:821733/100‑1 (MQ=255)
                 gAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGcc                                                                                <  1:62972/100‑1 (MQ=255)
                 gAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGcc                                                                                <  2:3112396/100‑1 (MQ=255)
                   aCAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCtg                                                                              >  2:959211/1‑100 (MQ=255)
                   aCAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtg                                                                              >  1:1891184/1‑100 (MQ=255)
                   aCAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtg                                                                              <  1:1852762/100‑1 (MQ=255)
                     aaTCTGGAAAAGAGGCTGTATCTACTCAGTCCCCCTCGCGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCCACCGCCtgc                                                                             <  1:155213/99‑2 (MQ=255)
                     aaTATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtgat                                                                            >  1:2391523/1‑100 (MQ=255)
                     aaTATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtgat                                                                            <  2:1928020/100‑1 (MQ=255)
                      atatGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtgatg                                                                           >  1:2931644/1‑100 (MQ=255)
                        atGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGaa                                                                         >  1:835645/1‑100 (MQ=255)
                        aggggaaTGAGGCGTTATATCTTCTGTTCTCTGGATGGGTCGTTGGCGGCAATATAACCCGCTCTGCGCCTTCAGTCTCGGCGGTCAAACGCCtg                                                                              <  1:334468/90‑1 (MQ=255)
                         tGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAAc                                                                        >  2:1570904/1‑100 (MQ=255)
                          ggCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACt                                                                       <  2:2590193/100‑1 (MQ=255)
                          ggCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACt                                                                       <  2:2631796/100‑1 (MQ=255)
                           gCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTg                                                                      <  2:1977602/100‑1 (MQ=255)
                            gatctaGGCTGTATATAGTCAGTACTCTTCGTGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGc                                                                     <  2:2829539/95‑1 (MQ=255)
                            cAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGc                                                                     >  1:2141300/1‑100 (MQ=255)
                                gAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCCCACGCCTGATGAACTGCTAcg                                                                 >  2:281732/1‑100 (MQ=255)
                                 aGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGcagatc                                                                >  1:2829516/1‑95 (MQ=255)
                                  ggCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTAcggc                                                               >  1:2151950/1‑100 (MQ=255)
                                   gCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTAcggcg                                                              <  2:1598885/100‑1 (MQ=255)
                                   gCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTAcggcg                                                              >  1:2157519/1‑100 (MQ=255)
                                    cTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTAcggcgg                                                             <  2:2638564/100‑1 (MQ=255)
                                       tatataTTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTc                                                          <  2:2367292/100‑1 (MQ=255)
                                       tatataTTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTc                                                          >  1:967181/1‑100 (MQ=255)
                                       tatataTTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTc                                                          <  2:334732/100‑1 (MQ=255)
                                         tatatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCAc                                                        >  2:2961807/1‑100 (MQ=255)
                                           tatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTc                                                      <  1:2467619/100‑1 (MQ=255)
                                            atTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCa                                                     >  1:512280/1‑100 (MQ=255)
                                               cAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAAcc                                                  >  1:337753/1‑100 (MQ=255)
                                                 gTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAg                                                <  1:657714/100‑1 (MQ=255)
                                                  tACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGt                                               >  1:3148037/1‑100 (MQ=255)
                                                     tctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATg                                            <  2:213552/100‑1 (MQ=255)
                                                      ctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATgg                                           >  1:505271/1‑100 (MQ=255)
                                                      ctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATgg                                           >  2:923450/1‑100 (MQ=255)
                                                      ctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATgg                                           >  2:2240974/1‑100 (MQ=255)
                                                       ttCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGc                                          >  1:3356024/1‑100 (MQ=255)
                                                       ttCGGGGATCGTTAGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGc                                          >  2:2298046/1‑100 (MQ=255)
                                                        tCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGc                                          <  1:1929470/99‑1 (MQ=255)
                                                          ggggATCGTTGGCGGCAATATGACCCGCTCTTCCCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTcc                                       <  2:752714/100‑1 (MQ=255)
                                                             gATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTAc                                    <  2:974602/100‑1 (MQ=255)
                                                              aTCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTAc                                    <  1:1543231/99‑1 (MQ=255)
                                                              aTCGTTAGCGGCAATATCACCCGCTCTTCGCCAGCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATgg                                           <  1:688014/92‑1 (MQ=255)
                                                               tCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAg                                  >  2:2621673/1‑100 (MQ=255)
                                                               tCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAg                                  <  2:593240/100‑1 (MQ=255)
                                                                cGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGc                                 <  2:951023/100‑1 (MQ=255)
                                                                cGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGc                                 <  2:3206955/100‑1 (MQ=255)
                                                                 gTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCa                                >  2:722868/1‑100 (MQ=255)
                                                                 gTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCa                                >  1:1412553/1‑100 (MQ=255)
                                                                  ttGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCAc                               <  2:967184/100‑1 (MQ=255)
                                                                     gcggcAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCCGTTCACTCAACCAGTATGGCTCCTACAGCACTgc                            <  2:387924/100‑1 (MQ=255)
                                                                     gcggcAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTgc                            <  1:2767710/100‑1 (MQ=255)
                                                                      cggcAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTgcg                           <  2:2307316/100‑1 (MQ=255)
                                                                      cggcAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTgcg                           >  1:2094569/1‑100 (MQ=255)
                                                                      cggcAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTgcg                           >  1:2259237/1‑100 (MQ=255)
                                                                      cggcAATACCACCCGCTCTTCGCCACCCGTCTCGGCGGTCAACCGCCTGATGAACTGCTCCCGCGGTTCACTCAACCAGTATGGCTCCTCCAGCACTg                             <  1:2703966/98‑1 (MQ=255)
                                                                          aaTATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCaaa                       <  1:3009626/100‑1 (MQ=255)
                                                                          aaTATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCaaa                       <  2:538188/100‑1 (MQ=255)
                                                                              tGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTc                   >  1:1557125/1‑100 (MQ=255)
                                                                               gACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTcc                  <  2:2847284/100‑1 (MQ=255)
                                                                               gACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGACCTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTcc                  <  1:2068603/100‑1 (MQ=255)
                                                                                 cccGCTCTTCGCCATCAGTCTCGGCGGTCACACGCCTGTTGGACTGATCCGGCGGTTCACTCTACTAGTATGGCTCCTACAGCACTGCGCAAATTTCCtg                <  1:281731/100‑1 (MQ=255)
                                                                                  ccGCTCTTCGCCATCAGTCTCTGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAACTTTCCtgt               <  1:1437952/100‑1 (MQ=255)
                                                                                  ccGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCtgt               <  2:316513/100‑1 (MQ=255)
                                                                                  ccGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCtgt               <  2:2835013/100‑1 (MQ=255)
                                                                                   cGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCtgtg              <  2:242811/100‑1 (MQ=255)
                                                                                     ctctTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGcc            <  2:512282/100‑1 (MQ=255)
                                                                                      tctTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCa           <  2:857017/100‑1 (MQ=255)
                                                                                         tCGCCATCAGCCTCGGCGGTCAAACGTCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGc        >  2:2595754/1‑100 (MQ=255)
                                                                                            ccATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCt     >  2:2768579/1‑100 (MQ=255)
                                                                                              aTCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTac   <  2:803005/100‑1 (MQ=255)
                                                                                              aTCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTac   <  1:1393318/100‑1 (MQ=255)
                                                                                               tCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTaca  <  2:2280244/100‑1 (MQ=255)
                                                                                               tCAGCCACTGCGGTCAAACGCCTGATGAACTGCTTCGGCGGGTCACTCAACCAGTATGGCTCCTACAGCACTGAGCAAATTTCCTGTGCCATGCCCTaa   <  2:1577047/99‑2 (MQ=255)
                                                                                                   |                                                                                               
GCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACA  >  NC_000913/1658493‑1658687

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