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,7050CT95.8% 762.7 / 18.6 235G51D (GGT→GAT) slyDFKBP‑type peptidyl prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base C (5/5);  new base T (121/104);  total (126/109)
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.91e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 235 total aligned reads displayed.

TTCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTA  >  CP009273/3471560‑3471849
                                                                                                                                                 |                                                                                                                                                
ttCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCt                                                                                                                                                >  2:899940/1‑148 (MQ=255)
 tCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCtt                                                                                                                                               <  1:89293/148‑1 (MQ=255)
   aGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTcc                                                                                                                                             <  1:612293/148‑1 (MQ=255)
        ggAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAgcgc                                                                                                                                        >  2:1832977/1‑148 (MQ=255)
        ggAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACAAAAAATTTGTCGCCACATTCATGATCTGCAAgcgc                                                                                                                                        >  1:126009/1‑148 (MQ=255)
          aaaCGCATACCTACCTGCAGTTCATCAACGCCCATAGATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCg                                                                                                                                      >  2:2889512/1‑148 (MQ=255)
          aaaCGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCg                                                                                                                                      >  1:2324240/1‑148 (MQ=255)
           aaCGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGt                                                                                                                                     >  1:684724/1‑148 (MQ=255)
           aaCGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGt                                                                                                                                     <  2:1720970/148‑1 (MQ=255)
            aCGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGtt                                                                                                                                    >  2:1265353/1‑148 (MQ=255)
             cGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGttt                                                                                                                                   <  1:1362517/148‑1 (MQ=255)
             cGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGAACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACAACAAATTTGTAGCCAACTTCAAGATCTTCCAGCGCCGttt                                                                                                                                   >  1:2652654/1‑148 (MQ=255)
                aTACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCa                                                                                                                                <  1:2291513/148‑1 (MQ=255)
                  acctacctGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAgg                                                                                                                              >  1:681634/1‑148 (MQ=255)
                   cctacctGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGc                                                                                                                             >  1:768635/1‑148 (MQ=255)
                     tacctGCAGTTCATCAACGCCCATAAATAAGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGAGCCGTTTCCAGGCCa                                                                                                                           >  2:2449922/1‑148 (MQ=255)
                      acctGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCag                                                                                                                          <  2:1056850/148‑1 (MQ=255)
                       cctGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCaga                                                                                                                         >  1:1971317/1‑148 (MQ=255)
                        ctGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCagag                                                                                                                        >  2:2668934/1‑148 (MQ=255)
                        ctGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCagag                                                                                                                        <  2:3174683/148‑1 (MQ=255)
                         tGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAgaga                                                                                                                       <  1:200364/148‑1 (MQ=255)
                           cAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATc                                                                                                                     <  2:2397956/148‑1 (MQ=255)
                           cAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATc                                                                                                                     >  2:664910/1‑148 (MQ=255)
                            aGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCa                                                                                                                    <  2:416406/148‑1 (MQ=255)
                             gTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAg                                                                                                                   <  2:574117/148‑1 (MQ=255)
                              ttaatcaACGCCCATAAATACGTCGTTAGGAAAACGTTGCACCAGGTGTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACAACAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGCCCAGAGAACAgg                                                                                                                  >  2:3431288/4‑148 (MQ=255)
                              ttCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAgg                                                                                                                  >  1:3350651/1‑148 (MQ=255)
                               tcatcaACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAggg                                                                                                                 >  1:119526/1‑148 (MQ=255)
                               tcatcaACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAggg                                                                                                                 >  2:3134163/1‑148 (MQ=255)
                                 atcaACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGaa                                                                                                               <  2:251181/148‑1 (MQ=255)
                                           aTAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCAt                                                                                                     >  1:542698/1‑148 (MQ=255)
                                           aTAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCACAGATCAGGGAACCGTGACCAt                                                                                                     >  1:3441144/1‑148 (MQ=255)
                                           aTAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTAAACAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCAt                                                                                                     >  2:293967/1‑148 (MQ=255)
                                            tAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATg                                                                                                    >  1:454050/1‑148 (MQ=255)
                                            tAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTACAGGCCACAGATCAGGGAACAGTGACCAtt                                                                                                    >  1:1681152/1‑147 (MQ=255)
                                             aaaTACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCACATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGc                                                                                                   <  2:2741216/148‑1 (MQ=255)
                                              aaTACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCa                                                                                                  >  2:1605246/1‑148 (MQ=255)
                                               aTACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAg                                                                                                 >  1:2175788/1‑148 (MQ=255)
                                                 aCGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGt                                                                                               >  1:825841/1‑148 (MQ=255)
                                                 aCGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGt                                                                                               >  2:384388/1‑148 (MQ=255)
                                                    tCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGt                                                                                            <  1:717054/148‑1 (MQ=255)
                                                       ttAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCa                                                                                         <  2:2634854/148‑1 (MQ=255)
                                                          ggAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAgcg                                                                                      >  2:684830/1‑148 (MQ=255)
                                                           gAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAgcgg                                                                                     <  1:1389250/148‑1 (MQ=255)
                                                           gAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAgcgg                                                                                     >  2:2769740/1‑148 (MQ=255)
                                                             acacGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAgcggcg                                                                                   >  2:2139711/1‑148 (MQ=255)
                                                                cGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCAc                                                                                >  2:772698/1‑148 (MQ=255)
                                                                       ccAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCgg                                                                         <  1:1492467/148‑1 (MQ=255)
                                                                       ccAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCgg                                                                         <  1:3082105/148‑1 (MQ=255)
                                                                       ccAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCgg                                                                         <  1:3019645/148‑1 (MQ=255)
                                                                         aGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGgag                                                                       >  2:2775240/1‑148 (MQ=255)
                                                                         aGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGgag                                                                       >  2:3115208/1‑148 (MQ=255)
                                                                         aGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTGGGCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAATCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGgag                                                                       <  1:1572541/148‑1 (MQ=255)
                                                                          ggTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGgaga                                                                      <  1:1433291/148‑1 (MQ=255)
                                                                            ttttCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACt                                                                    >  2:1138848/1‑148 (MQ=255)
                                                                             tttCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACtc                                                                   <  2:1556600/148‑1 (MQ=255)
                                                                              ttCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACtca                                                                  <  1:2057617/148‑1 (MQ=255)
                                                                               tcgtcgTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACtcat                                                                 <  2:1265941/148‑1 (MQ=255)
                                                                               tcgtcgTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACtcat                                                                 <  2:1401311/148‑1 (MQ=255)
                                                                                  tcgtACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATcaa                                                              >  1:3420966/1‑148 (MQ=255)
                                                                                      aCTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATcaaccaa                                                          >  2:1937397/1‑148 (MQ=255)
                                                                                      aCTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATcaaccaa                                                          >  1:797453/1‑148 (MQ=255)
                                                                                             gTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCg                                                   <  2:3407922/148‑1 (MQ=255)
                                                                                             gTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCg                                                   <  2:2116334/148‑1 (MQ=255)
                                                                                              tAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGt                                                  >  1:2557354/1‑148 (MQ=255)
                                                                                              tAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGt                                                  >  2:3091274/1‑148 (MQ=255)
                                                                                              tAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGt                                                  <  1:2889512/148‑1 (MQ=255)
                                                                                                aGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtct                                                <  1:2525840/148‑1 (MQ=255)
                                                                                                aGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtct                                                <  2:3213832/148‑1 (MQ=255)
                                                                                                aGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtct                                                >  1:1903977/1‑148 (MQ=255)
                                                                                                 gCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtctt                                               >  2:2881117/1‑148 (MQ=255)
                                                                                                 gCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtctt                                               <  1:435468/148‑1 (MQ=255)
                                                                                                   gtcgtTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtcttct                                             <  1:900036/148‑1 (MQ=255)
                                                                                                      gtTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgta                                          <  2:440298/148‑1 (MQ=255)
                                                                                                               aaCAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGa                                 <  2:2191077/148‑1 (MQ=255)
                                                                                                                aCAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGAt                                >  1:1728322/1‑148 (MQ=255)
                                                                                                                aCAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGAt                                >  1:3114138/1‑148 (MQ=255)
                                                                                                                 cAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATa                               <  2:2175788/148‑1 (MQ=255)
                                                                                                                 cAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATa                               >  1:1704571/1‑148 (MQ=255)
                                                                                                                    cGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGc                            <  1:502951/148‑1 (MQ=255)
                                                                                                                        aTCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAgg                        >  1:3503079/1‑148 (MQ=255)
                                                                                                                         tCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGc                       <  1:418705/148‑1 (MQ=255)
                                                                                                                            aaTTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGa                    <  1:418919/148‑1 (MQ=255)
                                                                                                                            aaTTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCAAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGGCTTCTGTACGTACCTGATAGGCCAGGCTGa                    >  2:3372306/1‑148 (MQ=255)
                                                                                                                               ttGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGacca                 >  1:2711558/1‑148 (MQ=255)
                                                                                                                                 gTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGaccacc               <  2:2216705/148‑1 (MQ=255)
                                                                                                                                  tCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGAccacca              <  2:1971317/148‑1 (MQ=255)
                                                                                                                                  tCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGAccacca              <  1:664910/148‑1 (MQ=255)
                                                                                                                                   cGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAg             <  1:2152190/148‑1 (MQ=255)
                                                                                                                                     ccAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGt           <  1:3272059/148‑1 (MQ=255)
                                                                                                                                        aCTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCATGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACAACCAGGTCtt        >  2:1681369/1‑148 (MQ=255)
                                                                                                                                        aCTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCtt        <  2:681634/148‑1 (MQ=255)
                                                                                                                                         cTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCttt       <  1:2463075/148‑1 (MQ=255)
                                                                                                                                         cTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCttt       >  2:2082754/1‑148 (MQ=255)
                                                                                                                                         cTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCttt       <  2:2962264/148‑1 (MQ=255)
                                                                                                                                           tCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCATGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTg     <  2:52739/148‑1 (MQ=255)
                                                                                                                                           tCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTg     <  2:2324240/148‑1 (MQ=255)
                                                                                                                                              tGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTa  >  1:520876/1‑148 (MQ=255)
                                                                                                                                              tGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTa  <  1:2298822/148‑1 (MQ=255)
                                                                                                                                              tGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTa  >  1:1050351/1‑148 (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: