Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,471,787 C→A 14.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,7870CA14.3% 912.3 / 105.4 280D24Y (GAT→TAT) slyDFKBP‑type peptidyl prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base C (114/125);  new base A (15/25);  total (129/151)
Fisher's exact test for biased strand distribution p-value = 3.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

Only 100 of 295 total aligned reads displayed.

TCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCGTGTACTCACTAGCTTAAAAAA  >  CP009273/3471639‑3471935
                                                                                                                                                    |                                                                                                                                                    
tcgtcgTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACtcatc                                                                                                                                                      >  2:4790450/1‑149 (MQ=255)
  gtcgtACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATcaa                                                                                                                                                    <  1:3001094/149‑1 (MQ=255)
        cTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAaca                                                                                                                                              <  2:2281039/149‑1 (MQ=255)
         tGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAacac                                                                                                                                             >  2:2839512/1‑149 (MQ=255)
         tGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAacac                                                                                                                                             <  2:188186/149‑1 (MQ=255)
           aCCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAacacac                                                                                                                                           <  2:4986207/149‑1 (MQ=255)
           aCCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAacacac                                                                                                                                           >  2:3998861/1‑149 (MQ=255)
             cGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCg                                                                                                                                         >  2:443346/1‑149 (MQ=255)
               tAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtc                                                                                                                                       >  2:4240993/1‑149 (MQ=255)
                aaGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtct                                                                                                                                      <  1:1029735/149‑1 (MQ=255)
                aaGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtct                                                                                                                                      <  2:2664965/149‑1 (MQ=255)
                 aGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtctt                                                                                                                                     >  2:4927962/1‑149 (MQ=255)
                 aGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtctt                                                                                                                                     >  1:920398/1‑149 (MQ=255)
                 aGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATAAACCAACACACCGtctt                                                                                                                                     >  1:2085559/1‑149 (MQ=255)
                  gCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtcttc                                                                                                                                    <  2:4761951/149‑1 (MQ=255)
                      cgtTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATAAACCAACACACCGTCTTCTgta                                                                                                                                >  1:4579855/1‑149 (MQ=255)
                        ttCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgtacg                                                                                                                              >  2:3917353/1‑149 (MQ=255)
                        ttCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgtacg                                                                                                                              >  1:74732/1‑149 (MQ=255)
                          cgcgCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgtacgta                                                                                                                            <  2:1224522/149‑1 (MQ=255)
                          cgcgCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgtacgta                                                                                                                            >  2:4556238/1‑149 (MQ=255)
                          cgcgCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgtacgta                                                                                                                            <  1:2788712/149‑1 (MQ=255)
                            cgcCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTAcc                                                                                                                          <  1:1682968/149‑1 (MQ=255)
                            cgcCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTAcc                                                                                                                          <  1:472829/149‑1 (MQ=255)
                               cAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGa                                                                                                                       >  1:3675208/1‑149 (MQ=255)
                                aaCAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGAt                                                                                                                      <  2:3038055/149‑1 (MQ=255)
                                     cGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATAAACCAACACACCGTCTTCTGTACGTACCTGATAGGcc                                                                                                                 <  2:1600088/149‑1 (MQ=255)
                                      gACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCa                                                                                                                <  1:2523358/149‑1 (MQ=255)
                                      gACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCa                                                                                                                >  1:3697536/1‑149 (MQ=255)
                                      gACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCa                                                                                                                <  2:3178791/149‑1 (MQ=255)
                                       aCATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAg                                                                                                               >  2:5187578/1‑149 (MQ=255)
                                        cATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAgg                                                                                                              <  2:3134415/149‑1 (MQ=255)
                                         aTCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCCGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGc                                                                                                             >  1:1469818/1‑149 (MQ=255)
                                         aTCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGc                                                                                                             <  2:1356554/149‑1 (MQ=255)
                                           cAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTg                                                                                                           <  1:1973352/149‑1 (MQ=255)
                                            aaaTTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATAAACCAACACACCGTCTTCTGTACGTACCGGATAGGCCAGGCTGa                                                                                                          >  1:2374187/1‑149 (MQ=255)
                                             aaTTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGac                                                                                                         >  2:714282/1‑149 (MQ=255)
                                             aaTTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATAAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGac                                                                                                         >  2:2546927/1‑149 (MQ=255)
                                               tttGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGacca                                                                                                       >  1:3941114/1‑149 (MQ=255)
                                                ttGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGaccac                                                                                                      <  1:4855201/149‑1 (MQ=255)
                                                  gTCGCCATCTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCTGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGAccacca                                                                                                    <  1:2923158/149‑1 (MQ=255)
                                                   tCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCCGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAg                                                                                                   >  2:4613764/1‑149 (MQ=255)
                                                   tCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAg                                                                                                   <  2:338360/149‑1 (MQ=255)
                                                   tCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATAAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAg                                                                                                   <  2:3803937/149‑1 (MQ=255)
                                                   tCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACACACCGGAGACACATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAg                                                                                                   >  1:375727/1‑149 (MQ=255)
                                                       cAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATAAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCt                                                                                               <  1:949203/149‑1 (MQ=255)
                                                         aCTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCGGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCttt                                                                                             >  2:2903487/1‑149 (MQ=255)
                                                            tCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGc                                                                                          <  1:1620551/149‑1 (MQ=255)
                                                            tCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGc                                                                                          >  1:1813620/1‑149 (MQ=255)
                                                             cATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCt                                                                                         <  2:373260/149‑1 (MQ=255)
                                                              aTGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTa                                                                                        <  1:4092823/149‑1 (MQ=255)
                                                              aTGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTa                                                                                        >  2:1940610/1‑149 (MQ=255)
                                                               tGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTAc                                                                                       <  1:4629219/149‑1 (MQ=255)
                                                               tGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTAc                                                                                       >  1:4812650/1‑149 (MQ=255)
                                                               tGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTAc                                                                                       <  1:2014077/149‑1 (MQ=255)
                                                                   cTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTc                                                                                   <  1:2124283/149‑1 (MQ=255)
                                                                    ttCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCa                                                                                  <  1:1883572/149‑1 (MQ=255)
                                                                     tCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATAAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCAt                                                                                 >  2:2560011/1‑149 (MQ=255)
                                                                         gcgcCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGata                                                                             <  1:1128203/149‑1 (MQ=255)
                                                                          cgcCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGatat                                                                            <  1:4038388/149‑1 (MQ=255)
                                                                             cGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATAtctc                                                                         <  2:618306/149‑1 (MQ=255)
                                                                              gTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTcc                                                                        <  2:563166/149‑1 (MQ=255)
                                                                              gTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTcc                                                                        <  2:175359/149‑1 (MQ=255)
                                                                              gTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTcc                                                                        <  2:2470002/149‑1 (MQ=255)
                                                                               tttCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCt                                                                       >  1:2589101/1‑149 (MQ=255)
                                                                                 tCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGa                                                                     <  1:1851023/149‑1 (MQ=255)
                                                                                   cAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGc                                                                   >  2:3586345/1‑149 (MQ=255)
                                                                                      gCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATg                                                                <  1:152960/149‑1 (MQ=255)
                                                                                       ccAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATgg                                                               <  1:5190675/149‑1 (MQ=255)
                                                                                        cAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATggg                                                              <  2:1674034/149‑1 (MQ=255)
                                                                                             aTCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGgaagaa                                                         <  2:2389110/149‑1 (MQ=255)
                                                                                             aTCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGgaagaa                                                         <  1:2054360/149‑1 (MQ=255)
                                                                                              tCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAAt                                                        <  2:2188663/149‑1 (MQ=255)
                                                                                                 gggAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATAAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGt                                                     >  1:2466321/1‑149 (MQ=255)
                                                                                                   gAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTgg                                                   <  1:3789188/149‑1 (MQ=255)
                                                                                                          gACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATAAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGAt                                            <  1:583633/149‑1 (MQ=255)
                                                                                                            ccATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTg                                          >  1:1446390/1‑149 (MQ=255)
                                                                                                            ccATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTg                                          >  1:1290121/1‑149 (MQ=255)
                                                                                                             cATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATAAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGt                                         <  1:692622/149‑1 (MQ=255)
                                                                                                             cATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATAAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGt                                         <  1:5059313/149‑1 (MQ=255)
                                                                                                              aTGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTa                                        <  1:4556238/149‑1 (MQ=255)
                                                                                                               tGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAg                                       <  1:1945387/149‑1 (MQ=255)
                                                                                                                   ggTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGa                                   >  2:720056/1‑149 (MQ=255)
                                                                                                                    gTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGaa                                  >  1:3540417/1‑149 (MQ=255)
                                                                                                                       gTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGACTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTc                               <  2:109030/149‑1 (MQ=255)
                                                                                                                             cggcgCACTCACCGGAGACTCATAAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGc                         <  2:1692826/149‑1 (MQ=255)
                                                                                                                              ggcgCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGcc                        >  1:1721972/1‑149 (MQ=255)
                                                                                                                                 gcACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCgtg                     >  2:2813224/1‑149 (MQ=255)
                                                                                                                                    cTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCGTGTAc                  >  1:2317910/1‑149 (MQ=255)
                                                                                                                                      cACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCGTGTACTc                <  1:3286994/149‑1 (MQ=255)
                                                                                                                                      cACCGGAGACTCATCAACCAACACACCGTCTTCTGGACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCGTGTACTc                <  2:4524261/149‑1 (MQ=255)
                                                                                                                                       aCCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCGTGTACTCa               >  1:1169673/1‑149 (MQ=255)
                                                                                                                                        ccGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCGTGTACTCAc              >  1:1018524/1‑149 (MQ=255)
                                                                                                                                         cGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCGTGTACTCACt             >  2:1507695/1‑149 (MQ=255)
                                                                                                                                          ggAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCGTGTACTCACTa            >  2:2753645/1‑149 (MQ=255)
                                                                                                                                          ggAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCGTGTACTCACTa            <  2:32426/149‑1 (MQ=255)
                                                                                                                                           gagaCTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCGTGTACTCACTAg           <  1:4799086/149‑1 (MQ=255)
                                                                                                                                           gagaCTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCGTGTACTCACTAg           >  1:4421863/1‑149 (MQ=255)
                                                                                                                                              aCTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCGTGTACTCACTAGCtt        >  1:2346950/1‑149 (MQ=255)
                                                                                                                                              aCTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCGTGTACTCACTAGCtt        <  2:3461325/149‑1 (MQ=255)
                                                                                                                                                tcatcaACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCGTGTACTCACTAGCTTaa      <  2:1826133/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACTTTCATGATATCTCCTGAGCATGGGAAGAATAGTGGCGCAGATTGTAGCGGAATTCTGCAGCCGTGTACTCACTAGCTTAAAAAA  >  CP009273/3471639‑3471935

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: ‑