Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 104,352 A→G E124G (GAG→GGG)  ftsA → ATP‑binding cell division FtsK recruitment protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913104,3520AG98.7% 271.4 / NA 77E124G (GAG→GGG) ftsAATP‑binding cell division FtsK recruitment protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (43/33);  total (44/33)
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

TCTGGCGCTTTCTGGTAAGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGAGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCG  >  NC_000913/104221‑104485
                                                                                                                                   |                                                                                                                                     
tCTGGCGCTTTCTGGTAAGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCgt                                                                                                                                <  1:836457/139‑1 (MQ=255)
    gcgcTTTCTGGTAAGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGtgc                                                                                                                            >  2:923574/1‑139 (MQ=255)
     cgcTTTCTGGTAAGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGtgtg                                                                                                                             <  1:1134421/137‑1 (MQ=255)
     cgcTTTCTGGTAAGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGtgtg                                                                                                                             >  2:1134421/1‑137 (MQ=255)
     cgcTTTCTGGTAAGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGtgct                                                                                                                           >  2:635467/1‑139 (MQ=255)
      gcTTTCTGGTAAGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGATGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATCAGCATCGTGtgctg                                                                                                                          >  2:1355493/1‑139 (MQ=255)
      gcTTTCTGGTAAGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGtgctg                                                                                                                          >  1:61133/1‑139 (MQ=255)
      gcTTTCTGGTAAGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGtgctg                                                                                                                          >  2:285760/1‑139 (MQ=255)
        tttCTGGTAAGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCa                                                                                                                        >  2:760104/1‑139 (MQ=255)
          tCTGGTAAGCACATCAGCTGCCAGAATGAAATTGATATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCAtg                                                                                                                      >  2:483667/1‑139 (MQ=255)
              gTAAGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGAt                                                                                                                  <  1:807849/139‑1 (MQ=255)
                  gCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCg                                                                                                              >  2:1523928/1‑139 (MQ=255)
                  gCACATCAGCTGCCAGAATGAAATTGGTATAGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATc                                                                                                                                  >  2:740482/1‑119 (MQ=255)
                  gCACATCAGCTGCCAGAATGAAATTGGTATAGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATc                                                                                                                                  <  1:740482/119‑1 (MQ=255)
                    acaTCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCa                                                                                                            <  2:23841/139‑1 (MQ=255)
                      aTCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAag                                                                                                          >  1:493119/1‑139 (MQ=255)
                         aGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGt                                                                                                       <  1:830130/139‑1 (MQ=255)
                             gCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAACACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGc                                                                                                   >  1:729534/1‑139 (MQ=255)
                              ccAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGtgtg                                                                                                                             <  1:954029/112‑1 (MQ=255)
                              ccAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGtgtg                                                                                                                             >  2:954029/1‑112 (MQ=255)
                               cAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGa                                                                                                 <  1:285760/139‑1 (MQ=255)
                                   aTGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGa                                                                                             >  1:451857/1‑139 (MQ=255)
                                   aTGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGa                                                                                             <  1:760104/139‑1 (MQ=255)
                                        aTTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATc                                                                                        >  2:546022/1‑139 (MQ=255)
                                           ggTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAgg                                                                                     <  2:1341320/139‑1 (MQ=255)
                                            gTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTAACAGGa                                                                                    >  1:69578/1‑139 (MQ=255)
                                             tATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGaa                                                                                   <  1:515529/139‑1 (MQ=255)
                                                ggTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAggg                                                                                <  2:729534/139‑1 (MQ=255)
                                                 gTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGa                                                                               >  2:962785/1‑139 (MQ=255)
                                                 gTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGa                                                                               >  1:179888/1‑139 (MQ=255)
                                                    ccTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCa                                                                            >  1:1048957/1‑139 (MQ=255)
                                                      tATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTACTGCATGTGATCCCGCAAGAGTa                                                                                                      <  1:1357755/111‑1 (MQ=255)
                                                      tATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTACTGCATGTGATCCCGCAAGAGTa                                                                                                      >  2:1357755/1‑111 (MQ=255)
                                                       aTTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGc                                                                                                             >  2:1570977/1‑103 (MQ=255)
                                                       aTTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGc                                                                                                             <  1:1570977/103‑1 (MQ=255)
                                                            tGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGtgtg                                                                                                                             <  1:1314782/82‑1 (MQ=255)
                                                            tGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGtgtg                                                                                                                             >  2:1314782/1‑82 (MQ=255)
                                                            tGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGc                                                                                                             <  1:1157816/98‑1 (MQ=255)
                                                            tGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGc                                                                                                             >  2:1157816/1‑98 (MQ=255)
                                                                 aagaagTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAgg                                                               <  1:635467/139‑1 (MQ=255)
                                                                  agaagTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGa                                                                                                                   <  2:199538/86‑1 (MQ=255)
                                                                  agaagTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGa                                                                                                                   >  1:199538/1‑86 (MQ=255)
                                                                  agaagTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGa                                                                                             <  1:1361031/108‑1 (MQ=255)
                                                                  agaagTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGa                                                                                             >  2:1361031/1‑108 (MQ=255)
                                                                   gaagTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCg                                                                                                              <  1:596031/90‑1 (MQ=255)
                                                                   gaagTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCg                                                                                                              >  2:596031/1‑90 (MQ=255)
                                                                         aCGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCgg                                                       <  2:505273/139‑1 (MQ=255)
                                                                              agaagaTGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCaa                                                                                                           >  1:1474696/1‑82 (MQ=255)
                                                                              agaagaTGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCaa                                                                                                           <  2:1474696/82‑1 (MQ=255)
                                                                              agaagaTGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTAt                                                                                         >  2:1308015/1‑100 (MQ=255)
                                                                              agaagaTGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTAt                                                                                         <  1:1308015/100‑1 (MQ=255)
                                                                                   aTGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGAt                                             >  1:787864/1‑139 (MQ=255)
                                                                                    tgtgGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATg                                            <  2:451857/139‑1 (MQ=255)
                                                                                     gtgGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGAt                                                                                                >  1:398695/1‑86 (MQ=255)
                                                                                     gtgGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGAt                                                                                                <  2:398695/86‑1 (MQ=255)
                                                                                     gtgGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGc                                           >  1:1232230/1‑139 (MQ=255)
                                                                                         aaaaCGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGc                                       <  2:1048957/139‑1 (MQ=255)
                                                                                         aaaaCGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGc                                       <  2:1232230/139‑1 (MQ=255)
                                                                                         aaaaCGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCAGATGCAGGc                                       <  2:179888/139‑1 (MQ=255)
                                                                                           aaCGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCgtg                                                                                                                               <  1:1244128/49‑1 (MQ=255)
                                                                                           aaCGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCgtg                                                                                                                               >  2:1244128/1‑49 (MQ=255)
                                                                                             cgtcgtCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCTGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCaaaa                                   >  2:956988/1‑139 (MQ=255)
                                                                                             cgtcgtCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCaaaa                                   >  1:52328/1‑139 (MQ=255)
                                                                                             cgtcgtCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCaaaa                                   >  1:794563/1‑139 (MQ=255)
                                                                                             cgtcgtCCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGGATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCaaaa                                   >  2:620846/1‑139 (MQ=255)
                                                                                                  tcCATACCGCGAAATCGGTGCGTGTGCGCGATGGGCATCGTGCGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCa                              >  2:527323/1‑139 (MQ=255)
                                                                                                           cgAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCAtgt                                                                                                                     >  1:738474/1‑43 (MQ=255)
                                                                                                           cgAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCAtgt                                                                                                                     <  2:738474/43‑1 (MQ=255)
                                                                                                           cgAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATcac                     >  1:633648/1‑139 (MQ=255)
                                                                                                           cgAAATCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATcac                     >  1:46262/1‑139 (MQ=255)
                                                                                                             aaaTCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAgg                                                                                     <  2:1582799/73‑1 (MQ=255)
                                                                                                             aaaTCGGTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAgg                                                                                     >  1:1582799/1‑73 (MQ=255)
                                                                                                                   gTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGaa                                                                        <  1:1373206/80‑1 (MQ=255)
                                                                                                                   gTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGaa                                                                        >  2:1373206/1‑80 (MQ=255)
                                                                                                                   gTGCGTGTGCGCGATGGGCATCGTGTGCTGCATGTGATCCCACAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTcaca             >  1:993029/1‑139 (MQ=255)
                                                                                                                         gtgCGCGATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGata       <  2:52328/139‑1 (MQ=255)
                                                                                                                              cgATGGGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAGGTGCACCTGATCACATGTCACAACGATATGGCg  >  1:303079/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                                     
TCTGGCGCTTTCTGGTAAGCACATCAGCTGCCAGAATGAAATTGGTATGGTGCCTATTTCTGAAGAAGAAGTGACGCAAGAAGATGTGGAAAACGTCGTCCATACCGCGAAATCGGTGCGTGTGCGCGATGAGCATCGTGTGCTGCATGTGATCCCGCAAGAGTATGCGATTGACTATCAGGAAGGGATCAAGAATCCGGTAGGACTTTCGGGCGTGCGGATGCAGGCAAAAGTGCACCTGATCACATGTCACAACGATATGGCG  >  NC_000913/104221‑104485

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