Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,822,045 G→A L169L (CTG→CTA b1741 → putative excinuclease subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,822,0450GA91.5% 702.4 / 29.1 237L169L (CTG→CTAb1741putative excinuclease subunit
Reads supporting (aligned to +/- strand):  ref base G (13/7);  new base A (116/99);  total (130/107)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

Only 100 of 251 total aligned reads displayed.

CCCCCAACCTGTTTGGCCTGTTTGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTGCTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAGCATGCTGCGGTAAAGAGAGCCA  >  NC_000913/1821954‑1822134
                                                                                           |                                                                                         
 ccccAACCTGTTTGGCCTGTTTGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCctact                                                                                         >  1:872271/1‑93 (MQ=255)
 ccccAACCTGTTTGGCCTGTTTGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCctact                                                                                         <  2:1068252/93‑1 (MQ=255)
 ccccAACCTGTTTGGCCTGTTTGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCctact                                                                                         >  1:1004240/1‑93 (MQ=255)
  cccAACCTGTTTGGCCTGTTTGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCctactg                                                                                        >  2:6455763/1‑93 (MQ=255)
  cccAACCTGTTTGGCCTGTTTGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCctactg                                                                                        >  1:298311/1‑93 (MQ=255)
  cccAACCTGTTTGGCCTGTTTGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCctactg                                                                                        >  1:2243664/1‑93 (MQ=255)
  cccAACCTGTTTGGCCTGTTTACCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCctgctg                                                                                        >  2:4230947/1‑93 (MQ=255)
   ccAACCTGTTTGGCCTGTTTGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTgg                                                                                       >  1:4743607/1‑93 (MQ=255)
    cAACCTGTTTGGCCTGTTTGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTggg                                                                                      >  2:2352013/1‑93 (MQ=255)
     aaCCTGTTTGGCCTGTTTGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGa                                                                                     >  1:7927687/1‑93 (MQ=255)
     aaCCTGTTTGGCCTGTTTGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGa                                                                                     <  2:2531237/93‑1 (MQ=255)
           tttGGCCTGTTTGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGaa                                                                               >  2:3178436/1‑93 (MQ=255)
            ttGGCCTGTTTGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAAc                                                                              <  1:2049171/93‑1 (MQ=255)
            ttGGCCTGTTTGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAAc                                                                              <  1:8321007/93‑1 (MQ=255)
            ttGACCTGTTTGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAAc                                                                              <  1:9396741/93‑1 (MQ=255)
             tGGCCTGTTTGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAAcc                                                                             <  2:5690669/93‑1 (MQ=255)
              gtccTTTTTGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGTCCTTCTGGGACTGGTACCg                                                                            >  2:8977544/3‑93 (MQ=255)
               gCCTGTTTGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGGCTGGAACCGt                                                                           >  2:3124441/1‑93 (MQ=255)
               gCCTGTTTGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGt                                                                           <  1:1445273/93‑1 (MQ=255)
                ccTGTTTGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGtt                                                                          <  1:9163661/93‑1 (MQ=255)
                 cTGTTTGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTa                                                                         <  2:8518731/93‑1 (MQ=255)
                   gTTTGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGGCTGGAACCGTTAAg                                                                       >  2:5646207/1‑93 (MQ=255)
                    tttGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGt                                                                      >  2:8855298/1‑93 (MQ=255)
                    tttGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGt                                                                      >  2:6970132/1‑93 (MQ=255)
                    tttGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGt                                                                      <  2:3428319/93‑1 (MQ=255)
                    tttGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTc                                                                     >  2:1014294/1‑93 (MQ=255)
                          aaTAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGt                                                                >  1:1085765/1‑93 (MQ=255)
                          aaTAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGt                                                                >  1:1364530/1‑93 (MQ=255)
                            tAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCg                                                              >  2:7110076/1‑93 (MQ=255)
                              ggCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTg                                                            <  2:8439631/93‑1 (MQ=255)
                               gcgcgcAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGc                                                           >  2:576331/1‑93 (MQ=255)
                               gcgcgcAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGc                                                           <  2:7542141/93‑1 (MQ=255)
                               gcgcgcAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGc                                                           <  2:4459697/93‑1 (MQ=255)
                                cgcgcAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCa                                                          >  1:9263649/1‑93 (MQ=255)
                                 gcgcAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCAt                                                         <  2:4653057/93‑1 (MQ=255)
                                   gcAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGt                                                       <  1:8550176/93‑1 (MQ=255)
                                   gcAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTCGAACCGTTAAGTCGCGGTCGTGCATGt                                                       >  1:1510064/1‑93 (MQ=255)
                                    cAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGtt                                                      >  2:877204/1‑93 (MQ=255)
                                     aGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGttt                                                     <  2:9344814/93‑1 (MQ=255)
                                      gatttGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGtttt                                                    >  1:3743435/3‑93 (MQ=255)
                                       cTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTc                                                   <  1:4621005/93‑1 (MQ=255)
                                        tttGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTGCTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCg                                                  >  1:6976830/1‑93 (MQ=255)
                                        tttGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCg                                                  >  1:7287350/1‑93 (MQ=255)
                                           gCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTGCTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTc                                               >  1:2335001/1‑93 (MQ=255)
                                           gCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTc                                               >  2:543105/1‑93 (MQ=255)
                                           gCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTc                                               >  1:3410061/1‑93 (MQ=255)
                                           gCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTc                                               <  2:4053466/93‑1 (MQ=255)
                                             aaGCATTGCAGACCATCGCTGATGAACAAAAACTTTGGTATGGCCTACTGGGACTGGAACGGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAg                                             >  1:6117509/1‑93 (MQ=255)
                                                cATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGGCTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAgcgc                                          <  1:3039549/93‑1 (MQ=255)
                                                cATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAgcgc                                          >  2:367032/1‑93 (MQ=255)
                                                cATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAgcgc                                          <  1:2697309/93‑1 (MQ=255)
                                                 aTTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCt                                         >  1:5962292/1‑93 (MQ=255)
                                                    gCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTaaa                                      >  1:5018034/1‑93 (MQ=255)
                                                     cAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTaaaa                                     <  1:3233859/93‑1 (MQ=255)
                                                      aGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTGCTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAAc                                    >  1:7611179/1‑93 (MQ=255)
                                                      aGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAAc                                    >  1:1816339/1‑93 (MQ=255)
                                                      aGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAAc                                    >  2:7709049/1‑93 (MQ=255)
                                                       gACCATCGCTGATGAAGAAAAACTTTGTTATGGCCTGCTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACg                                   >  2:6706443/1‑93 (MQ=255)
                                                       gACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACg                                   >  1:9101189/1‑93 (MQ=255)
                                                       gACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACg                                   <  2:138012/93‑1 (MQ=255)
                                                       gACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACg                                   >  2:1701923/1‑93 (MQ=255)
                                                       gACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACg                                   >  1:1273592/1‑93 (MQ=255)
                                                       gACCATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGGGCATGTTTTCGTCCACCGCAAAAAct                                   >  1:4362789/1‑92 (MQ=255)
                                                         ccATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGtt                                 <  1:9864920/93‑1 (MQ=255)
                                                         ccATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGtt                                 >  1:7468135/1‑93 (MQ=255)
                                                         ccATCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGtt                                 >  1:6078777/1‑93 (MQ=255)
                                                            tCGCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTgcg                              >  1:6991553/1‑93 (MQ=255)
                                                             cGCTGATGAACAAAAACTTTGTTATGGCCTGCTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTgcgc                             <  2:8037404/93‑1 (MQ=255)
                                                              gCTGATGAACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGcc                            >  2:2092711/1‑93 (MQ=255)
                                                                tgatgaACAAAAACTTTGTTATGGCCTACTGGGGCTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCgg                          >  2:7450407/1‑93 (MQ=255)
                                                                tgatgaACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCgg                          <  2:9563419/93‑1 (MQ=255)
                                                                 tttgaACAAATACTTTGTTTTGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTATCGTTTAGCGCTTATACGTTGCGCCTGa                         >  1:8947394/3‑93 (MQ=255)
                                                                 gatgaACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGa                         >  1:2739817/1‑93 (MQ=255)
                                                                  atgaACAAAAACTTTGTTATGGCCTACTGGGGCTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAg                        >  1:28032/1‑93 (MQ=255)
                                                                  atgaACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAg                        >  2:2596210/1‑93 (MQ=255)
                                                                  atgaACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAg                        <  2:9200064/93‑1 (MQ=255)
                                                                   tgaACAAAAACTTTGTTATGGCCTGCTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAGc                       >  2:8913901/1‑93 (MQ=255)
                                                                   tgaACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAGc                       >  2:6153785/1‑93 (MQ=255)
                                                                    gaACAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAGCa                      >  2:4843008/1‑93 (MQ=255)
                                                                       cAAAAACTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAGCAtgc                   <  1:3646762/93‑1 (MQ=255)
                                                                         aaaaCTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAGCAtgctg                 >  2:7835348/1‑93 (MQ=255)
                                                                         aaaaCTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAGCAtgctg                 >  1:5869004/1‑93 (MQ=255)
                                                                           aaCTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAGCATGCTGCg               >  1:5377485/1‑93 (MQ=255)
                                                                            aCTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAGCATGCTGCgg              >  2:1961545/1‑93 (MQ=255)
                                                                             cTTTGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAGCATGCTGCGGt             <  1:7739016/93‑1 (MQ=255)
                                                                               ttGTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAGCATGCTGCGGTaa           <  2:5034952/93‑1 (MQ=255)
                                                                                 gTTATGGCCTGCTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAGCATGCTGCGGTAAag         >  1:8159591/1‑93 (MQ=255)
                                                                                 gTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAGCATGCTGCGGTAAag         <  1:9543197/93‑1 (MQ=255)
                                                                                 gTTATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAGCATGCTGCGGTAAag         >  1:4229521/1‑93 (MQ=255)
                                                                                  ttATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCTCCGGAGCATGCTGCGGTAAaga        <  1:6656066/93‑1 (MQ=255)
                                                                                  ttATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAGCATGCTGCGGTAAaga        <  1:2643553/93‑1 (MQ=255)
                                                                                  ttATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAGCATGCTGCGGTAAaga        <  1:6608612/93‑1 (MQ=255)
                                                                                   tATGGCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAGCATGCTGCGGTAAagag       >  2:1197057/1‑93 (MQ=255)
                                                                                    aTGGCCTGCTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAGCATGCTGCGGTAAagaga      >  2:6509119/1‑93 (MQ=255)
                                                                                      ggCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAGCATGCTGCGGTAAAGAGAGc    >  2:2367125/1‑93 (MQ=255)
                                                                                       gCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAGCATGCTGCGGTAAAGAGAGcc   <  1:6283152/93‑1 (MQ=255)
                                                                                       gCCTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAGCATGCTGCGGTAAAGAGAGcc   >  2:2212342/1‑93 (MQ=255)
                                                                                        ccTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAGCATGCTGCGGTAAAGAGAGCCa  >  2:8107616/1‑93 (MQ=255)
                                                                                        ccTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAGCATGCTGCGGTAAAGAGAGCCa  >  1:9955745/1‑93 (MQ=255)
                                                                                        ccTACTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCCCCGGAGCATGCTGCGGTAAAGAGAGCCa  <  1:6068335/93‑1 (MQ=255)
                                                                                           |                                                                                         
CCCCCAACCTGTTTGGCCTGTTTGCCAATAGGCGCGCAGCTTTGCAAGCATTGCAGACCATCGCTGATGAACAAAAACTTTGTTATGGCCTGCTGGGACTGGAACCGTTAAGTCGCGGTCGTGCATGTTTTCGTTCAGCGCTAAAACGTTGCGCCGGAGCATGCTGCGGTAAAGAGAGCCA  >  NC_000913/1821954‑1822134

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