Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,202,678 C→G 6.1% P263R (CCG→CGG)  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,6780CG 6.1% 345.2 / 11.2 98P263R (CCG→CGG) ccafused tRNA nucleotidyl transferase/2'3'‑cyclic phosphodiesterase/2'nucleotidase and phosphatase
Reads supporting (aligned to +/- strand):  ref base C (44/48);  new base G (3/3);  total (47/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

GCCCCTGCCAAGTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTC  >  NC_000913/3202542‑3202823
                                                                                                                                        |                                                                                                                                                 
gCCCCTGCCAAGTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTg                                                                                                                                       <  3:1009441/149‑1 (MQ=255)
   ccTGCCAAGTGGCATCCGGAAATCGATACGGGTATTCATACCTTATTGACGCTCTCAATGGCGGCGATTCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCGGCCAGAGTTCTGGcg                                                                                                                                    >  4:842333/1‑148 (MQ=255)
    cTGCCAAGTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCcg                                                                                                                                   >  2:959421/1‑149 (MQ=255)
     tGCCAAGTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCGGCCAGAGTTCTggcggc                                                                                                                                  <  3:842333/149‑1 (MQ=255)
     tGCCAAGTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCcgc                                                                                                                                  >  3:1352858/1‑149 (MQ=255)
      gCCAAGTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCCCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCGGCCAGAGTTCTGgcggcg                                                                                                                                 >  4:366239/1‑149 (MQ=255)
         aaGTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCCCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCGGCCAGAGTTCTGGCGGCGtca                                                                                                                              <  3:366239/149‑1 (MQ=255)
           gTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGtcatc                                                                                                                            <  1:1233307/149‑1 (MQ=255)
              gCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatg                                                                                                                         <  4:147028/149‑1 (MQ=255)
               cATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcacgg                                                                                                                        >  4:183234/1‑149 (MQ=255)
                 tCCGGAAATCGTTACGGGGATTCATACCTTAATGACGCTCTCGATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtc                                                                                                                      >  2:330008/1‑149 (MQ=255)
                 tCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtc                                                                                                                      >  4:850925/1‑149 (MQ=255)
                 tCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtc                                                                                                                      <  4:1352858/149‑1 (MQ=255)
                 tCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtc                                                                                                                      <  1:959421/149‑1 (MQ=255)
                             tACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCACGGTCATGGCCcggcgg                                                                                                          <  3:183234/149‑1 (MQ=255)
                                gggTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTcatcat                                                                                                                          <  3:110583/130‑1 (MQ=255)
                                gggTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTcatcat                                                                                                                          >  4:110583/1‑130 (MQ=255)
                                gggTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgtg                                                                                                       >  1:746153/1‑149 (MQ=255)
                                  gTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCcggc                                                                                                            >  4:874360/1‑142 (MQ=255)
                                  gTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCcggc                                                                                                            <  3:874360/142‑1 (MQ=255)
                                  gTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGtt                                                                                                     >  1:1328577/1‑149 (MQ=255)
                                    aTTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTaa                                                                                                   <  2:746153/149‑1 (MQ=255)
                                      tCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGt                                                                                                 >  2:808167/1‑149 (MQ=255)
                                      tCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGt                                                                                                 <  1:330008/149‑1 (MQ=255)
                                        aTACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCGCGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgt                                                                                                        <  2:987633/140‑1 (MQ=255)
                                        aTACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCGCGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgt                                                                                                        >  1:987633/1‑140 (MQ=255)
                                          aCCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGt                                                                                             <  1:808167/149‑1 (MQ=255)
                                           ccTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTg                                                                                            <  3:850925/149‑1 (MQ=255)
                                            cTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGt                                                                                             <  2:870831/147‑1 (MQ=255)
                                            cTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGt                                                                                             >  1:870831/1‑147 (MQ=255)
                                            cTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTgg                                                                                           >  1:842512/1‑149 (MQ=255)
                                             ttAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGa                                                                                          >  3:623823/1‑149 (MQ=255)
                                               aaTGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGtt                                                                                                     <  2:1339729/136‑1 (MQ=255)
                                               aaTGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGtt                                                                                                     >  1:1339729/1‑136 (MQ=255)
                                                 tGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGaacaa                                                                                      <  4:17056/149‑1 (MQ=255)
                                                   aCGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgt                                                                                                        >  3:503229/1‑129 (MQ=255)
                                                   aCGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgt                                                                                                        <  4:503229/129‑1 (MQ=255)
                                                     gCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGt                                                                                                 <  4:501148/134‑1 (MQ=255)
                                                     gCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGaacaa                                                                                      <  2:1022917/145‑1 (MQ=255)
                                                     gCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGaacaa                                                                                      >  1:1022917/1‑145 (MQ=255)
                                                     gCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGCTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGt                                                                                                 >  3:501148/1‑134 (MQ=255)
                                                       tctcAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGt                                                                                                 <  3:829238/132‑1 (MQ=255)
                                                       tctcAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGt                                                                                                 >  4:829238/1‑132 (MQ=255)
                                                         tcAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTaa                                                                                                   <  1:1076026/128‑1 (MQ=255)
                                                         tcAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTaa                                                                                                   >  2:1076026/1‑128 (MQ=255)
                                                         tcAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGATGACGCGGCCAGAGTTCTggcggc                                                                                                                                  <  2:547127/97‑1 (MQ=38)
                                                         tcAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGATGACGCGGCCAGAGTTCTggcggc                                                                                                                                  >  1:547127/1‑97 (MQ=38)
                                                              ggcggcGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTc                                                                         <  2:1328577/149‑1 (MQ=255)
                                                                 ggcgATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGaa                                                                                                                                            >  2:16774/1‑79 (MQ=255)
                                                                 ggcgATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGaa                                                                                                                                            <  1:16774/79‑1 (MQ=255)
                                                                   cgATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTa                                                                                                    <  3:925489/117‑1 (MQ=255)
                                                                   cgATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTa                                                                                                    >  4:925489/1‑117 (MQ=255)
                                                                       gCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGc                                                                >  2:1208422/1‑149 (MQ=255)
                                                                        cTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCcggcg                                                                                                           <  1:939890/105‑1 (MQ=255)
                                                                        cTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCcggcg                                                                                                           >  2:939890/1‑105 (MQ=255)
                                                                         tGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGaa                                                                                                                                            <  3:445181/71‑1 (MQ=255)
                                                                         tGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGaa                                                                                                                                            >  4:445181/1‑71 (MQ=255)
                                                                             tCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCCCGTCAACATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATg                                                          <  4:898413/149‑1 (MQ=255)
                                                                               cGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGaa                                                        >  4:850809/1‑149 (MQ=255)
                                                                                gCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGaa                                                                                         <  3:240946/115‑1 (MQ=255)
                                                                                gCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGaa                                                                                         >  4:240946/1‑115 (MQ=255)
                                                                                      cGATGTCCGTTTCGCAACTTTATGCCACGATCCCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTg                                                                                            <  3:297234/106‑1 (MQ=255)
                                                                                      cGATGTCCGTTTCGCAACTTTATGCCACGATCCCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTg                                                                                            >  4:297234/1‑106 (MQ=255)
                                                                                       gATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGa                                                <  1:1208422/149‑1 (MQ=255)
                                                                                          gTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGt                                                                                             >  1:682420/1‑101 (MQ=255)
                                                                                          gTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGt                                                                                             <  2:682420/101‑1 (MQ=255)
                                                                                          gTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGAttt                                             >  2:1406862/1‑149 (MQ=255)
                                                                                          gTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGAttt                                             >  2:841039/1‑149 (MQ=255)
                                                                                           tCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTcgc                                                  >  4:183342/1‑143 (MQ=255)
                                                                                           tCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTcgc                                                  <  3:183342/143‑1 (MQ=255)
                                                                                           tCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTa                                            >  3:1137292/1‑149 (MQ=255)
                                                                                            ccGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGaac                                                                                        <  3:40487/104‑1 (MQ=255)
                                                                                            ccGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGaac                                                                                        >  4:40487/1‑104 (MQ=255)
                                                                                            ccGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAg                                           <  2:842512/149‑1 (MQ=255)
                                                                                              gTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGc                                                                                <  3:1272422/110‑1 (MQ=255)
                                                                                              gTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGc                                                                                >  4:1272422/1‑110 (MQ=255)
                                                                                              gTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGcc                                         <  4:623823/149‑1 (MQ=255)
                                                                                                   gCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACt                                    >  3:661848/1‑149 (MQ=255)
                                                                                                     aaCTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCaaa                                                            >  2:759131/1‑123 (MQ=255)
                                                                                                     aaCTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCATGCCaaa                                                            <  1:759131/123‑1 (MQ=255)
                                                                                                       cTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGtt                                                                                                     <  3:1245488/80‑1 (MQ=255)
                                                                                                       cTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGtt                                                                                                     >  4:1245488/1‑80 (MQ=255)
                                                                                                          tATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCt                             <  4:1137292/149‑1 (MQ=255)
                                                                                                              ccACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGAc                                     >  4:1184107/1‑137 (MQ=255)
                                                                                                              ccACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGAc                                     <  3:1184107/137‑1 (MQ=255)
                                                                                                              ccACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGt                         <  3:850809/149‑1 (MQ=255)
                                                                                                               cACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACtggtgg                               >  1:913197/1‑142 (MQ=255)
                                                                                                               cACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACtggtgg                               <  2:913197/142‑1 (MQ=255)
                                                                                                                 cGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGTCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGt                                                                                             <  2:217758/78‑1 (MQ=255)
                                                                                                                 cGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGTCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGt                                                                                             >  1:217758/1‑78 (MQ=255)
                                                                                                                     ctcGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTAt                                                                                  >  3:597502/1‑85 (MQ=255)
                                                                                                                     ctcGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTAt                                                                                  <  4:597502/85‑1 (MQ=255)
                                                                                                                      tcGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACt                                                                                    >  4:598099/1‑82 (MQ=255)
                                                                                                                      tcGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACt                                                                                    <  3:598099/82‑1 (MQ=255)
                                                                                                                       cGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGAtc                <  1:1406862/149‑1 (MQ=255)
                                                                                                                       cGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGAtc                <  1:841039/149‑1 (MQ=255)
                                                                                                                                   gACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATc           >  1:289201/1‑142 (MQ=255)
                                                                                                                                   gACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATc           <  2:289201/142‑1 (MQ=255)
                                                                                                                                     cgccgcCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTc  <  4:661848/149‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                 
GCCCCTGCCAAGTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTC  >  NC_000913/3202542‑3202823

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