Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 667,899 C→A 100% D781Y (GAT→TAT)  leuS ← leucyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273667,8990CA98.7% 817.7 / 0.5 254D781Y (GAT→TAT) leuSleucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base A (98/152);  total (100/154)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 283 total aligned reads displayed.

GCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAGTTTGTTCA  >  CP009273/667754‑668042
                                                                                                                                                 |                                                                                                                                               
gCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATAtagc                                                                                                                                               >  2:3803313/1‑148 (MQ=255)
gCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATAtagc                                                                                                                                               >  1:1263758/1‑148 (MQ=255)
gCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATAtagc                                                                                                                                               <  2:1008765/148‑1 (MQ=255)
gCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATAtagc                                                                                                                                               <  1:2133665/148‑1 (MQ=255)
 cACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATAtagcc                                                                                                                                              <  1:271563/148‑1 (MQ=255)
 cACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATAtagcc                                                                                                                                              >  2:2970133/1‑148 (MQ=255)
     ttCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATAtagccttcg                                                                                                                                          >  2:2913888/1‑148 (MQ=255)
      tCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTTTCGTTCtagccttcgc                                                                                                                                         >  1:367948/1‑148 (MQ=255)
         cgAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATagccttcgcctt                                                                                                                                      >  1:2516820/1‑148 (MQ=255)
            aCCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCa                                                                                                                                   >  1:1378973/1‑148 (MQ=255)
             ccTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAg                                                                                                                                  <  1:585149/148‑1 (MQ=255)
                gTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTc                                                                                                                               >  2:1863471/1‑148 (MQ=255)
                gTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTc                                                                                                                               <  2:3695873/148‑1 (MQ=255)
                    ttcCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGc                                                                                                                           >  2:3720952/1‑148 (MQ=255)
                    ttcCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGc                                                                                                                           >  2:3720052/1‑148 (MQ=255)
                    ttcCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGc                                                                                                                           <  1:297294/148‑1 (MQ=255)
                        gTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCcaca                                                                                                                       <  2:302272/148‑1 (MQ=255)
                              tCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGa                                                                                                                 <  2:3374142/148‑1 (MQ=255)
                                 aCCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGc                                                                                                              >  1:2966952/1‑148 (MQ=255)
                                 aCCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGc                                                                                                              >  1:2828485/1‑148 (MQ=255)
                                  ccGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTTGAGTATTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCa                                                                                                             <  1:3014932/148‑1 (MQ=255)
                                  ccGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCa                                                                                                             >  1:1693993/1‑148 (MQ=255)
                                   cGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAg                                                                                                            <  2:2219620/148‑1 (MQ=255)
                                   cGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAg                                                                                                            <  1:3447112/148‑1 (MQ=255)
                                    ggAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGa                                                                                                           <  2:1790328/148‑1 (MQ=255)
                                    ggAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGa                                                                                                           <  1:2750015/148‑1 (MQ=255)
                                     gAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGAt                                                                                                          >  2:350101/1‑148 (MQ=255)
                                       aCGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGAtgt                                                                                                        >  1:3261820/1‑148 (MQ=255)
                                       aCGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGAtgt                                                                                                        <  1:3357346/148‑1 (MQ=255)
                                       aCGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGAtgt                                                                                                        >  1:1358470/1‑148 (MQ=255)
                                        cGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGAtgtg                                                                                                       >  1:1649627/1‑148 (MQ=255)
                                         ggTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGc                                                                                                      >  1:1804910/1‑148 (MQ=255)
                                           tGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCgg                                                                                                    >  2:3400899/1‑148 (MQ=255)
                                           tGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCgg                                                                                                    <  1:3491900/148‑1 (MQ=255)
                                                ttGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGa                                                                                               <  1:3507738/148‑1 (MQ=255)
                                                ttGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGa                                                                                               <  2:1595863/148‑1 (MQ=255)
                                                ttGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGa                                                                                               <  1:1654344/148‑1 (MQ=255)
                                                   gCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACg                                                                                            >  2:278689/1‑148 (MQ=255)
                                                    cACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATATGCGGGGTGAACgg                                                                                           >  1:1821312/1‑148 (MQ=255)
                                                     aCGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACggg                                                                                          <  1:784944/148‑1 (MQ=255)
                                                      cGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGACCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGt                                                                                         >  1:224807/1‑148 (MQ=255)
                                                           tttACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGc                                                                                    <  1:903332/148‑1 (MQ=255)
                                                            ttACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCa                                                                                   <  1:166538/148‑1 (MQ=255)
                                                             tACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCAt                                                                                  <  2:3120924/148‑1 (MQ=255)
                                                             tACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCAt                                                                                  <  1:1437380/148‑1 (MQ=255)
                                                              aCCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTACGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGACCGGGTTAAGCATa                                                                                 >  2:53298/1‑148 (MQ=255)
                                                               ccGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATAc                                                                                >  1:2727612/1‑148 (MQ=255)
                                                                cGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACg                                                                               >  1:2323533/1‑148 (MQ=255)
                                                                 gTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACgg                                                                              <  2:3805069/148‑1 (MQ=255)
                                                                 gTTAACCTGCACCAAGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACgg                                                                              >  2:140754/1‑148 (MQ=255)
                                                                    aaCCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACa                                                                           >  2:661390/1‑148 (MQ=255)
                                                                    aaCCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACa                                                                           <  2:1936283/148‑1 (MQ=255)
                                                                    aaCCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACa                                                                           >  2:745851/1‑148 (MQ=255)
                                                                    aaCCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACa                                                                           <  2:3428067/148‑1 (MQ=255)
                                                                         gCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGc                                                                      <  2:1684998/148‑1 (MQ=255)
                                                                         gCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGc                                                                      <  1:1336086/148‑1 (MQ=255)
                                                                               cGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCcagcag                                                                <  2:2442361/148‑1 (MQ=255)
                                                                               cGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCcagcag                                                                >  1:519442/1‑148 (MQ=255)
                                                                               cGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCcagcag                                                                <  1:1627907/148‑1 (MQ=255)
                                                                               cGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCcagcag                                                                >  1:2710697/1‑148 (MQ=255)
                                                                                 aCCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTg                                                              <  2:932393/148‑1 (MQ=255)
                                                                                    aGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCtt                                                           <  1:2391749/148‑1 (MQ=255)
                                                                                     gCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTc                                                          <  2:132068/148‑1 (MQ=255)
                                                                                     gCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTc                                                          <  1:2847013/148‑1 (MQ=255)
                                                                                      cGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTcc                                                         <  1:2149289/148‑1 (MQ=255)
                                                                                      cGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCCTTGTCGATATAGCCATCGCCTTTCAGTTCCTTCCACAGCGTGAATCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTcc                                                         >  1:378726/1‑148 (MQ=255)
                                                                                        tGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTg                                                       >  2:1835174/1‑148 (MQ=255)
                                                                                             tCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCa                                                  >  1:876614/1‑148 (MQ=255)
                                                                                              cTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCag                                                 <  1:1531691/148‑1 (MQ=255)
                                                                                              cTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCAAGCAGTGCTTCCTGCATCag                                                 <  1:2439449/148‑1 (MQ=255)
                                                                                              cTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCAAGCAGTGCTTCCTGCATCag                                                 <  1:1696450/148‑1 (MQ=255)
                                                                                                 ccaccaTCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAgc                                              >  1:3418698/1‑148 (MQ=255)
                                                                                                  caccaTCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAgcg                                             <  1:3038967/148‑1 (MQ=255)
                                                                                                    ccaTCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGcgcg                                           >  1:2637982/1‑148 (MQ=255)
                                                                                                    ccaTCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGcgcg                                           >  2:897644/1‑148 (MQ=255)
                                                                                                       tCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATc                                        <  2:40321/148‑1 (MQ=255)
                                                                                                       tCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATc                                        <  2:2917532/148‑1 (MQ=255)
                                                                                                         gCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCt                                      >  2:3534545/1‑148 (MQ=255)
                                                                                                          cTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTg                                     <  1:2409651/148‑1 (MQ=255)
                                                                                                           tttttCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGc                                    <  2:2008761/148‑1 (MQ=255)
                                                                                                            ttttCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCt                                   <  2:1518232/148‑1 (MQ=255)
                                                                                                                   cAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATc                            <  2:2710697/148‑1 (MQ=255)
                                                                                                                   cAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCGGCATCAGAGCGCGATCCTGCTCGCCATc                            >  2:1128041/1‑148 (MQ=255)
                                                                                                                     gCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCgg                          <  1:986912/148‑1 (MQ=255)
                                                                                                                     gCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCgg                          >  1:1001372/1‑148 (MQ=255)
                                                                                                                     gCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCgg                          >  1:1001506/1‑148 (MQ=255)
                                                                                                                      cAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGt                         <  2:2089524/148‑1 (MQ=255)
                                                                                                                       aaCCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGtt                        <  2:2851498/148‑1 (MQ=255)
                                                                                                                         ccGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTgg                      <  2:2409273/148‑1 (MQ=255)
                                                                                                                           ggCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTg                    <  1:2958487/148‑1 (MQ=255)
                                                                                                                           ggCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTg                    <  1:1741917/148‑1 (MQ=255)
                                                                                                                           ggCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTg                    <  1:2563395/148‑1 (MQ=255)
                                                                                                                           ggCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTg                    <  1:1088199/148‑1 (MQ=255)
                                                                                                                             ccACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCt                  <  2:519442/148‑1 (MQ=255)
                                                                                                                                  gcgcGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGc             <  2:2618131/148‑1 (MQ=255)
                                                                                                                                    gcgTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCa           <  2:453852/148‑1 (MQ=255)
                                                                                                                                     cgTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCagag                                               <  1:480634/111‑1 (MQ=255)
                                                                                                                                        tGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAGttt       >  1:1634183/1‑148 (MQ=255)
                                                                                                                                          tCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAGTTTGt     <  2:1254557/148‑1 (MQ=255)
                                                                                                                                           cGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAGTTTGtt    <  2:3199940/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAGTTTGTTCA  >  CP009273/667754‑668042

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: