Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,471,705 C→T 94.5% 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,7050CT94.5% 523.1 / 17.0 164G51D (GGT→GAT) slyDFKBP‑type peptidyl prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base C (6/3);  new base T (78/77);  total (84/80)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 165 total aligned reads displayed.

GTTTCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTA  >  CP009273/3471558‑3471849
                                                                                                                                                   |                                                                                                                                                
gTTTCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGAc                                                                                                                                                  <  2:49888/148‑1 (MQ=255)
  ttCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCt                                                                                                                                                >  1:3229533/1‑148 (MQ=255)
  ttCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCt                                                                                                                                                >  1:1520756/1‑148 (MQ=255)
   tCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCtt                                                                                                                                               >  2:1531985/1‑148 (MQ=255)
     aGCCAGGAAACGCATACCTACCTGCAGGTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTcc                                                                                                                                             <  2:3077221/148‑1 (MQ=255)
       ccAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAg                                                                                                                                           <  1:1979595/148‑1 (MQ=255)
          ggAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAgcgc                                                                                                                                        >  1:1023331/1‑148 (MQ=255)
           gAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGcc                                                                                                                                       >  1:3011509/1‑148 (MQ=255)
            aaaCGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCg                                                                                                                                      >  1:2300452/1‑148 (MQ=255)
             aaCGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGt                                                                                                                                     <  2:1548935/148‑1 (MQ=255)
             aaCGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGt                                                                                                                                     <  2:2844479/148‑1 (MQ=255)
             aaCGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGt                                                                                                                                     <  2:3189414/148‑1 (MQ=255)
              aCGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGtt                                                                                                                                    <  1:3253212/148‑1 (MQ=255)
              aCGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGtt                                                                                                                                    <  1:123906/148‑1 (MQ=255)
               cGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGttt                                                                                                                                   <  2:2568441/148‑1 (MQ=255)
               cGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGttt                                                                                                                                   <  2:1487678/148‑1 (MQ=255)
                gCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTc                                                                                                                                  <  1:2278126/148‑1 (MQ=255)
                gCATACCTACCTGCAGGTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTc                                                                                                                                  <  2:2536367/148‑1 (MQ=255)
                  aTACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCa                                                                                                                                <  2:1941974/148‑1 (MQ=255)
                   tacctaccTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCATACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAg                                                                                                                               >  1:2467017/1‑148 (MQ=255)
                    acctacctGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAgg                                                                                                                              <  2:103457/148‑1 (MQ=255)
                     cctacctGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGc                                                                                                                             >  1:1855840/1‑148 (MQ=255)
                      ctacctGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGcc                                                                                                                            <  1:1683648/148‑1 (MQ=255)
                       tacctGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCa                                                                                                                           >  2:602830/1‑148 (MQ=255)
                        acctGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCag                                                                                                                          >  2:1860577/1‑148 (MQ=255)
                        acctGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCag                                                                                                                          <  2:1670773/148‑1 (MQ=255)
                        acctGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTACAGGCCag                                                                                                                          >  1:450054/1‑148 (MQ=255)
                           tGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAgaga                                                                                                                       <  2:1456367/148‑1 (MQ=255)
                             cAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATc                                                                                                                     >  2:1261639/1‑148 (MQ=255)
                               gTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAg                                                                                                                   >  1:1141912/1‑148 (MQ=255)
                                  catcaACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGa                                                                                                                <  2:3270554/148‑1 (MQ=255)
                                    tcaACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAAc                                                                                                              >  1:2874216/1‑148 (MQ=255)
                                     caACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGAAATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAAcc                                                                                                             <  1:897939/148‑1 (MQ=255)
                                      aaCGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCg                                                                                                            >  1:2298271/1‑148 (MQ=255)
                                            cATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCa                                                                                                      <  2:1343767/148‑1 (MQ=255)
                                              tAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATg                                                                                                    >  1:899517/1‑148 (MQ=255)
                                                   aCGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGt                                                                                               >  2:792808/1‑148 (MQ=255)
                                                   aCGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGt                                                                                               >  2:2280958/1‑148 (MQ=255)
                                                      tCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGt                                                                                            >  2:2852489/1‑148 (MQ=255)
                                                       cTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTc                                                                                           <  1:2234685/148‑1 (MQ=255)
                                                       cTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTc                                                                                           <  1:2134964/148‑1 (MQ=255)
                                                         ttAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCa                                                                                         >  1:1565697/1‑148 (MQ=255)
                                                         ttAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCa                                                                                         >  2:409623/1‑148 (MQ=255)
                                                             gAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAgg                                                                                                <  2:2585070/137‑1 (MQ=255)
                                                               acacGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAgcggcg                                                                                   >  2:1281979/1‑148 (MQ=255)
                                                                   gTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACt                                                                               <  2:3291948/148‑1 (MQ=255)
                                                                     tGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCa                                                                             <  1:1206888/148‑1 (MQ=255)
                                                                      gCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCAc                                                                            <  1:2783096/148‑1 (MQ=255)
                                                                       cACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCAcc                                                                           >  1:2591632/1‑148 (MQ=255)
                                                                         ccAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCgg                                                                         >  1:1762942/1‑148 (MQ=255)
                                                                         ccAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCgg                                                                         <  1:3054402/148‑1 (MQ=255)
                                                                         ccAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCgg                                                                         >  2:912810/1‑148 (MQ=255)
                                                                          cAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGga                                                                        >  1:2476865/1‑148 (MQ=255)
                                                                           aGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTACAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCATCGGCGCACTCACCGgag                                                                       >  1:2851892/1‑148 (MQ=255)
                                                                               tttCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACtc                                                                   <  2:2824142/148‑1 (MQ=255)
                                                                                ttCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACtca                                                                  >  2:2405616/1‑148 (MQ=255)
                                                                                 tcgtcgTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACtcat                                                                 <  2:951791/148‑1 (MQ=255)
                                                                                   gtcgtACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACtcatca                                                               <  2:1141912/148‑1 (MQ=255)
                                                                                   gtcgtACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACtcatca                                                               <  1:3286934/148‑1 (MQ=255)
                                                                                   gtcgtACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACtcatca                                                               <  1:796595/148‑1 (MQ=255)
                                                                                       tACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATcaacca                                                           >  1:1153738/1‑148 (MQ=255)
                                                                                       tACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATcaacca                                                           >  1:1153868/1‑148 (MQ=255)
                                                                                        aCTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATcaaccaa                                                          >  2:1845557/1‑148 (MQ=255)
                                                                                         cTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATcaaccaac                                                         <  1:1468087/148‑1 (MQ=255)
                                                                                                 aaGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtc                                                 >  1:2896383/1‑148 (MQ=255)
                                                                                                  aGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCACCACACCGtct                                                <  2:326156/148‑1 (MQ=255)
                                                                                                   gCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtctt                                               <  2:1023331/148‑1 (MQ=255)
                                                                                                        gtTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgta                                          <  2:2407483/148‑1 (MQ=255)
                                                                                                           cgcgCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgtacgt                                       <  2:2874048/148‑1 (MQ=255)
                                                                                                             cgcCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgtacgtac                                     <  2:1832821/148‑1 (MQ=255)
                                                                                                             cgcCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgtacgtac                                     <  2:158665/148‑1 (MQ=255)
                                                                                                             cgcCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgtacgtac                                     <  1:1845267/148‑1 (MQ=255)
                                                                                                                cAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTg                                  <  2:1000166/148‑1 (MQ=255)
                                                                                                                cAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTg                                  <  1:3075365/148‑1 (MQ=255)
                                                                                                                  aCAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGAt                                <  2:1762942/148‑1 (MQ=255)
                                                                                                                   cAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATa                               >  2:1828221/1‑148 (MQ=255)
                                                                                                                   cAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATa                               >  2:178666/1‑148 (MQ=255)
                                                                                                                   cAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATa                               >  2:178669/1‑148 (MQ=255)
                                                                                                                      cGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGc                            <  2:968537/148‑1 (MQ=255)
                                                                                                                      cGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGc                            >  2:884093/1‑148 (MQ=255)
                                                                                                                      cGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGc                            <  2:671420/148‑1 (MQ=255)
                                                                                                                       gACATCAACTTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGcc                           <  2:525272/148‑1 (MQ=255)
                                                                                                                        aCATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCa                          <  1:1531985/148‑1 (MQ=255)
                                                                                                                        aCATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCa                          <  1:1531947/148‑1 (MQ=255)
                                                                                                                          aTCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAgg                        <  1:2257795/148‑1 (MQ=255)
                                                                                                                          aTCAAATTTGGCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCTCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAgg                        <  2:1728887/148‑1 (MQ=255)
                                                                                                                            cAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCt                      >  2:1794026/1‑148 (MQ=255)
                                                                                                                            cAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCt                      >  1:161093/1‑148 (MQ=255)
                                                                                                                            cAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACACCAGACACATCAACAAACACACCGTCATCTCTACGTACCTGATAGCCCAggat                      >  2:2440957/1‑146 (MQ=255)
                                                                                                                              aaTTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGa                    <  1:665087/148‑1 (MQ=255)
                                                                                                                               aTTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGac                   >  2:1884189/1‑148 (MQ=255)
                                                                                                                                       ccAACTTCCTGATCGTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGt           <  1:531345/148‑1 (MQ=255)
                                                                                                                                       ccAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGt           >  2:2277062/1‑148 (MQ=255)
                                                                                                                                         aaCTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATAAGGGAACCGTGAACATTCAGGTAGTCCAGCGGCGCACTCACCCGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATATGCCAGGCTGACCACCACGTCt         >  1:6402/1‑148 (MQ=255)
                                                                                                                                          aCTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCtt        >  1:1813231/1‑148 (MQ=255)
                                                                                                                                          aCTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACAACGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCtt        <  1:602830/148‑1 (MQ=255)
                                                                                                                                             tCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTg     <  2:2300452/148‑1 (MQ=255)
                                                                                                                                             tCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTg     <  1:2800649/148‑1 (MQ=255)
                                                                                                                                              cATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCAGATAGGCCAGGCTGACCACCAGGTCtttttc    >  1:1975339/1‑146 (MQ=255)
                                                                                                                                                tGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTa  >  2:3211221/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GTTTCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTA  >  CP009273/3471558‑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: