Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,843,138 G→A P453S (CCG→TCG)  hycE ← large subunit of hydrogenase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,843,1380GA99.5% 730.1 / ‑6.0 202P453S (CCG→TCG) hycElarge subunit of hydrogenase 3
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (79/122);  minor base T (0/1);  total (79/123)
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.99e-01

Only 100 of 211 total aligned reads displayed.

GGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCGCAGCCC  >  NC_000913/2843048‑2843226
                                                                                          |                                                                                        
 gATATCATCGCCGCGCGGCGCTGCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAg                                                                                       <  2:2649283/93‑1 (MQ=255)
   tatCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGt                                                                                     <  2:3158452/93‑1 (MQ=255)
   tatCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGt                                                                                     <  1:3579748/93‑1 (MQ=255)
   tatCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGt                                                                                     >  1:1733878/1‑93 (MQ=255)
    atcatcGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGtt                                                                                    <  2:2345006/93‑1 (MQ=255)
     tcatcGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTa                                                                                   >  1:1096212/1‑93 (MQ=255)
     tcatcGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTa                                                                                   >  1:1409198/1‑93 (MQ=255)
     tcatcGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTa                                                                                   <  2:1977316/93‑1 (MQ=255)
       atcGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATc                                                                                 <  1:1061432/93‑1 (MQ=255)
       atcGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATc                                                                                 <  2:2113824/93‑1 (MQ=255)
       atcGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATc                                                                                 <  2:3519550/93‑1 (MQ=255)
        tcGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATcc                                                                                <  2:1882255/93‑1 (MQ=255)
        tcGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGATATcc                                                                                >  2:3702759/1‑93 (MQ=255)
          cccgcGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAg                                                                              <  2:3374636/92‑1 (MQ=255)
            cgcgcgGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGAc                                                                            >  2:3806673/1‑93 (MQ=255)
                 ggcgCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAg                                                                       >  1:911852/1‑93 (MQ=255)
                  gcgcTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGt                                                                      <  1:2447781/93‑1 (MQ=255)
                    gcTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCg                                                                    <  1:4288419/93‑1 (MQ=255)
                    gcTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCg                                                                    <  1:3557381/93‑1 (MQ=255)
                    gcTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCg                                                                    >  1:3377494/1‑93 (MQ=255)
                    gcTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCg                                                                    <  1:2512824/93‑1 (MQ=255)
                    gcTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCg                                                                    >  2:4284691/1‑93 (MQ=255)
                     cTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGa                                                                   >  1:3445863/1‑93 (MQ=255)
                      ttCGGCAAAGCCCAGCGCAGAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGAt                                                                  >  2:4584788/1‑93 (MQ=255)
                       tCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCTACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTGGAta                                                                 >  2:3102169/1‑92 (MQ=255)
                        cGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCa                                                                <  2:4375595/93‑1 (MQ=255)
                        cGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCa                                                                <  1:1477598/93‑1 (MQ=255)
                        cGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCa                                                                <  1:2002292/93‑1 (MQ=255)
                        cGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCa                                                                <  2:1388127/93‑1 (MQ=255)
                         ggCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCAt                                                               <  1:2200272/93‑1 (MQ=255)
                          gCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGGGGGCCACCCGACAGGTTATCCAGACCGGAGTCGATCAtt                                                              >  2:2313271/1‑92 (MQ=255)
                            aaaGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGtt                                                            >  2:3749686/1‑93 (MQ=255)
                            aaaGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGtt                                                            >  2:2971414/1‑93 (MQ=255)
                            aaaGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGtc                                                            >  1:1501836/1‑92 (MQ=255)
                             aaGCCCATCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTc                                                           >  2:4454756/1‑93 (MQ=255)
                             aaGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTc                                                           <  1:525760/93‑1 (MQ=255)
                             aaGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTc                                                           >  2:1553415/1‑93 (MQ=255)
                              aGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCa                                                          <  2:1703706/93‑1 (MQ=255)
                              aGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCa                                                          <  2:1225556/93‑1 (MQ=255)
                              aGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCa                                                          <  2:4078694/93‑1 (MQ=255)
                                cccAGCGCAAAGCGGTGCGGAATGTAGGTACAGCCCTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAgc                                                        <  1:2278236/93‑1 (MQ=255)
                                cccAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAgc                                                        <  2:689696/93‑1 (MQ=255)
                                cccAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAgc                                                        <  2:1067481/93‑1 (MQ=255)
                                  cAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAgcgc                                                      <  1:3212232/93‑1 (MQ=255)
                                    gcgcAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCgg                                                    >  1:4507099/1‑93 (MQ=255)
                                     cgcAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGt                                                   >  1:2101179/1‑93 (MQ=255)
                                     cgcAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGt                                                   <  2:332832/93‑1 (MQ=255)
                                      gcAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGta                                                  >  2:1323202/1‑93 (MQ=255)
                                      gcAAAGCGGGGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCAATCATGTTCAGCGCGGta                                                  <  2:1272104/93‑1 (MQ=255)
                                        aacgCGGTGCGTAATGTAAGTAAAGCCTTCCACGATCAGTGGGCCACCTGTCAGGTTATCCAGTCCGTAGTTGATCATGTTCAGGGCGGtata                                                <  2:4152774/90‑1 (MQ=255)
                                        aaaGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGAGAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGtata                                                <  2:3669855/93‑1 (MQ=255)
                                        aaaGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGtata                                                <  1:4188505/93‑1 (MQ=255)
                                         aaGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAg                                               >  1:491335/1‑93 (MQ=255)
                                           gCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGAc                                             >  1:1322793/1‑93 (MQ=255)
                                           gCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGAc                                             <  2:2450560/93‑1 (MQ=255)
                                           gCGGTGCGGAAGGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGAc                                             <  2:2796308/93‑1 (MQ=255)
                                            cGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACt                                            <  1:4048907/93‑1 (MQ=255)
                                            cGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGAAGTGGATCATGTTCAGCGCGGTATAGACt                                            >  2:3237143/1‑93 (MQ=255)
                                              gTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTc                                          <  2:1178481/93‑1 (MQ=255)
                                              gTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTc                                          <  2:1024296/93‑1 (MQ=255)
                                              gTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTc                                          <  1:4346178/93‑1 (MQ=255)
                                                   gAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGa                                     <  2:1895718/93‑1 (MQ=255)
                                                   gAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGACCATGTTCAGCGCGGTATAGACTTCGTTGa                                     <  1:272234/93‑1 (MQ=255)
                                                    aaTGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGAt                                    <  1:4605181/93‑1 (MQ=255)
                                                    aaTGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGAt                                    <  2:3034954/93‑1 (MQ=255)
                                                    aaTGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGAt                                    <  2:3921284/93‑1 (MQ=255)
                                                     aTGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATa                                   >  1:1615825/1‑93 (MQ=255)
                                                     aTGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATa                                   <  1:2956664/93‑1 (MQ=255)
                                                      tGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATAc                                  <  2:342765/93‑1 (MQ=255)
                                                      tGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATAc                                  >  2:854452/1‑93 (MQ=255)
                                                      tGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATAc                                  >  1:1070863/1‑93 (MQ=255)
                                                       gTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACg                                 <  2:3573605/93‑1 (MQ=255)
                                                         aGGTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCa                               <  1:1227569/93‑1 (MQ=255)
                                                           gTAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACt                             <  2:3649216/93‑1 (MQ=255)
                                                            tAAAGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACtt                            >  2:411658/1‑93 (MQ=255)
                                                             aaaGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACttt                           <  1:1145477/93‑1 (MQ=255)
                                                             aaaGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACttt                           >  1:3868383/1‑93 (MQ=255)
                                                              aaGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTc                          >  1:2749128/1‑93 (MQ=255)
                                                              aaGCCTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTc                          >  2:2846546/1‑93 (MQ=255)
                                                                  cTTCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGAc                      <  1:1585322/93‑1 (MQ=255)
                                                                  cTTCCACCATCAGTGGGCCACCCGACAGGATATCCAGAACGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGAc                      >  1:467161/1‑93 (MQ=255)
                                                                    tCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACgg                    >  2:2236612/1‑93 (MQ=255)
                                                                    tCCACCATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACgg                    >  1:1495109/1‑93 (MQ=255)
                                                                    tCCACAATCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGCCTTCGTTGATACGCACTTTCAGACgg                    <  2:1983932/93‑1 (MQ=255)
                                                                       accaTCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGaa                 >  1:3583229/1‑93 (MQ=255)
                                                                        ccaTCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGaaa                <  1:3804051/93‑1 (MQ=255)
                                                                        ccaTCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGaaa                <  2:2366658/93‑1 (MQ=255)
                                                                          aTCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATa              <  2:1354968/93‑1 (MQ=255)
                                                                          aTCAGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATa              >  2:317730/1‑93 (MQ=255)
                                                                          aTCAGTGGGACACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATa              >  2:1002113/1‑93 (MQ=255)
                                                                            cAGTGGGCCACCCGTCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAAc            <  1:650778/93‑1 (MQ=255)
                                                                             aGTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGTGAAATAAct           >  1:728503/1‑92 (MQ=255)
                                                                              gTGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGt          <  1:2610103/93‑1 (MQ=255)
                                                                               tgtgCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTc         <  2:775467/90‑1 (MQ=255)
                                                                               tGGGCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTc         <  2:4199066/93‑1 (MQ=255)
                                                                                gggCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCg        >  1:3411644/1‑93 (MQ=255)
                                                                                gggCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCg        >  1:2378712/1‑93 (MQ=255)
                                                                                 ggCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCGc       <  2:1704788/93‑1 (MQ=255)
                                                                                  gCCACCCGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCGCa      <  2:3853454/93‑1 (MQ=255)
                                                                                      cccGACAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCGCAGccc  <  1:4600935/93‑1 (MQ=255)
                                                                                          |                                                                                        
GGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCGCAGCCC  >  NC_000913/2843048‑2843226

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