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,7050CT97.8% 604.4 / 2.8 186G51D (GGT→GAT) slyDFKBP‑type peptidyl prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (105/76);  total (107/79)
Fisher's exact test for biased strand distribution p-value = 3.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 186 total aligned reads displayed.

TTCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTG  >  CP009273/3471560‑3471846
                                                                                                                                                 |                                                                                                                                             
ttCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCt                                                                                                                                             >  1:108205/1‑148 (MQ=255)
    gCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCa                                                                                                                                         >  2:2572941/1‑148 (MQ=255)
        ggAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAgcgc                                                                                                                                     <  2:322876/148‑1 (MQ=255)
        ggAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAgcgc                                                                                                                                     <  1:1092743/148‑1 (MQ=255)
          aaaCGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACTAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCg                                                                                                                                   <  2:493908/148‑1 (MQ=255)
          aaaCGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCg                                                                                                                                   <  1:3183974/148‑1 (MQ=255)
          aaaCGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCg                                                                                                                                   <  1:585088/148‑1 (MQ=255)
          aaaCGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCg                                                                                                                                   >  2:3188190/1‑148 (MQ=255)
            aCGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGtt                                                                                                                                 >  2:1807345/1‑148 (MQ=255)
               cATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTcc                                                                                                                              >  1:2565950/1‑148 (MQ=255)
                aTACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCa                                                                                                                             >  1:393941/1‑148 (MQ=255)
                      acctGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCag                                                                                                                       >  2:2964882/1‑148 (MQ=255)
                       cctGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCaca                                                                                                                      >  1:1292390/1‑146 (MQ=255)
                         tGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAgaga                                                                                                                    >  1:3043478/1‑148 (MQ=255)
                              ttCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAgg                                                                                                               >  2:1235244/1‑148 (MQ=255)
                              ttCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAgg                                                                                                               >  1:1268811/1‑148 (MQ=255)
                              ttCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAgg                                                                                                               >  1:2491175/1‑148 (MQ=255)
                                catcaACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGa                                                                                                             <  1:663965/148‑1 (MQ=255)
                                  tcaACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAAc                                                                                                           >  1:2860939/1‑148 (MQ=255)
                                   caACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAAcc                                                                                                          >  1:1453616/1‑148 (MQ=255)
                                      cGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTg                                                                                                       <  1:881590/148‑1 (MQ=255)
                                      cGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTg                                                                                                       >  2:873219/1‑148 (MQ=255)
                                        cccATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAATTTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGAc                                                                                                     >  1:2163136/1‑148 (MQ=255)
                                         ccATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGAcc                                                                                                    >  2:3266652/1‑148 (MQ=255)
                                         ccATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGAcc                                                                                                    <  2:2806119/148‑1 (MQ=255)
                                          cATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCACATTTGGCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCa                                                                                                   <  1:1236096/148‑1 (MQ=255)
                                             aaaTACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGc                                                                                                >  1:2700172/1‑148 (MQ=255)
                                             aaaTACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGc                                                                                                >  2:3018862/1‑148 (MQ=255)
                                             aaaTACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGc                                                                                                >  2:3234785/1‑148 (MQ=255)
                                                 aCGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGt                                                                                            <  2:1743005/148‑1 (MQ=255)
                                                    tCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGt                                                                                         <  1:2509978/148‑1 (MQ=255)
                                                    tCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGt                                                                                         <  2:3403046/148‑1 (MQ=255)
                                                    tCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGt                                                                                         >  2:3646290/1‑148 (MQ=255)
                                                    tCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGt                                                                                         >  2:3646755/1‑148 (MQ=255)
                                                     cTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTc                                                                                        <  1:807495/148‑1 (MQ=255)
                                                       ttAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCa                                                                                      >  2:3213765/1‑148 (MQ=255)
                                                       ttAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCa                                                                                      >  1:1807795/1‑148 (MQ=255)
                                                       ttAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCa                                                                                      >  2:1067553/1‑148 (MQ=255)
                                                        tAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAg                                                                                     <  2:1942015/148‑1 (MQ=255)
                                                          ggAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAgcg                                                                                   >  2:1284813/1‑148 (MQ=255)
                                                          ggAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAgcg                                                                                   >  2:1181499/1‑148 (MQ=255)
                                                           gAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAgcgg                                                                                  >  2:136413/1‑148 (MQ=255)
                                                                  ttGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTc                                                                           <  1:3347691/148‑1 (MQ=255)
                                                                    gCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCAc                                                                         <  2:1380282/148‑1 (MQ=255)
                                                                       ccAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCgg                                                                      >  1:1311136/1‑148 (MQ=255)
                                                                       ccAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCgg                                                                      >  1:572666/1‑148 (MQ=255)
                                                                        cAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGga                                                                     >  1:3481342/1‑148 (MQ=255)
                                                                        cAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGga                                                                     >  1:3481357/1‑148 (MQ=255)
                                                                         aGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGgag                                                                    >  2:3484090/1‑148 (MQ=255)
                                                                          ggTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGgaga                                                                   >  2:2552002/1‑148 (MQ=255)
                                                                          ggTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGgaga                                                                   >  2:1923228/1‑148 (MQ=255)
                                                                          ggTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTACATGCCCGAGATCAGGGAACCGTGACCAAGCATGTAGTCCAGCGGCGCACTCACCCgaga                                                                   >  1:3524427/1‑148 (MQ=255)
                                                                              ttCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACtca                                                               <  2:1077482/148‑1 (MQ=255)
                                                                               tcgtcgTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACtcat                                                              <  1:795324/148‑1 (MQ=255)
                                                                                 gtcgtACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACtcatca                                                            <  1:450986/148‑1 (MQ=255)
                                                                                  tcgtACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATcaa                                                           >  1:2971646/1‑148 (MQ=255)
                                                                                   cgtACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATcaac                                                          >  2:554337/1‑148 (MQ=255)
                                                                                          aCCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAacaca                                                   >  2:2470572/1‑148 (MQ=255)
                                                                                           ccGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAacacac                                                  <  1:1078040/148‑1 (MQ=255)
                                                                                           ccGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAacacac                                                  >  1:1514094/1‑148 (MQ=255)
                                                                                              tAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGt                                               <  1:1877951/148‑1 (MQ=255)
                                                                                              tAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGt                                               >  2:2095354/1‑148 (MQ=255)
                                                                                                aGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtct                                             <  1:1154934/148‑1 (MQ=255)
                                                                                                 gCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCATCGCCGTTTCCAGGCCAGAGATAAGGGAACCGTGACCATCCACGTAGTCCCGCGGCGCACTCACCGGAGACTCATCAACAAACACAACGtcat                                            >  2:2766524/1‑146 (MQ=255)
                                                                                                 gCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtctt                                            >  2:2371076/1‑148 (MQ=255)
                                                                                                  cgtcgtTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTACAGGCCCCAGATCACGGAACCGTGACCAAGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCCtcttc                                           >  1:3388485/1‑148 (MQ=255)
                                                                                                    tcgtTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGAACATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTg                                         >  1:1820980/1‑148 (MQ=255)
                                                                                                      gtTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgta                                       <  2:1655788/148‑1 (MQ=255)
                                                                                                       ttCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgtac                                      >  2:1011577/1‑148 (MQ=255)
                                                                                                        tCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgtacg                                     <  2:350810/148‑1 (MQ=255)
                                                                                                        tCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgtacg                                     <  2:1849401/148‑1 (MQ=255)
                                                                                                        tCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCAAGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgtacg                                     >  2:129842/1‑148 (MQ=255)
                                                                                                         cgcgCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgtacgt                                    <  2:2171021/148‑1 (MQ=255)
                                                                                                           cgcCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgtacgtac                                  <  1:678599/148‑1 (MQ=255)
                                                                                                               aaCAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCATGCCAGAGATCAGGGAAACGTGACCAAGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATAAACCAACACACCGTCTTGTGTAAGTACCTGa                              >  1:2892859/1‑148 (MQ=255)
                                                                                                               aaCAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGa                              >  1:832062/1‑148 (MQ=255)
                                                                                                               aaCAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCGTCTGTACGTACCTGa                              >  1:580004/1‑148 (MQ=255)
                                                                                                                aCAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCCACACACCGTCTTCTGTACGTACCTGAt                             >  1:733059/1‑148 (MQ=255)
                                                                                                                aCAGCGACATCAAATTTGACGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGAAAAGGGAAACGTGACAATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCAACAACCAACACACCCTCTACAGTACGTAACAGAt                             >  2:1109257/1‑148 (MQ=255)
                                                                                                                 cAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATa                            <  2:115333/148‑1 (MQ=255)
                                                                                                                 cAGCGACATCAAATTTGTCGCCAACTTCATGATCATCCAGCGCCGTTTCCACCCAAGATATCACCGAAACCTCACCCTGCAGGTAGTCCAGCGCCGCACTCACCCCACACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATa                            >  1:1411460/1‑148 (MQ=255)
                                                                                                                      aCATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCa                       >  1:1802685/1‑148 (MQ=255)
                                                                                                                       cATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAg                      >  1:2197566/1‑148 (MQ=255)
                                                                                                                        aTCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAgg                     >  2:1530081/1‑148 (MQ=255)
                                                                                                                        aTCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAgg                     >  2:3117948/1‑148 (MQ=255)
                                                                                                                         tCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGc                    <  1:3224855/148‑1 (MQ=255)
                                                                                                                         tCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTGTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCCGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGc                    >  1:649739/1‑148 (MQ=255)
                                                                                                                          cAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCCTGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCt                   >  2:487681/1‑148 (MQ=255)
                                                                                                                               ttGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGacca              <  1:1633750/148‑1 (MQ=255)
                                                                                                                                     ccAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGGCGTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGt        >  1:2022040/1‑148 (MQ=255)
                                                                                                                                       aaCTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCt      <  2:118447/148‑1 (MQ=255)
                                                                                                                                        aCTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCCGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCCGGTGtt     >  2:660250/1‑148 (MQ=255)
                                                                                                                                        aCTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCtt     <  2:3378590/148‑1 (MQ=255)
                                                                                                                                        aCTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCtt     >  1:3690822/1‑148 (MQ=255)
                                                                                                                                         cTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCttt    <  1:1181499/148‑1 (MQ=255)
                                                                                                                                         cTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCttt    <  2:2664860/148‑1 (MQ=255)
                                                                                                                                         cTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCttt    >  1:1920474/1‑148 (MQ=255)
                                                                                                                                         cTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCttt    >  1:3471522/1‑148 (MQ=255)
                                                                                                                                         cTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTACAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCttt    <  1:3086386/148‑1 (MQ=255)
                                                                                                                                           tCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTg  >  2:950921/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
TTCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTG  >  CP009273/3471560‑3471846

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