Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,202,685 A→G 17.3% E265E (GAA→GAG cca → fused tRNA nucleotidyl transferase/2'3'‑cyclic phosphodiesterase/2'nucleotidase and phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,202,6850AG17.3% 429.7 / 81.0 198E265E (GAA→GAGccafused tRNA nucleotidyl transferase/2'3'‑cyclic phosphodiesterase/2'nucleotidase and phosphatase
Reads supporting (aligned to +/- strand):  ref base A (108/55);  new base G (13/21);  total (121/77)
Fisher's exact test for biased strand distribution p-value = 3.36e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-01

Only 100 of 202 total aligned reads displayed.

AATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCC  >  NC_000913/3202589‑3202784
                                                                                                |                                                                                                   
  tGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTg                                                                                                >  1:3768285/1‑100 (MQ=255)
   gACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAActct                                                                                                 >  1:418227/1‑98 (MQ=255)
     ctctctCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGcc                                                                                             <  1:73409/98‑1 (MQ=255)
     cGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGcc                                                                                             >  2:1939921/1‑100 (MQ=255)
     cGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGcc                                                                                             >  2:112338/1‑100 (MQ=255)
       ctctCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCcg                                                                                            >  1:3899639/1‑99 (MQ=255)
         ctcAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCGGCCAGAGTTCTGGCGGCGt                                                                                         <  2:100027/100‑1 (MQ=38)
         ctcAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCcgcg                                                                                          >  2:2629369/1‑99 (MQ=255)
           cAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCACCTTTATGCCACGATCTCGGTAAAGGGCTGACGCGGCCAGAGTTCTGGCGGCGtca                                                                                       <  1:3358800/100‑1 (MQ=38)
           cAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCGGCCAGAGTTCTGGCGGCGtca                                                                                       <  2:4059247/100‑1 (MQ=38)
           cAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCGGCCAGAGTTCTGGCGGCGtca                                                                                       <  2:2131441/100‑1 (MQ=38)
             aTGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGGCGCGGCCAGAGTTCTGGCGGCGtcagc                                                                                     <  2:3177969/100‑3 (MQ=37)
             aTGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCGGCCAGAGTTCTGGCGGCGtcagc                                                                                     <  2:911888/100‑3 (MQ=38)
             aTGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCGGCCAGAGTTCTGGCGGCGtcagc                                                                                     <  2:656007/100‑3 (MQ=38)
             aTGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCGGCCAGAGTTCTGGCGGCGtcagc                                                                                     <  2:1678802/100‑3 (MQ=38)
             aTGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCGGCCAGAGTTCTGGCGGCGtcagc                                                                                     <  2:156235/100‑3 (MQ=38)
             aTGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtcatc                                                                                     >  2:197693/1‑100 (MQ=255)
              tGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCGGCCAGAGTTCTGGCGGCGtcagcc                                                                                    <  2:2359192/100‑4 (MQ=37)
              tGGCGGCAATGCTGAGTCCGCAGGTCTATGCCTGTCTCGCAGCTTTATGCCACGATCTCGGTAAAGGGCTGACGCGGCCAGGGTTCTGGCGGCGtcagc                                                                                     <  1:211914/99‑3 (MQ=19)
                gcggcgATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCGGCCAGAGTTCTGGCGGCGtcagc                                                                                     <  1:3729244/97‑3 (MQ=37)
                gcggcgATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatg                                                                                  >  1:322992/1‑100 (MQ=255)
                gcggcgATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatg                                                                                  >  1:4255119/1‑100 (MQ=255)
                   gcgATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtc                                                                               <  2:2881021/100‑1 (MQ=255)
                   gcgATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtc                                                                               >  1:1507130/1‑100 (MQ=255)
                   gcgATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtc                                                                               >  1:1589453/1‑100 (MQ=255)
                       tGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtcatgg                                                                           >  2:2032721/1‑100 (MQ=255)
                         cTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGcc                                                                         >  2:1372606/1‑100 (MQ=255)
                          tGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGccc                                                                        >  2:3204248/1‑100 (MQ=255)
                             gTCCGCAGGTCGATGTCCGTTTCGCAGCTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCcggc                                                                     <  1:2665772/100‑1 (MQ=255)
                             gTCCGCAGGTCGATGTCCGTGTCTCAACGTTATGCCACGAGGTCTGTAAAGGGTTGACGCCACCATAACGCTGGCCGCGTCATCATGAGCATGGGCcggc                                                                     >  1:1980614/1‑100 (MQ=255)
                              tCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCcggcg                                                                    <  1:1393497/100‑1 (MQ=255)
                                cGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCggg                                                                  <  2:2533508/100‑1 (MQ=255)
                                 gtaGGTCTTTGTTCGTTTCTCATCTTTATGCCATGTTTTCTGTAAAGGGCTGACGCCGCCAGTACTCTGGCCGCGTCTTCTTTGTCATGGCTCGGCggg                                                                  <  1:1107992/97‑1 (MQ=255)
                                 gCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATGTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgt                                                                 >  1:2350831/1‑100 (MQ=255)
                                 gCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgt                                                                 <  2:4144/100‑1 (MQ=255)
                                 gCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgt                                                                 <  2:552450/100‑1 (MQ=255)
                                  cAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgtg                                                                <  2:415776/100‑1 (MQ=255)
                                  cAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgtg                                                                >  1:3533373/1‑100 (MQ=255)
                                    ggTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGtt                                                              >  1:1704462/1‑100 (MQ=255)
                                    ggTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGtt                                                              >  2:3952457/1‑100 (MQ=255)
                                     ggcaaagTCCGTTCCGCAGCTGTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCggg                                                                  <  1:1720490/89‑1 (MQ=255)
                                        gATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAGCTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGt                                                          >  1:2829659/1‑100 (MQ=255)
                                        gATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGt                                                          >  2:2946844/1‑100 (MQ=255)
                                         aTGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGtt                                                         <  2:3900306/100‑1 (MQ=255)
                                         aTGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGtt                                                         <  2:3507123/100‑1 (MQ=255)
                                            tCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGt                                                      >  2:867865/1‑100 (MQ=255)
                                             ccGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTg                                                     <  1:930105/100‑1 (MQ=255)
                                             ccGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTg                                                     >  1:1619626/1‑100 (MQ=255)
                                             ccGTTTCGCAACTTTATGCCACGATCGCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTg                                                     <  1:2549240/100‑1 (MQ=255)
                                               gTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTATGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGa                                                   >  2:1465211/1‑100 (MQ=255)
                                                 ttCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGaac                                                 >  1:1890910/1‑100 (MQ=255)
                                                   cGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGaacaa                                               >  1:1019643/1‑100 (MQ=255)
                                                       aCTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTAt                                           >  2:1253970/1‑100 (MQ=255)
                                                        cTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATg                                          >  1:747668/1‑100 (MQ=255)
                                                          ttATGCCATGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGcc                                        >  1:474207/1‑100 (MQ=255)
                                                          ttATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGcc                                        <  1:3045960/100‑1 (MQ=255)
                                                            aTGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAg                                      >  2:1779672/1‑100 (MQ=255)
                                                             tGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGc                                     >  2:3955029/1‑100 (MQ=255)
                                                             tGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGc                                     >  2:1523387/1‑100 (MQ=255)
                                                              gCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCg                                    <  1:2522472/100‑1 (MQ=255)
                                                              gCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCg                                    >  1:3325706/1‑100 (MQ=255)
                                                               ccACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCAGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGt                                   >  1:2364082/1‑100 (MQ=255)
                                                                cACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTc                                  >  1:698202/1‑100 (MQ=255)
                                                                  cGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTg                                <  2:281355/100‑1 (MQ=255)
                                                                   gATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTgc                               >  1:490566/1‑100 (MQ=255)
                                                                    aTCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGGGTTAAGTTAGTGGAACAACTATGCCAGCGTCTgc                               >  1:3977813/1‑99 (MQ=255)
                                                                      ctcGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGcgcg                            >  2:3044747/1‑100 (MQ=255)
                                                                      ctcGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGcgcg                            >  2:2586719/1‑100 (MQ=255)
                                                                       tcGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGt                           >  2:1689479/1‑100 (MQ=255)
                                                                        cGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCTCGTCATCATGGCCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTg                          <  1:3933137/100‑1 (MQ=255)
                                                                        cGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTg                          <  2:2515672/100‑1 (MQ=255)
                                                                          gTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGcc                        >  1:1221281/1‑100 (MQ=255)
                                                                              aGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCaaa                     >  2:435401/1‑99 (MQ=255)
                                                                              aGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAAt                    <  2:3436173/100‑1 (MQ=255)
                                                                               gggCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCCAATg                   >  1:856223/1‑100 (MQ=255)
                                                                                 gCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGaa                 <  2:2788372/100‑1 (MQ=255)
                                                                                  cTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGaaa                <  1:3460520/100‑1 (MQ=255)
                                                                                   tGACGCGGCCAGAGTTCTGGCGGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAAt               >  1:2252809/1‑100 (MQ=38)
                                                                                   tGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAAt               <  2:1704469/100‑1 (MQ=255)
                                                                                    gACGCGGCCAGAGTTCTGGCGGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAAtt              >  1:795494/1‑100 (MQ=38)
                                                                                    gACGCGGCCAGAGTTCTGGCGGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAAtt              >  1:3405253/1‑100 (MQ=38)
                                                                                     aCGCGGCCAGAGTTCTGGCGGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTc             >  1:2757773/1‑100 (MQ=38)
                                                                                      cgcggcCAGAGTTCTGGCGGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGTCAGCGTCTGCGCGTGCCAAATGAAATTcg            >  1:3780500/1‑100 (MQ=255)
                                                                                      cgcggcCAGAGTTCTGGCGGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTcg            >  2:2199485/1‑100 (MQ=255)
                                                                                      cgcggcCAGAGTTCTGGCGGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTcg            >  2:1717912/1‑100 (MQ=255)
                                                                                      cgcggcCAGAGTTCTGGCGGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTcg            >  1:4339944/1‑100 (MQ=255)
                                                                                      cgcggcCAGAGTTCTGGCGGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTcg            >  2:1180485/1‑100 (MQ=255)
                                                                                      cgccgcCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTcg            >  2:1315555/1‑100 (MQ=255)
                                                                                      ccgatctaGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTcg            <  2:1495447/93‑1 (MQ=255)
                                                                                        ccgccAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTcgcg          >  2:1069574/1‑100 (MQ=255)
                                                                                        ccgccAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTcgcg          <  1:18237/100‑1 (MQ=255)
                                                                                        ccgccAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTcgcg          >  1:2440540/1‑100 (MQ=255)
                                                                                         cgctAGATCTCTGGCCGCGTCATCATGGTCATGGCCCGGCTGGTGTTAAGTTAGTGGAACAACTATGCCAGCGGCGGCGCGTGCCAAATGAAATTCTCGa         <  1:2625448/100‑1 (MQ=255)
                                                                                         cgccAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTcgcg          >  2:716926/1‑99 (MQ=255)
                                                                                         cgccAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGa         >  2:3587762/1‑100 (MQ=255)
                                                                                           ccAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGAt        >  1:952350/1‑99 (MQ=255)
                                                                                            cAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGAttt      >  2:3828261/1‑100 (MQ=255)
                                                                                            cAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGAttt      >  2:1589369/1‑100 (MQ=255)
                                                                                             aGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGAGAtcgg     >  1:1495441/1‑97 (MQ=255)
                                                                                                aCTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGcc  >  1:3609497/1‑100 (MQ=255)
                                                                                                |                                                                                                   
AATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCC  >  NC_000913/3202589‑3202784

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