Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,847,139 A→G G42G (GGT→GGC)  hycC ← membrane‑spanning protein of hydrogenase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,847,1390AG99.4% 531.2 / ‑4.9 166G42G (GGT→GGC) hycCmembrane‑spanning protein of hydrogenase 3
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (64/101);  total (64/102)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.56e-01

Only 100 of 175 total aligned reads displayed.

GCACATCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTACCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACAAACCA  >  NC_000913/2847051‑2847228
                                                                                        |                                                                                         
gCACATCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCaa                                                                                       <  2:2261334/93‑1 (MQ=255)
 cacaTCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAAc                                                                                      <  1:3428001/93‑1 (MQ=255)
 cacaTCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAAc                                                                                      <  1:4463653/93‑1 (MQ=255)
 cacaTCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAAc                                                                                      <  2:2113394/93‑1 (MQ=255)
  acaTCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAAcc                                                                                     <  1:4495469/93‑1 (MQ=255)
      cGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCgcgc                                                                                 >  2:527691/1‑93 (MQ=255)
         aGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGcgccgc                                                                              >  2:618995/1‑93 (MQ=255)
          gCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCgccgcc                                                                             >  2:4067796/1‑93 (MQ=255)
          gCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCgccgcc                                                                             >  2:3137167/1‑93 (MQ=255)
            ttACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCaa                                                                           >  1:2449358/1‑93 (MQ=255)
            ttACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCaa                                                                           >  1:1938514/1‑93 (MQ=255)
            ttACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCaa                                                                           >  2:525360/1‑93 (MQ=255)
             tACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAAt                                                                          <  2:923450/93‑1 (MQ=255)
              aCCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCGAtt                                                                         >  2:814623/1‑93 (MQ=255)
              aCCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAAtt                                                                         <  1:4592572/93‑1 (MQ=255)
               ccAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTc                                                                        >  1:1711123/1‑93 (MQ=255)
                cggCGACAGCGCACCGCTCAGGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTcc                                                                       <  1:842283/91‑1 (MQ=255)
                 aGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCa                                                                      >  2:1665672/1‑93 (MQ=255)
                   cGACAGGGCACCGCGCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGc                                                                    <  1:1097762/93‑1 (MQ=255)
                   cGACAGCGCTCCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGc                                                                    <  2:1424335/93‑1 (MQ=255)
                    gACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTCTACAGAATGCCAACCGCCCCGCCAATTCCAGca                                                                   >  2:3626199/1‑92 (MQ=255)
                     aCAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTa                                                                  <  2:2694232/93‑1 (MQ=255)
                     aCAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTa                                                                  >  2:1163677/1‑93 (MQ=255)
                      cAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTAt                                                                 <  1:4205114/93‑1 (MQ=255)
                      cAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGAGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTAt                                                                 <  2:4296330/93‑1 (MQ=255)
                       aGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTACCAACCGCGCCGCCAATTCCAGCTATc                                                                <  2:2249031/93‑1 (MQ=255)
                         cgcACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCa                                                              <  2:3731127/93‑1 (MQ=255)
                          gcACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAg                                                             >  1:3072108/1‑93 (MQ=255)
                             ccGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCa                                                          >  1:624622/1‑93 (MQ=255)
                                cTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGGCAAATCCAGCTATCCAAACACtt                                                       >  2:1238704/1‑93 (MQ=255)
                                        acccGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCttt                                               <  1:409518/91‑1 (MQ=255)
                                        aaCCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATACAGCCACTTAACGCttt                                               >  1:2095419/1‑93 (MQ=255)
                                         aCCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCtttt                                              >  1:442713/1‑93 (MQ=255)
                                         aCCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCtttt                                              <  1:1063009/93‑1 (MQ=255)
                                          ccGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCttttt                                             <  1:1397334/93‑1 (MQ=255)
                                          ccGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCttttt                                             <  1:2259764/93‑1 (MQ=255)
                                          ccGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCttttt                                             <  1:2918467/93‑1 (MQ=255)
                                           cgcgCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTg                                            <  1:3888115/93‑1 (MQ=255)
                                            gcgcCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTga                                           <  2:4325936/93‑1 (MQ=255)
                                              gcCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaa                                         <  2:3072108/93‑1 (MQ=255)
                                               ccAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaa                                        <  1:3401734/93‑1 (MQ=255)
                                               ccAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaa                                        <  1:588951/93‑1 (MQ=255)
                                                cagtcagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaag                                       <  1:2368879/93‑1 (MQ=255)
                                                cagtcagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaag                                       >  2:476650/1‑93 (MQ=255)
                                                cagtcagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaag                                       >  1:9810/1‑93 (MQ=255)
                                                 agtcagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaaga                                      <  1:1311602/93‑1 (MQ=255)
                                                 agtcagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaaga                                      <  1:1831720/93‑1 (MQ=255)
                                                 agtcagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaaga                                      <  2:234476/93‑1 (MQ=255)
                                                 agtcagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaaga                                      >  1:4368540/1‑93 (MQ=255)
                                                  gtcagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGGCGCCAATTCTAGCTATCCAGCCACTTAACGCTTTTTgaaaagaa                                     >  2:2993553/1‑93 (MQ=255)
                                                  gtcagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaagaa                                     <  1:3652548/93‑1 (MQ=255)
                                                  gtcagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaagaa                                     >  2:3089601/1‑93 (MQ=255)
                                                  gtcagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaagaa                                     <  1:3316547/93‑1 (MQ=255)
                                                   tcagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaagaaa                                    <  2:3509040/93‑1 (MQ=255)
                                                   tcagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaagaaa                                    >  1:3713469/1‑93 (MQ=255)
                                                    cagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaagaaaa                                   <  1:4054660/93‑1 (MQ=255)
                                                     agtACAGTGATGCCCGCGGCTGCCGTATACAGCCTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGaaaagaaaag                                  <  1:1541876/93‑1 (MQ=255)
                                                     agtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGaaaagaaaag                                  <  1:1349906/93‑1 (MQ=255)
                                                       tACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAagag                                <  2:2660671/93‑1 (MQ=255)
                                                       tACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAagag                                <  2:761943/93‑1 (MQ=255)
                                                        tcAGTAATGCCCGCGCCTGCCGTGTCCAGACTGCCACCCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAgaga                               <  2:2861362/92‑1 (MQ=255)
                                                        aCATTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAgaga                               <  1:2406473/93‑1 (MQ=255)
                                                         cAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGaa                              <  1:3806566/93‑1 (MQ=255)
                                                          aGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGaaa                             >  1:3667531/1‑93 (MQ=255)
                                                          aGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGaaa                             <  2:736406/93‑1 (MQ=255)
                                                           gTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAAt                            <  2:1717665/93‑1 (MQ=255)
                                                           gTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAAt                            <  2:1533165/93‑1 (MQ=255)
                                                           gTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAAt                            <  2:2361437/93‑1 (MQ=255)
                                                            tGAAGCCCGCGGCTGCCGTATGCAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCCCTTAACGCTTTTTGAAAAGAAAAGAGAAATg                           <  1:1208728/93‑1 (MQ=255)
                                                            tGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATg                           <  1:4321512/93‑1 (MQ=255)
                                                                  cccgggCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCAATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGaa                     <  2:1485910/89‑1 (MQ=255)
                                                                  ccGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGaa                     >  2:4598354/1‑93 (MQ=255)
                                                                  ccGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGaa                     <  1:4542548/93‑1 (MQ=255)
                                                                  cagcgGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGaa                     <  1:223005/91‑1 (MQ=255)
                                                                   cgcgGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAAc                    <  1:138809/93‑1 (MQ=255)
                                                                   cgcgGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAAc                    <  2:3195055/93‑1 (MQ=255)
                                                                     cgGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAg                  <  1:3451295/93‑1 (MQ=255)
                                                                      ggCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgc                 <  2:4459268/93‑1 (MQ=255)
                                                                      ggCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgc                 >  2:68514/1‑93 (MQ=255)
                                                                       gCTGCCGTATAGAGACTCCAAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgcg                <  2:3988092/93‑1 (MQ=255)
                                                                       gCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgcg                <  1:196598/93‑1 (MQ=255)
                                                                        cTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgcgg               >  2:2107916/1‑93 (MQ=255)
                                                                         tGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTTACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgcggc              >  2:1878211/1‑93 (MQ=255)
                                                                         tGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgcggc              <  2:3708158/93‑1 (MQ=255)
                                                                          gCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgcggcg             >  2:1127851/1‑93 (MQ=255)
                                                                           ccGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgcggcgg            <  1:2862606/93‑1 (MQ=255)
                                                                            cGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgcggcggc           <  1:4623046/93‑1 (MQ=255)
                                                                             gTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgcggcggcg          >  2:4155910/1‑93 (MQ=255)
                                                                              tataCAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGa         >  2:3277520/1‑93 (MQ=255)
                                                                              tataCAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGa         >  1:3536507/1‑93 (MQ=255)
                                                                                taCAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACa       >  2:2409105/1‑93 (MQ=255)
                                                                                  cAGACTGCCAACCGCGCCGCTAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACaaa     <  1:2097443/93‑1 (MQ=255)
                                                                                  cAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACaaa     <  1:705682/93‑1 (MQ=255)
                                                                                  cAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACaaa     >  1:525897/1‑93 (MQ=255)
                                                                                  cAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACaaa     <  1:4223046/93‑1 (MQ=255)
                                                                                    gACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACAAAcc   >  2:3023998/1‑93 (MQ=255)
                                                                                    gACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACAAAcc   <  2:2980066/93‑1 (MQ=255)
                                                                                    gACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACAAAcc   >  1:906228/1‑93 (MQ=255)
                                                                                     aCTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACAAAcca  <  2:3968106/93‑1 (MQ=255)
                                                                                     aCTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACAAAcca  >  2:452923/1‑93 (MQ=255)
                                                                                        |                                                                                         
GCACATCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTACCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACAAACCA  >  NC_000913/2847051‑2847228

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