Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,471,787 C→A 7.3% D24Y (GAT→TAT)  slyD ← FKBP‑type peptidyl prolyl cis‑trans isomerase (rotamase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,471,7870CA 7.3% 929.7 / 37.8 245D24Y (GAT→TAT) slyDFKBP‑type peptidyl prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base C (98/129);  new base A (10/8);  total (108/137)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 255 total aligned reads displayed.

TCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCGTGTACTCACTAGCTTAAAA  >  CP009273/3471639‑3471933
                                                                                                                                                    |                                                                                                                                                  
   tcgtACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATcaac                                                                                                                                                 >  2:2892523/1‑149 (MQ=255)
         tGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAacac                                                                                                                                           >  1:553195/1‑149 (MQ=255)
           aCCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAacacac                                                                                                                                         <  2:558697/149‑1 (MQ=255)
                aaGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtct                                                                                                                                    <  1:4477333/149‑1 (MQ=255)
                aaGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATAAACCAACACACCGtct                                                                                                                                    <  2:4620281/149‑1 (MQ=255)
                 aGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtctt                                                                                                                                   >  2:4501961/1‑149 (MQ=255)
                   cgtcgtTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtcttct                                                                                                                                 <  1:2249607/149‑1 (MQ=255)
                     tcgtTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgt                                                                                                                               <  1:1349442/149‑1 (MQ=255)
                     tcgtTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgt                                                                                                                               <  2:3149528/149‑1 (MQ=255)
                      cgtTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgta                                                                                                                              <  1:2411195/149‑1 (MQ=255)
                         tCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgtacgt                                                                                                                           >  2:1737835/1‑149 (MQ=255)
                          cgcgCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGATCCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgtacgta                                                                                                                          <  2:3848648/149‑1 (MQ=255)
                           gcgcCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgtacgtac                                                                                                                         <  1:3659986/149‑1 (MQ=255)
                               cAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGa                                                                                                                     >  1:3503412/1‑149 (MQ=255)
                                aaCAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGAt                                                                                                                    <  2:4290980/149‑1 (MQ=255)
                                aaCAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTACAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGAGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTATGTACGTACCTGAt                                                                                                                    >  1:151178/1‑149 (MQ=255)
                                  cAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAg                                                                                                                  <  1:3278482/149‑1 (MQ=255)
                                  cAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGGGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAg                                                                                                                  <  2:25616/149‑1 (MQ=255)
                                    gCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGc                                                                                                                <  1:2977551/149‑1 (MQ=255)
                                    gCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGc                                                                                                                <  1:2523896/149‑1 (MQ=255)
                                     cGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGcc                                                                                                               <  1:2805007/149‑1 (MQ=255)
                                       aCATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAg                                                                                                             <  2:3941834/149‑1 (MQ=255)
                                        cATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAgg                                                                                                            <  1:3354367/149‑1 (MQ=255)
                                         aTCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGc                                                                                                           >  2:893781/1‑149 (MQ=255)
                                         aTCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGc                                                                                                           <  1:3922336/149‑1 (MQ=255)
                                         aTCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGc                                                                                                           <  1:2535397/149‑1 (MQ=255)
                                          tCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCt                                                                                                          >  1:3757053/1‑149 (MQ=255)
                                                     gCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGt                                                                                               >  1:17813/1‑149 (MQ=255)
                                                     gCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGt                                                                                               >  2:4230998/1‑149 (MQ=255)
                                                     gCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGt                                                                                               >  1:2316443/1‑149 (MQ=255)
                                                      ccAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTc                                                                                              >  2:1242498/1‑149 (MQ=255)
                                                       cAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCt                                                                                             <  1:1984816/149‑1 (MQ=255)
                                                        aaCTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCtt                                                                                            >  2:1976742/1‑149 (MQ=255)
                                                         aCTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATAAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCttt                                                                                           >  2:3255706/1‑149 (MQ=255)
                                                          cTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCtttt                                                                                          <  1:4319372/149‑1 (MQ=255)
                                                            tCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGc                                                                                        >  1:1101409/1‑149 (MQ=255)
                                                              aTGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTa                                                                                      >  1:3016544/1‑149 (MQ=255)
                                                               tGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATAAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCCCCAGGTCTTTTGCTAc                                                                                     >  2:1247710/1‑149 (MQ=255)
                                                                 aCCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACtt                                                                                   <  1:1929521/149‑1 (MQ=255)
                                                                 aCCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACtt                                                                                   >  1:2378229/1‑149 (MQ=255)
                                                                   cTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTc                                                                                 <  1:4309790/149‑1 (MQ=255)
                                                                    ttCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCa                                                                                <  1:3096984/149‑1 (MQ=255)
                                                                    ttCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCa                                                                                <  1:2242383/149‑1 (MQ=255)
                                                                    ttCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCa                                                                                <  1:1639377/149‑1 (MQ=255)
                                                                     tCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACCCCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCAt                                                                               <  2:2654432/149‑1 (MQ=255)
                                                                             cGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATAtctc                                                                       >  2:1649048/1‑149 (MQ=255)
                                                                             cGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATAtctc                                                                       <  2:2837781/149‑1 (MQ=255)
                                                                              gTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTcc                                                                      >  2:741317/1‑149 (MQ=255)
                                                                               tttCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCt                                                                     <  2:1199599/149‑1 (MQ=255)
                                                                                     ggCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCAt                                                               >  1:1696248/1‑149 (MQ=255)
                                                                                     ggCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCAt                                                               <  2:667692/149‑1 (MQ=255)
                                                                                      gCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATg                                                              <  2:389535/149‑1 (MQ=255)
                                                                                       ccAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATgg                                                             >  1:1043256/1‑149 (MQ=255)
                                                                                        cAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATggg                                                            >  1:1923604/1‑149 (MQ=255)
                                                                                        cAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATggg                                                            >  1:3327985/1‑149 (MQ=255)
                                                                                        cAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATggg                                                            >  1:3327774/1‑149 (MQ=255)
                                                                                         agagATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGga                                                           >  1:402578/1‑149 (MQ=255)
                                                                                         agagATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGga                                                           >  2:3452280/1‑149 (MQ=255)
                                                                                         agagATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATAAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGga                                                           >  2:1547016/1‑149 (MQ=255)
                                                                                           agaTCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGgaag                                                         >  1:788267/1‑149 (MQ=255)
                                                                                            gaTCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGgaaga                                                        <  2:2684413/149‑1 (MQ=255)
                                                                                                aGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAg                                                    >  2:33876/1‑149 (MQ=255)
                                                                                                 gggAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGt                                                   >  2:357985/1‑149 (MQ=255)
                                                                                                 gggAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGt                                                   <  1:3866479/149‑1 (MQ=255)
                                                                                                   gAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTgg                                                 <  2:1052924/149‑1 (MQ=255)
                                                                                                   gAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTgg                                                 <  1:632922/149‑1 (MQ=255)
                                                                                                         tGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGa                                           >  1:942876/1‑149 (MQ=255)
                                                                                                         tGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGa                                           >  2:4342667/1‑149 (MQ=255)
                                                                                                         tGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATAAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGa                                           >  1:3070880/1‑149 (MQ=255)
                                                                                                            ccATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTg                                        >  1:896584/1‑149 (MQ=255)
                                                                                                             cATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGt                                       <  2:4194215/149‑1 (MQ=255)
                                                                                                              aTGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTa                                      <  2:296319/149‑1 (MQ=255)
                                                                                                              aTGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTa                                      <  2:2202984/149‑1 (MQ=255)
                                                                                                               tGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAg                                     <  1:1609759/149‑1 (MQ=255)
                                                                                                                gCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGc                                    <  2:3632582/149‑1 (MQ=255)
                                                                                                                    gTAGTCCAGCGGCGCACTCACCGGAGACTCATAAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGaa                                >  2:4011279/1‑149 (MQ=255)
                                                                                                                      aGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAAtt                              <  1:1185735/149‑1 (MQ=255)
                                                                                                                       gTCCAGCGGCGCACTCAGGGGTGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTc                             <  2:2423292/149‑1 (MQ=255)
                                                                                                                       gTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTc                             <  1:4237966/149‑1 (MQ=255)
                                                                                                                       gTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTc                             <  2:1581151/149‑1 (MQ=255)
                                                                                                                       gTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTc                             <  2:1106192/149‑1 (MQ=255)
                                                                                                                       gTCCAGCGGCGCACTCACCGGAGACTCATAAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTc                             >  1:2280836/1‑149 (MQ=255)
                                                                                                                        tCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCt                            >  1:2189036/1‑149 (MQ=255)
                                                                                                                          cAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGc                          <  2:2803267/149‑1 (MQ=255)
                                                                                                                            gcggcgCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAg                        >  2:2262972/1‑149 (MQ=255)
                                                                                                                            gcggcgCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAg                        <  1:1606636/149‑1 (MQ=255)
                                                                                                                              ggcgCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGcc                      <  1:2358448/149‑1 (MQ=255)
                                                                                                                               gcgcACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCg                     <  1:4575059/149‑1 (MQ=255)
                                                                                                                                  cACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCgtgt                  >  2:4371621/1‑149 (MQ=255)
                                                                                                                                  cACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCgtgt                  <  2:4601564/149‑1 (MQ=255)
                                                                                                                                    cTCACCGGAGACTCATAAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCGTGTAc                >  1:4523370/1‑149 (MQ=255)
                                                                                                                                      cACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCGTGTACTc              >  1:3514838/1‑149 (MQ=255)
                                                                                                                                       aCCGGAGACTCACCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCGTGTACTCa             <  1:4446484/149‑1 (MQ=255)
                                                                                                                                        ccGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCGTGTACTCAc            <  1:1423177/149‑1 (MQ=255)
                                                                                                                                          ggAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCGTGTACTCACTa          <  1:2396963/149‑1 (MQ=255)
                                                                                                                                           gagaCTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCGTGTACTCACTAg         <  2:1490698/149‑1 (MQ=255)
                                                                                                                                             gaCTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCGTGTACTCACTAGCt       <  1:36030/149‑1 (MQ=255)
                                                                                                                                               cTCAGCAACCAACACACCGTCGTCTGTACGTACCTGAGAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGACTAGTGGCGCAGATTGGAGCGGAATTCTGCAGCCGTGTACTCACGAGCTTa     >  2:3077460/1‑149 (MQ=255)
                                                                                                                                                tcatcaACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCGTGTACTCACTAGCTTaa    >  1:1261456/1‑149 (MQ=255)
                                                                                                                                                  atcaACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCGTGTACTCACTAGCTTaaaa  >  2:129355/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCGTGTACTCACTAGCTTAAAA  >  CP009273/3471639‑3471933

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