Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,913,862 T→C D381G (GAT→GGT)  atpD ← membrane‑bound ATP synthase, F1 sector, beta‑subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,913,8620TC99.7% 909.9 / ‑4.8 286D381G (GAT→GGT) atpDmembrane‑bound ATP synthase, F1 sector, beta‑subunit
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (131/154);  total (132/154)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.79e-01

Only 100 of 295 total aligned reads displayed.

AATACTTCTGCCACGAAGAACGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCATCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCACGCGCGGTGTCGTAGTGTTCCTGA  >  NC_000913/3913772‑3913951
                                                                                          |                                                                                         
 aTACTTCTGCCACGAAGAACGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCa                                                                                        <  2:577760/93‑1 (MQ=255)
 aTACTTCTGCCACGAAGAACGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCa                                                                                        >  2:3753006/1‑93 (MQ=255)
    gttCTGCCACGAAGAACGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATAc                                                                                     <  1:2046335/92‑1 (MQ=255)
     ttCTGCCACGAAGAACGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATAcc                                                                                    <  1:403009/93‑1 (MQ=255)
     ttCTGCCACGAAGAACGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATAcc                                                                                    >  1:2910144/1‑93 (MQ=255)
     ttCTGCCACGAAGAACGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATAcc                                                                                    >  1:2870587/1‑93 (MQ=255)
     ttCTGCCACGAAGAACGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATAcc                                                                                    >  1:2198710/1‑93 (MQ=255)
           cACGAAGAACGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGGGCTACCACCAGTTTGTCTTCTTCAGACAGCTAACCCATACCCAGGAt                                                                              >  2:1761327/1‑93 (MQ=255)
           cACGAAGAACGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGAt                                                                              >  1:4531619/1‑93 (MQ=255)
           cACGAAGAACGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGAt                                                                              >  2:2737786/1‑93 (MQ=255)
           cACGAAGAACGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGAt                                                                              >  2:280918/1‑93 (MQ=255)
            aCGAAGAACGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATg                                                                             <  2:2202547/93‑1 (MQ=255)
             cGAAGAACGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATgg                                                                            >  1:1201521/1‑93 (MQ=255)
             cGAAGAACGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATgg                                                                            <  1:4283635/93‑1 (MQ=255)
             cGAAGAACGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATgg                                                                            <  2:3589300/93‑1 (MQ=255)
             cGAAGAACGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATgg                                                                            <  2:4413209/93‑1 (MQ=255)
              gaagaaCGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGc                                                                           <  2:4019133/93‑1 (MQ=255)
              gaagaaCGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGc                                                                           >  1:668999/1‑93 (MQ=255)
                agaaCGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCga                                                                         <  1:1641446/93‑1 (MQ=255)
                 gaaCGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCgat                                                                        <  2:2098035/93‑1 (MQ=255)
                 gaaCGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCgat                                                                        <  1:1777236/93‑1 (MQ=255)
                  aaCGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGTTGGCgatg                                                                       >  2:2357476/1‑93 (MQ=255)
                  aaCGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCgatg                                                                       >  1:3356064/1‑93 (MQ=255)
                   aCGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCgatga                                                                      >  1:542290/1‑93 (MQ=255)
                   aCGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCgatga                                                                      <  1:3018689/93‑1 (MQ=255)
                   aCGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCgatga                                                                      <  1:3915044/93‑1 (MQ=255)
                     ggCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGatgatg                                                                    <  2:2327366/93‑1 (MQ=255)
                      gCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGTTGGCGATGATGt                                                                   <  2:4303886/93‑1 (MQ=255)
                      gCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGt                                                                   <  2:3551446/93‑1 (MQ=255)
                        tGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCt                                                                 <  1:3762833/93‑1 (MQ=255)
                         gggACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCtt                                                                >  1:4114402/1‑93 (MQ=255)
                          ggACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCttt                                                               <  1:911699/93‑1 (MQ=255)
                           gACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTc                                                              <  2:1905787/93‑1 (MQ=255)
                           gACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTc                                                              <  2:2042749/93‑1 (MQ=255)
                            aCAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCa                                                             >  2:488341/1‑93 (MQ=255)
                             cAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAg                                                            <  2:2891854/93‑1 (MQ=255)
                              aGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGt                                                           <  2:216733/93‑1 (MQ=255)
                               ggAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGtt                                                          <  2:133983/93‑1 (MQ=255)
                               ggAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGtt                                                          >  1:2457392/1‑93 (MQ=255)
                                gAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGACGGCGATGATGTCTTTCAGTTc                                                         <  1:1319326/93‑1 (MQ=255)
                                 aaGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTcc                                                        >  2:1733601/1‑93 (MQ=255)
                                     gcTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGAt                                                    <  1:2714846/93‑1 (MQ=255)
                                     gcTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCGTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGAt                                                    <  1:1698946/93‑1 (MQ=255)
                                      cTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATa                                                   >  2:1386599/1‑93 (MQ=255)
                                       tGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACATTTCACCCATCCCCAGGATGGCGATGATGTCTTTCAGTTCCTGATaa                                                  <  2:4523662/93‑1 (MQ=255)
                                       tGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATaa                                                  >  2:338346/1‑93 (MQ=255)
                                       tGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATaa                                                  >  1:3109538/1‑93 (MQ=255)
                                       tGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATaa                                                  <  2:2484274/93‑1 (MQ=255)
                                        ggATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAAc                                                 <  2:1469907/93‑1 (MQ=255)
                                        ggATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAAc                                                 >  1:4311014/1‑93 (MQ=255)
                                         gATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACg                                                >  1:72242/1‑93 (MQ=255)
                                          aTCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGt                                               >  2:668156/1‑93 (MQ=255)
                                          aTCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGt                                               <  2:195039/93‑1 (MQ=255)
                                           tCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGtt                                              >  1:2748466/1‑93 (MQ=255)
                                           tCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGtt                                              <  2:1640564/93‑1 (MQ=255)
                                           tCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGtt                                              >  2:4070132/1‑93 (MQ=255)
                                             ttACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGc                                            >  1:3448970/1‑93 (MQ=255)
                                             ttACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGc                                            <  2:3469132/93‑1 (MQ=255)
                                              tACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATAGCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCa                                           <  1:3293104/93‑1 (MQ=255)
                                              tACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCa                                           >  2:2524787/1‑93 (MQ=255)
                                               aCGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAg                                          <  2:863491/93‑1 (MQ=255)
                                               aCGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAg                                          <  2:3147615/93‑1 (MQ=255)
                                                 gAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGa                                        <  2:3240983/93‑1 (MQ=255)
                                                  aGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGAt                                       <  1:948929/93‑1 (MQ=255)
                                                      cgcgCTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGAc                                   <  1:1732544/93‑1 (MQ=255)
                                                       gcgcTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACt                                  >  1:1022803/1‑93 (MQ=255)
                                                         gcTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGa                                >  1:639567/1‑93 (MQ=255)
                                                         gcTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGa                                >  2:3314310/1‑93 (MQ=255)
                                                          cTACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGaa                               <  2:975016/93‑1 (MQ=255)
                                                           tACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAAc                              >  2:2647442/1‑93 (MQ=255)
                                                           tACCACCAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAAc                              <  1:2006695/93‑1 (MQ=255)
                                                            accaccAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACg                             >  2:4148918/1‑93 (MQ=255)
                                                            accaccAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACg                             >  2:2371198/1‑93 (MQ=255)
                                                            accaccAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACg                             >  2:3052839/1‑93 (MQ=255)
                                                            accaccAGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACg                             >  2:1710456/1‑93 (MQ=255)
                                                              caccaGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGcc                           >  2:682099/1‑93 (MQ=255)
                                                               accaGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCa                          <  2:2408775/93‑1 (MQ=255)
                                                               accaGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCa                          >  1:185320/1‑93 (MQ=255)
                                                               accaGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCa                          <  2:3356277/93‑1 (MQ=255)
                                                               accaGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCa                          <  2:2924733/93‑1 (MQ=255)
                                                                ccaTTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATTACGTTGCAGGCTGGACTGAACGGCAc                         >  1:2244589/1‑93 (MQ=255)
                                                                ccaGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCAc                         <  2:2990005/93‑1 (MQ=255)
                                                                  aGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCAcgc                       >  2:2271272/1‑93 (MQ=255)
                                                                  aGTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCAcgc                       <  2:3897894/93‑1 (MQ=255)
                                                                   gTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCAcgcg                      >  1:3629284/1‑93 (MQ=255)
                                                                   gTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCAcgcg                      >  1:76401/1‑93 (MQ=255)
                                                                   gTTTGTCTTCTTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCAcgcg                      >  2:2471389/1‑93 (MQ=255)
                                                                        tcttctTCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCACGCGCGgtg                 >  1:4643768/1‑93 (MQ=255)
                                                                         cttcttCAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCACGCGCGgtgt                >  2:1546315/1‑93 (MQ=255)
                                                                           tcttcAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCACGCGCGGTGTCg              >  1:2381278/1‑93 (MQ=255)
                                                                             ttcAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCACGCGCGGTGTCGTa            >  1:1970197/1‑93 (MQ=255)
                                                                             ttcAGACAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCACGCGCGGTGTCGTa            >  1:1290620/1‑93 (MQ=255)
                                                                                  aCAGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCACGCGCGGTGTCGTAGTGtt       <  1:1525652/93‑1 (MQ=255)
                                                                                    aGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCACGCGCGGTGTCGTAGTGTTcc     <  1:4400466/93‑1 (MQ=255)
                                                                                    aGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCACGCGCGGTGTCGTAGTGTTcc     <  1:2981297/93‑1 (MQ=255)
                                                                                    aGTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCACGCGCGGTGTCGTAGTGTTcc     >  1:3076417/1‑93 (MQ=255)
                                                                                     gTTCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCACGCGCGGTGTCGTAGTGTTCCt    <  1:1372291/93‑1 (MQ=255)
                                                                                      ttCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCACGCGCGGTGTCGTAGTGTTCCTg   >  2:1179909/1‑93 (MQ=255)
                                                                                      ttCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCACGCGCGGTGTCGTAGTGTTCCTg   >  2:127518/1‑93 (MQ=255)
                                                                                       tCACCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCACGCGCGGTGTCGTAGTGTTCCTGa  <  2:1881017/93‑1 (MQ=255)
                                                                                          |                                                                                         
AATACTTCTGCCACGAAGAACGGCTGGGACAGGAAGCGCTGGATCTTACGAGCACGCGCTACCACCAGTTTGTCTTCTTCAGACAGTTCATCCATACCCAGGATGGCGATGATGTCTTTCAGTTCCTGATAACGTTGCAGGATGGACTGAACGCCACGCGCGGTGTCGTAGTGTTCCTGA  >  NC_000913/3913772‑3913951

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