Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,471,705 C→T 100% G51D (GGT→GAT)  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,7050CT96.0% 639.2 / 12.6 201G51D (GGT→GAT) slyDFKBP‑type peptidyl prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base C (4/4);  new base T (92/101);  total (96/105)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 201 total aligned reads displayed.

TTCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTA  >  CP009273/3471560‑3471849
                                                                                                                                                 |                                                                                                                                                
ttCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCt                                                                                                                                                <  2:3118633/148‑1 (MQ=255)
  cAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTc                                                                                                                                              <  1:1268458/148‑1 (MQ=255)
    gCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCa                                                                                                                                            >  1:634302/1‑148 (MQ=255)
    gCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCa                                                                                                                                            >  2:2629587/1‑148 (MQ=255)
    gCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCa                                                                                                                                            >  1:3551739/1‑148 (MQ=255)
         gAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGcc                                                                                                                                       >  2:1006000/1‑148 (MQ=255)
          aaaCGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCg                                                                                                                                      <  2:2874381/148‑1 (MQ=255)
          aaaCGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCg                                                                                                                                      <  2:1540344/148‑1 (MQ=255)
           aaCGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGt                                                                                                                                     <  2:1072466/148‑1 (MQ=255)
           aaCGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGt                                                                                                                                     <  1:718603/148‑1 (MQ=255)
                aTACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCa                                                                                                                                <  2:445987/148‑1 (MQ=255)
                aTACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCa                                                                                                                                <  2:2240665/148‑1 (MQ=255)
                 tacctaccTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAg                                                                                                                               >  1:2280507/1‑148 (MQ=255)
                   cctacctGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGc                                                                                                                             >  1:2731255/1‑148 (MQ=255)
                      acctGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCag                                                                                                                          >  1:2460390/1‑148 (MQ=255)
                       cctGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCaga                                                                                                                         >  1:3503446/1‑148 (MQ=255)
                       cctGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCaga                                                                                                                         >  2:276538/1‑148 (MQ=255)
                        ctGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCagag                                                                                                                        >  1:1063279/1‑148 (MQ=255)
                        ctGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTAACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCagag                                                                                                                        <  1:2585633/148‑1 (MQ=255)
                         tGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCCgaga                                                                                                                       >  1:1567782/1‑148 (MQ=255)
                         tGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAgaga                                                                                                                       >  2:471559/1‑148 (MQ=255)
                         tGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAgaga                                                                                                                       >  1:2777459/1‑148 (MQ=255)
                           cAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATc                                                                                                                     >  1:2534664/1‑148 (MQ=255)
                            aGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCa                                                                                                                    <  1:2543233/148‑1 (MQ=255)
                             gTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACGTCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAg                                                                                                                   >  1:1825532/1‑148 (MQ=255)
                             gTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAg                                                                                                                   <  1:851289/148‑1 (MQ=255)
                               tcatcaACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAggg                                                                                                                 >  2:1386379/1‑148 (MQ=255)
                                catcaACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCTTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGa                                                                                                                <  2:2952293/148‑1 (MQ=255)
                                 atcaACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGaa                                                                                                               <  2:1341594/148‑1 (MQ=255)
                                 atcaACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGaa                                                                                                               >  1:351647/1‑148 (MQ=255)
                                 atcaACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGaa                                                                                                               <  1:336966/148‑1 (MQ=255)
                                 atcaACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGaa                                                                                                               <  1:336552/148‑1 (MQ=255)
                                 atcaACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGaa                                                                                                               <  2:1686237/148‑1 (MQ=255)
                                    aaCGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCg                                                                                                            <  1:2494074/148‑1 (MQ=255)
                                    aaCGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCg                                                                                                            <  2:698603/148‑1 (MQ=255)
                                            tAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCTACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATg                                                                                                    <  2:1847729/148‑1 (MQ=255)
                                            tAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATg                                                                                                    >  1:179561/1‑148 (MQ=255)
                                            tAAATACGTCTTTAGGAACACGGTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCTACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATg                                                                                                    <  2:1859243/148‑1 (MQ=255)
                                             aaaTACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGACCCGTGACCATGc                                                                                                   >  2:968017/1‑148 (MQ=255)
                                                 aCGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGt                                                                                               >  1:3158559/1‑148 (MQ=255)
                                                 aCGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGt                                                                                               <  2:3416650/148‑1 (MQ=255)
                                                 aCGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGt                                                                                               <  1:992968/148‑1 (MQ=255)
                                                    tCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGt                                                                                            <  2:2385249/148‑1 (MQ=255)
                                                    tCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGt                                                                                            <  2:745881/148‑1 (MQ=255)
                                                    tCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGt                                                                                            >  1:1464337/1‑148 (MQ=255)
                                                         aGGAACACGTTGCACCAGGTGTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAgc                                                                                       <  2:2123966/148‑1 (MQ=255)
                                                          ggAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAgcg                                                                                      >  1:3079616/1‑148 (MQ=255)
                                                           gAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAgcgg                                                                                     >  1:2905203/1‑148 (MQ=255)
                                                              cacGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGAACAGGGAACCGTGACCATGCAGGTAGTCCAGCGgcgc                                                                                  >  1:2243267/1‑148 (MQ=255)
                                                               acGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCa                                                                                 <  2:401850/148‑1 (MQ=255)
                                                                cGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCAc                                                                                >  2:3053199/1‑148 (MQ=255)
                                                                   tGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCa                                                                             >  2:2412480/1‑148 (MQ=255)
                                                                    gCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCAc                                                                            <  2:1460081/148‑1 (MQ=255)
                                                                       ccAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCgg                                                                         <  2:2524938/148‑1 (MQ=255)
                                                                        cAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGga                                                                        <  2:2089831/148‑1 (MQ=255)
                                                                          ggTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCCGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGgaga                                                                      >  1:364434/1‑148 (MQ=255)
                                                                          ggTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGgaga                                                                      >  2:576205/1‑148 (MQ=255)
                                                                          ggTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGgaga                                                                      <  2:923772/148‑1 (MQ=255)
                                                                            ttttCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACt                                                                    >  1:664696/1‑148 (MQ=255)
                                                                             tttCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACtc                                                                   >  1:2029902/1‑148 (MQ=255)
                                                                                  tcgtACTGACCGTAAGCGTCGTTCTCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGGCGTTTCCAGGCCAGAGATCAGGGAACAGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCTTcaa                                                              <  2:1016547/148‑1 (MQ=255)
                                                                                  tcgtACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATcaa                                                              <  2:3017753/148‑1 (MQ=255)
                                                                                  tcgtACTGACCGGAAGCGTCGTTCGCGCCAACAGCGACATCAAATGTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATcaa                                                              <  2:2498132/148‑1 (MQ=255)
                                                                                       cTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCCGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATcaaccaac                                                         >  1:804186/1‑148 (MQ=255)
                                                                                       cTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATcaaccaac                                                         <  2:1762351/148‑1 (MQ=255)
                                                                                       cTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATcaaccaac                                                         >  1:3085932/1‑148 (MQ=255)
                                                                                          aCCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAacaca                                                      <  1:2110977/148‑1 (MQ=255)
                                                                                           ccGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAacacac                                                     >  2:3505189/1‑148 (MQ=255)
                                                                                              tAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGt                                                  <  2:1107257/148‑1 (MQ=255)
                                                                                               aaGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtc                                                 >  1:112935/1‑148 (MQ=255)
                                                                                               aaGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtc                                                 >  1:3509327/1‑148 (MQ=255)
                                                                                                aGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCNGGGAACCGTGACCATGCANGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtct                                                <  2:3318795/148‑1 (MQ=255)
                                                                                                 gCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtctt                                               <  2:1931762/148‑1 (MQ=255)
                                                                                                        tCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgt                                           <  1:939322/145‑1 (MQ=255)
                                                                                                           cgcCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgtacgtac                                     <  2:2460390/148‑1 (MQ=255)
                                                                                                            gcCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTAcc                                    <  1:738375/148‑1 (MQ=255)
                                                                                                               aaCAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGa                                 <  1:1043590/148‑1 (MQ=255)
                                                                                                                aCAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCCGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGAt                                >  1:3292656/1‑148 (MQ=255)
                                                                                                                aCAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGAt                                >  1:3379937/1‑148 (MQ=255)
                                                                                                                aCAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGAt                                <  2:2797354/148‑1 (MQ=255)
                                                                                                                aCAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGAt                                >  2:622357/1‑148 (MQ=255)
                                                                                                                aCAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGAt                                >  2:1196475/1‑148 (MQ=255)
                                                                                                                  aGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAg                              >  2:2862475/1‑148 (MQ=255)
                                                                                                                  aGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAg                              >  2:2674078/1‑148 (MQ=255)
                                                                                                                     gACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGcc                           <  2:789935/148‑1 (MQ=255)
                                                                                                                     gACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGcc                           <  1:1382482/148‑1 (MQ=255)
                                                                                                                      aCATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCa                          <  2:2905203/148‑1 (MQ=255)
                                                                                                                        aTCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAgg                        <  1:3122195/148‑1 (MQ=255)
                                                                                                                           aaaTTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTg                     <  2:2777420/148‑1 (MQ=255)
                                                                                                                            aaTTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGa                    <  2:3503303/148‑1 (MQ=255)
                                                                                                                             aTTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGac                   >  2:2992360/1‑148 (MQ=255)
                                                                                                                             aTTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGac                   >  2:2829034/1‑148 (MQ=255)
                                                                                                                             aTTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTAAGTACCTGATAGGCCAGGCTGac                   >  1:1217037/1‑148 (MQ=255)
                                                                                                                               ttGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGacca                 >  2:725756/1‑148 (MQ=255)
                                                                                                                                     ccAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGt           <  2:1063183/148‑1 (MQ=255)
                                                                                                                                      cAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTc          <  1:276538/148‑1 (MQ=255)
                                                                                                                                      cAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTc          <  1:2986915/148‑1 (MQ=255)
                                                                                                                                        aCTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCtt        <  1:853942/148‑1 (MQ=255)
                                                                                                                                         cTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCttt       <  1:2156552/148‑1 (MQ=255)
                                                                                                                                            cATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGc    <  1:50096/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTA  >  CP009273/3471560‑3471849

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: