Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,872,502 T→C D314G (GAT→GGT)  cysN ← ATP‑sulfurylase subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,872,5020TC99.5% 694.0 / NA 213D314G (GAT→GGT) cysNATP‑sulfurylase subunit 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (109/103);  minor base A (1/0);  total (110/103)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 226 total aligned reads displayed.

ATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACA  >  NC_000913/2872412‑2872591
                                                                                          |                                                                                         
atccatACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCAccc                                                                                         >  2:1004065/1‑93 (MQ=255)
atccatACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCAccc                                                                                         <  2:1928076/93‑1 (MQ=255)
atccatACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCAccc                                                                                         >  1:2047852/1‑93 (MQ=255)
atccatACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGATGATGTCGATCTCAccc                                                                                         >  1:5631590/1‑93 (MQ=255)
 tccatACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCg                                                                                        >  1:780016/1‑93 (MQ=255)
   catACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTc                                                                                      <  1:1624458/93‑1 (MQ=255)
    atACCACATCCACCGACGCGCTCTGCCCCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCa                                                                                     <  2:6614692/93‑1 (MQ=255)
    atACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCa                                                                                     <  1:3838335/93‑1 (MQ=255)
    atACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCa                                                                                     >  2:4866827/1‑93 (MQ=255)
     tACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAg                                                                                    >  2:6683240/1‑93 (MQ=255)
       ccACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGcc                                                                                  >  1:1485682/1‑92 (MQ=255)
       ccACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCa                                                                                  >  1:2719671/1‑93 (MQ=255)
        cacaTCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCAc                                                                                 <  1:5414960/93‑1 (MQ=255)
         acaTCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCAcc                                                                                >  2:3017499/1‑93 (MQ=255)
         acaTCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCAcc                                                                                >  1:64666/1‑93 (MQ=255)
          caTCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCa                                                                               <  1:1001012/93‑1 (MQ=255)
           aTCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAg                                                                              <  2:2929234/93‑1 (MQ=255)
            tCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGAGCTCACCCGTCAGCAGCAgg                                                                             >  2:5266254/1‑93 (MQ=255)
             ccACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAggg                                                                            <  1:4493387/93‑1 (MQ=255)
              cgccGGCGCCCCCTGCCCCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGt                                                                           <  2:4412732/91‑1 (MQ=255)
               cccGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTg                                                                          <  2:3328494/92‑1 (MQ=255)
                ccgccgCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTCATGTCGATCTCACCCGTCAGCACCAGGGTGa                                                                         <  1:5342546/93‑1 (MQ=255)
                ccGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGa                                                                         >  2:5176565/1‑93 (MQ=255)
                ccGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGa                                                                         >  1:3969013/1‑93 (MQ=255)
                ccGACGCCCTCTGCCCCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGa                                                                         <  2:5317273/93‑1 (MQ=255)
                 cGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGAt                                                                        >  1:1352514/1‑93 (MQ=255)
                 cGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGAt                                                                        <  2:6002091/93‑1 (MQ=255)
                  gACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATc                                                                       >  1:5918433/1‑93 (MQ=255)
                   aCGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCg                                                                      >  2:6557491/1‑93 (MQ=255)
                    cgcgCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGc                                                                     >  1:4321257/1‑93 (MQ=255)
                    cgcgCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGc                                                                     <  1:330265/93‑1 (MQ=255)
                     gcgcTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCt                                                                    <  2:3878744/93‑1 (MQ=255)
                     gcgcTCTCCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCt                                                                    <  1:1421756/93‑1 (MQ=255)
                      cgcTCTGCGCCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCtt                                                                   <  2:1505178/93‑1 (MQ=255)
                       gcTCTGCACCGCCGGTAACGGTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTtc                                                                  <  1:3511677/93‑1 (MQ=255)
                         tctGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTtctc                                                                >  1:4131448/1‑93 (MQ=255)
                         tctGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTtctc                                                                >  1:3343379/1‑93 (MQ=255)
                          ctGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTcc                                                               <  2:6204057/93‑1 (MQ=255)
                          ctGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTcc                                                               <  2:4488558/93‑1 (MQ=255)
                            gCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCgg                                                             >  1:5827855/1‑93 (MQ=255)
                             cACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGc                                                            >  1:5733350/1‑93 (MQ=255)
                              aCCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCa                                                           >  2:3428477/1‑93 (MQ=255)
                                cgccgGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCaaa                                                         >  2:2571733/1‑93 (MQ=255)
                                cgccgGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCaaa                                                         <  2:3847123/93‑1 (MQ=255)
                                cgccgGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCaaa                                                         <  1:4744201/93‑1 (MQ=255)
                                    ggTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGct                                                     >  2:777541/1‑93 (MQ=255)
                                    ggTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGct                                                     >  2:1090518/1‑93 (MQ=255)
                                    ggTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGct                                                     <  1:5039979/93‑1 (MQ=255)
                                     gTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGCGTGATCGCTTCTCCGGCAAAGGctt                                                    >  1:3591784/1‑93 (MQ=255)
                                      tAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGcttc                                                   <  1:2226492/93‑1 (MQ=255)
                                      tAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGcttc                                                   <  1:1869368/93‑1 (MQ=255)
                                       aaCGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGcttct                                                  <  2:3248900/93‑1 (MQ=255)
                                       aaCGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGcttct                                                  <  1:4854838/93‑1 (MQ=255)
                                       aaCGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGcttct                                                  <  2:4518094/93‑1 (MQ=255)
                                        aCGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAACGGcttctt                                                 <  1:6521733/93‑1 (MQ=255)
                                        aCGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGcttctt                                                 >  2:3652239/1‑93 (MQ=255)
                                        aCGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGcttctt                                                 <  2:4940847/93‑1 (MQ=255)
                                        aCGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGcttctt                                                 <  2:406843/93‑1 (MQ=255)
                                        aCGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGcttctt                                                 >  1:2506819/1‑93 (MQ=255)
                                          gCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTcg                                               >  2:1204240/1‑93 (MQ=255)
                                          gCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTcg                                               >  1:4287805/1‑93 (MQ=255)
                                           cTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTcgc                                              <  2:5806627/93‑1 (MQ=255)
                                           cTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTcgc                                              >  2:4229791/1‑93 (MQ=255)
                                           cTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTcgc                                              <  2:5906521/93‑1 (MQ=255)
                                            ttCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTcgcg                                             >  2:249837/1‑93 (MQ=255)
                                            ttCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTcgcg                                             >  1:6341444/1‑93 (MQ=255)
                                               gTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATc                                          <  1:1277331/93‑1 (MQ=255)
                                               gTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATc                                          >  2:354737/1‑93 (MQ=255)
                                                tCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCa                                         >  2:3984990/1‑93 (MQ=255)
                                                tCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCa                                         >  2:5571385/1‑93 (MQ=255)
                                                 cTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCAc                                        >  2:2612575/1‑93 (MQ=255)
                                                    ccgccAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCAt                                     <  1:2363152/93‑1 (MQ=255)
                                                         agcagcagATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCAGCaaaa                                <  2:1246361/93‑1 (MQ=255)
                                                             gcagATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCa                            >  2:2896051/1‑93 (MQ=255)
                                                             gcagATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCa                            <  2:1327045/93‑1 (MQ=255)
                                                              cagATCGCAACGGCTGATGTCGATCTCTCCCGTCAGCCCCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCAc                           <  1:4061269/93‑1 (MQ=255)
                                                               agATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACg                          >  2:2354601/1‑93 (MQ=255)
                                                               agATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACg                          >  2:183238/1‑93 (MQ=255)
                                                                gATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGa                         <  1:2724385/93‑1 (MQ=255)
                                                                 aTCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGAt                        <  1:120084/93‑1 (MQ=255)
                                                                      cACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCcgcg                   >  2:5120271/1‑93 (MQ=255)
                                                                      cACGGCTGATGTCGATCTCACCCGTCAGCACCAGGCTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCcgcg                   >  1:5282216/1‑93 (MQ=255)
                                                                       ccGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCcgcgc                  <  2:242645/92‑1 (MQ=255)
                                                                       aCGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCcgcgc                  >  1:6704142/1‑93 (MQ=255)
                                                                       aCGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCcgcgc                  >  1:4613443/1‑93 (MQ=255)
                                                                       aCGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCcgcgc                  >  1:597970/1‑93 (MQ=255)
                                                                         ggCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGa                >  1:3683540/1‑93 (MQ=255)
                                                                            tGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGt             <  2:5052921/93‑1 (MQ=255)
                                                                            tGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGt             >  1:3074382/1‑93 (MQ=255)
                                                                            tGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGt             <  1:502963/93‑1 (MQ=255)
                                                                             gATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGtt            <  1:2906537/93‑1 (MQ=255)
                                                                             gATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGtt            <  1:3031927/93‑1 (MQ=255)
                                                                               tGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTg          <  2:4807798/93‑1 (MQ=255)
                                                                               tGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTg          <  1:6227348/93‑1 (MQ=255)
                                                                                 tCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGAt        >  2:3594440/1‑93 (MQ=255)
                                                                                 tCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGAt        >  2:6382418/1‑93 (MQ=255)
                                                                                 tCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGAt        <  2:6457772/93‑1 (MQ=255)
                                                                                   gATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTc      <  2:5436937/93‑1 (MQ=255)
                                                                                    aTCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTcc     >  1:118654/1‑93 (MQ=255)
                                                                                       tcACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCcaca  >  1:865309/1‑93 (MQ=255)
                                                                                          |                                                                                         
ATCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACA  >  NC_000913/2872412‑2872591

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