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,5020TC94.7% 382.5 / 7.9 132D314G (GAT→GGT) cysNATP‑sulfurylase subunit 1
Reads supporting (aligned to +/- strand):  ref base T (3/4);  new base C (68/57);  total (71/61)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 138 total aligned reads displayed.

TCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACC  >  NC_000913/2872413‑2872593
                                                                                         |                                                                                           
tccatACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCg                                                                                          >  1:3349461/1‑93 (MQ=255)
tccatACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCg                                                                                          >  1:1161319/1‑93 (MQ=255)
  catACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTc                                                                                        >  1:4669367/1‑93 (MQ=255)
  catACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTc                                                                                        >  1:2688753/1‑93 (MQ=255)
  catACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTc                                                                                        >  1:2215944/1‑93 (MQ=255)
  catACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTc                                                                                        <  2:4414491/93‑1 (MQ=255)
  catACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTc                                                                                        >  1:1363073/1‑93 (MQ=255)
   atACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCa                                                                                       >  2:4038214/1‑93 (MQ=255)
       cacaTCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCAc                                                                                   >  2:1649020/1‑93 (MQ=255)
          aTCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAg                                                                                >  2:17045/1‑93 (MQ=255)
          aTCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGACACGGCTGATGTCGATCTCACACGTCAGCACCAg                                                                                >  2:3101034/1‑93 (MQ=255)
           tCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAgg                                                                               >  2:2142817/1‑93 (MQ=255)
           tCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAgg                                                                               >  1:1269477/1‑93 (MQ=255)
           tCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAgg                                                                               >  2:3451602/1‑93 (MQ=255)
             cACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAgggg                                                                             >  1:4501899/1‑92 (MQ=255)
             cACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGt                                                                             <  2:2435306/93‑1 (MQ=255)
             cACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGt                                                                             <  1:2834719/93‑1 (MQ=255)
              aCCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTg                                                                            >  2:719276/1‑93 (MQ=255)
               ccGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCATCGTGa                                                                           >  1:1747610/1‑93 (MQ=255)
               ccGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGGGa                                                                           >  2:3486346/1‑93 (MQ=255)
                cGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGAt                                                                          >  1:57167/1‑93 (MQ=255)
                 gACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGTAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATc                                                                         >  1:3359789/1‑93 (MQ=255)
                   cgcgcgCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGc                                                                       <  2:200174/93‑1 (MQ=255)
                    gcgcTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCt                                                                      <  2:515638/93‑1 (MQ=255)
                      ccactGCCCCGCCGGTAACGCTTCGTCGGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTtc                                                                    <  1:3861191/90‑1 (MQ=255)
                        tctGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTtctc                                                                  >  1:1126273/1‑93 (MQ=255)
                        tctGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTtctc                                                                  >  2:1482155/1‑93 (MQ=255)
                        tctGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTtctc                                                                  >  1:4085391/1‑93 (MQ=255)
                        cctGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTtctc                                                                  <  2:2404433/92‑1 (MQ=255)
                         ctGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTcc                                                                 <  2:3838233/93‑1 (MQ=255)
                         ctGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTcc                                                                 <  1:2636795/93‑1 (MQ=255)
                             aCCGCCGGTAACGCTTCGTCTGCCGCCAGCAACAGATCGCCACGGCTGATGTCGATCTCACCAGTCAGCACCAGGGTGATCGCTTCTCCGGCa                                                             >  1:544279/1‑93 (MQ=255)
                              ccgccgGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCaa                                                            >  2:1590252/1‑93 (MQ=255)
                              ccgccgGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCaa                                                            <  2:3579647/93‑1 (MQ=255)
                                gccgGTAACGCTTCGTCTGCCGCCAGCATCAGATCGCCACGGCTGATGGCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAg                                                          <  2:3700647/93‑1 (MQ=255)
                                gccgGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAg                                                          <  1:4522758/93‑1 (MQ=255)
                                gccgGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAg                                                          <  2:1234221/93‑1 (MQ=255)
                                gccgGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAg                                                          <  1:1929062/93‑1 (MQ=255)
                                  cgGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGc                                                        >  2:3617112/1‑93 (MQ=255)
                                    gTAACGCTTCGTCTGCCGCCCGCCGCAGATAGCCACGGGTGGTGTCGATCTCACCCGCCAGCACCAGCGTGCTAGCCTCCCCCGCAAAGGcct                                                      >  2:1115371/1‑91 (MQ=255)
                                    gTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGctt                                                      >  2:1115236/1‑93 (MQ=255)
                                     tAACGCTTCGTCTGCCGTCAGCAGCAGATAGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGcttc                                                     >  2:3327789/1‑93 (MQ=255)
                                     tAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGcttc                                                     <  2:3408659/93‑1 (MQ=255)
                                     tAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGcttc                                                     <  2:492407/93‑1 (MQ=255)
                                      aaCGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGcttct                                                    <  2:350577/93‑1 (MQ=255)
                                      aaCGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGcttct                                                    <  2:1455280/93‑1 (MQ=255)
                                      aaCGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGcttct                                                    <  2:866623/93‑1 (MQ=255)
                                        cGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCttcttc                                                  <  2:2245866/93‑1 (MQ=255)
                                         gCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTcg                                                 >  2:2933472/1‑93 (MQ=255)
                                         gCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGCGATCGCTTCTCCGGCAAAGGCTTCTTcg                                                 >  1:2205726/1‑93 (MQ=255)
                                         gCTTCGTCTGCCGACAGCAGCAGATCGCCACGGCTGATGTCAATCTCACCCGTCAGCACCAGGCTGATAGCTTATCCGGCAAAGGCTTCTTcg                                                 >  2:3136016/1‑93 (MQ=255)
                                          cTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTcgc                                                >  1:3327533/1‑93 (MQ=255)
                                          cTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTcgc                                                >  2:2898743/1‑93 (MQ=255)
                                             cGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGAt                                             <  1:1099118/93‑1 (MQ=255)
                                             cGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGAt                                             <  1:3040473/93‑1 (MQ=255)
                                             cGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGAt                                             >  1:3939526/1‑93 (MQ=255)
                                              gTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATc                                            >  2:3041381/1‑93 (MQ=255)
                                               tCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCa                                           <  1:2836017/93‑1 (MQ=255)
                                               tCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGAGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCa                                           >  1:2583479/1‑93 (MQ=255)
                                                 tGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCAcc                                         <  2:2410287/93‑1 (MQ=255)
                                                  gccgccAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCa                                        <  1:4361001/93‑1 (MQ=255)
                                                  gccgccAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCa                                        <  1:3097858/93‑1 (MQ=255)
                                                     gccAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCa                                     >  2:1286285/1‑93 (MQ=255)
                                                      ccAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCaa                                    <  2:865411/93‑1 (MQ=255)
                                                      ccAGCAGCAGATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCaa                                    >  1:1708443/1‑93 (MQ=255)
                                                       cagcagcagATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCCACATCaca                                   >  1:979906/1‑91 (MQ=255)
                                                        agcagcagATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCaaaa                                  >  2:4522664/1‑93 (MQ=255)
                                                        agcagcagATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCaaaa                                  <  1:529389/93‑1 (MQ=255)
                                                         gcagcagATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAg                                 <  1:3805522/93‑1 (MQ=255)
                                                         gcagcagATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAg                                 >  2:1688818/1‑93 (MQ=255)
                                                            gcagATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGCTCACCATCAAAAGTCa                              >  1:2363732/1‑93 (MQ=255)
                                                            gcagATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCa                              >  1:2053364/1‑93 (MQ=255)
                                                            gcagATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCa                              >  2:1951758/1‑93 (MQ=255)
                                                            gcagATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCa                              >  2:193935/1‑93 (MQ=255)
                                                             cagATCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCAc                             <  1:1483393/93‑1 (MQ=255)
                                                              agATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTAACg                            >  2:3728901/1‑93 (MQ=255)
                                                               aaTCGCCCCGCCTGATGTCGATCTCACCCCTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGa                           <  1:2974065/92‑1 (MQ=255)
                                                                 tCGCCGCGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATc                         <  1:2494683/93‑1 (MQ=255)
                                                                 tCGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATc                         <  2:4004530/93‑1 (MQ=255)
                                                                  cGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTTATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATcc                        <  2:2329467/93‑1 (MQ=255)
                                                                  cGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATcc                        <  2:123963/93‑1 (MQ=255)
                                                                  cGCCACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATcc                        <  2:4481517/93‑1 (MQ=255)
                                                                  cGCCACGGCTCATGTCGATCTTACCCGTCAGCACCAGGGTGCTCGCTTCTCCGGTAAAGGCTTCTTCCCGATCACCATCAAAAGTCACGATcc                        <  2:4295936/93‑1 (MQ=255)
                                                                    ccACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCcgc                      <  1:4105541/93‑1 (MQ=255)
                                                                    ccACGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCcgc                      <  1:864013/93‑1 (MQ=255)
                                                                      aCGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCcgcgc                    >  1:3518260/1‑93 (MQ=255)
                                                                      aCGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCcgcgc                    >  1:3561683/1‑93 (MQ=255)
                                                                       cGGCTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCcgcgcg                   >  1:2709401/1‑93 (MQ=255)
                                                                          cTGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGACCCGCGCGACg                >  2:1305470/1‑93 (MQ=255)
                                                                           tGATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGt               <  2:1444466/93‑1 (MQ=255)
                                                                           tGATGTCGATCTCACCCGTCAGCACCAGGCTGATCGCTTCTCCGGCAAAAGCTTCTTCGCGATCACCATCAAAAGTCACGATCCCCGCGACGt               >  1:4442783/1‑93 (MQ=255)
                                                                            gATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGtt              <  2:2115223/93‑1 (MQ=255)
                                                                            gATGTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGtt              <  2:1571603/93‑1 (MQ=255)
                                                                               gTCGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGa           <  2:1148323/93‑1 (MQ=255)
                                                                                 cGATCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGAtt         <  2:4326603/93‑1 (MQ=255)
                                                                                   aTCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTcc       <  1:3799241/93‑1 (MQ=255)
                                                                                   aTCTCACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTcc       >  2:133931/1‑93 (MQ=255)
                                                                                    tctcACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCca      <  2:3665189/93‑1 (MQ=255)
                                                                                      tcACCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCcaca    >  1:2726100/1‑93 (MQ=255)
                                                                                        aTCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGCCGTTTGATTCCACAcc  >  2:2330603/1‑93 (MQ=255)
                                                                                         |                                                                                           
TCCATACCACATCCACCGACGCGCTCTGCACCGCCGGTAACGCTTCGTCTGCCGCCAGCAGCAGATCGCCACGGCTGATGTCGATCTCATCCGTCAGCACCAGGGTGATCGCTTCTCCGGCAAAGGCTTCTTCGCGATCACCATCAAAAGTCACGATCCGCGCGACGTTTGATTCCACACC  >  NC_000913/2872413‑2872593

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