Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,715,643 G→T 100% E778* (GAA→TAA)  pka → protein lysine acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,715,6430GT98.1% 968.5 / 7.4 285E778* (GAA→TAA) pkaprotein lysine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/5);  new base T (128/150);  total (130/155)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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 286 total aligned reads displayed.

CGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCT  >  CP009273/2715496‑2715790
                                                                                                                                                   |                                                                                                                                                   
 gCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTaaa                                                                                                                                                  <  2:2796758/150‑1 (MQ=255)
      aTTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCg                                                                                                                                             >  1:3400210/1‑150 (MQ=255)
      aTTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCg                                                                                                                                             >  1:4095658/1‑150 (MQ=255)
        tttGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGtt                                                                                                                                           <  2:1274885/150‑1 (MQ=255)
         ttGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGttt                                                                                                                                          <  2:352251/150‑1 (MQ=255)
          tGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTg                                                                                                                                         >  2:568648/1‑150 (MQ=255)
           gCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGt                                                                                                                                        >  1:420541/1‑150 (MQ=255)
              tgaagtTGAGCCACAACTTCATCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGc                                                                                                                                     <  1:315429/149‑1 (MQ=255)
              agaagaTGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGc                                                                                                                                     <  1:4615677/150‑1 (MQ=255)
               gaagaTGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCg                                                                                                                                    <  1:1534946/150‑1 (MQ=255)
                 agaTGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGt                                                                                                                                  >  2:3153249/1‑150 (MQ=255)
                 agaTGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGt                                                                                                                                  >  2:3328188/1‑150 (MQ=255)
                   aTGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTac                                                                                                                                <  1:4501611/150‑1 (MQ=255)
                   aTGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTac                                                                                                                                <  2:678658/150‑1 (MQ=255)
                    tGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTacg                                                                                                                               <  2:2888236/150‑1 (MQ=255)
                      aGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTacgac                                                                                                                             <  1:1056095/150‑1 (MQ=255)
                       gCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTacgacg                                                                                                                            >  1:4138798/1‑150 (MQ=255)
                        ccACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTCCGATCGGTAAATGGCGTTTGTAGCGGTACGACGt                                                                                                                           >  2:1001675/1‑150 (MQ=255)
                        ccACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGt                                                                                                                           >  1:1791582/1‑150 (MQ=255)
                                    cAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACg                                                                                                               <  2:4597972/150‑1 (MQ=255)
                                     aaTTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACgg                                                                                                              >  2:3242317/1‑150 (MQ=255)
                                       ttCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGaa                                                                                                            <  1:238855/150‑1 (MQ=255)
                                         cATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAaga                                                                                                          >  2:2692833/1‑150 (MQ=255)
                                          aTTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAagag                                                                                                         <  2:18941/150‑1 (MQ=255)
                                            ttCGCGAGTCACCAAAGAAGATCTTTATTACCTCTACTTTAGCGAGATCAACTAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGAt                                                                                                       <  1:2778340/150‑1 (MQ=255)
                                               gcgAGTCACCAAATAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGATTTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCt                                                                                                    <  2:2204611/150‑1 (MQ=255)
                                               gcgAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCt                                                                                                    >  2:1708532/1‑150 (MQ=255)
                                               gcgAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCt                                                                                                    >  1:3275795/1‑150 (MQ=255)
                                               gcgAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCt                                                                                                    >  2:4133927/1‑150 (MQ=255)
                                                 gAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCg                                                                                                  >  2:267163/1‑150 (MQ=255)
                                                  aGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCgg                                                                                                 >  1:1963719/1‑150 (MQ=255)
                                                   gTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGc                                                                                                <  1:3908625/150‑1 (MQ=255)
                                                     cACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGt                                                                                              <  1:1853247/150‑1 (MQ=255)
                                                       ccAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCa                                                                                            >  2:4540774/1‑150 (MQ=255)
                                                       ccAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCa                                                                                            <  2:4174068/150‑1 (MQ=255)
                                                         aaaGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCAcg                                                                                          >  1:1136103/1‑150 (MQ=255)
                                                          aagaagATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCAcgc                                                                                         >  1:1365292/1‑150 (MQ=255)
                                                           agaagaTCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCAcgcg                                                                                        >  1:978532/1‑150 (MQ=255)
                                                           agaagaTCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCAcgcg                                                                                        <  1:3621123/150‑1 (MQ=255)
                                                           agaagaTCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCAcgcg                                                                                        <  1:3621285/150‑1 (MQ=255)
                                                             aagaTCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTg                                                                                      >  1:4013429/1‑150 (MQ=255)
                                                              agaTCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGc                                                                                     >  2:2241292/1‑150 (MQ=255)
                                                              agaTCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGc                                                                                     <  1:857510/150‑1 (MQ=255)
                                                                  cTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGAtt                                                                                 <  2:475040/150‑1 (MQ=255)
                                                                   tttATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGAttt                                                                                <  1:2229683/150‑1 (MQ=255)
                                                                   ttgattaCCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGAttt                                                                                <  2:1100279/147‑1 (MQ=255)
                                                                    ttattaCCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTc                                                                               <  1:4353872/150‑1 (MQ=255)
                                                                    ttattaCCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTc                                                                               >  1:852324/1‑150 (MQ=255)
                                                                     tattaCCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTcc                                                                              <  1:247125/150‑1 (MQ=255)
                                                                        taCCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGAt                                                                           <  1:2526140/150‑1 (MQ=255)
                                                                        taCCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGAt                                                                           <  1:2526162/150‑1 (MQ=255)
                                                                           cGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCt                                                                        <  1:359646/150‑1 (MQ=255)
                                                                            gCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTg                                                                       <  1:1065661/150‑1 (MQ=255)
                                                                            gCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTg                                                                       <  1:1279346/150‑1 (MQ=255)
                                                                               aCTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATa                                                                    <  2:4585769/150‑1 (MQ=255)
                                                                                    aGCGAGATCAACGAATTTACGCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATc                                                               <  1:3636627/150‑1 (MQ=255)
                                                                                    aGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATc                                                               >  2:3774566/1‑150 (MQ=255)
                                                                                        agaTCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATg                                                           <  2:3991200/150‑1 (MQ=255)
                                                                                        agaTCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATg                                                           >  2:3822333/1‑150 (MQ=255)
                                                                                           tCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCg                                                        <  2:890791/150‑1 (MQ=255)
                                                                                            cAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGa                                                       >  2:2939522/1‑150 (MQ=255)
                                                                                               cGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGaagt                                                    >  2:1227059/1‑148 (MQ=255)
                                                                                               cGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGAGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAAtt                                                    <  1:181487/150‑1 (MQ=255)
                                                                                                gAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAAttt                                                   <  1:2054716/150‑1 (MQ=255)
                                                                                                gAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAAttt                                                   >  1:2616969/1‑150 (MQ=255)
                                                                                                 aaTTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTg                                                  <  2:1401989/150‑1 (MQ=255)
                                                                                                 aaTTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTg                                                  <  2:1142114/150‑1 (MQ=255)
                                                                                                 aaTTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTg                                                  >  2:529918/1‑150 (MQ=255)
                                                                                                  aTTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTATCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGc                                                 <  1:1214941/150‑1 (MQ=255)
                                                                                                  aTTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGc                                                 >  1:3822533/1‑150 (MQ=255)
                                                                                                        ccATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACt                                           <  2:3837303/150‑1 (MQ=255)
                                                                                                        ccATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACt                                           <  2:3837248/150‑1 (MQ=255)
                                                                                                          aTGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAAGGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTgg                                         >  2:1186601/1‑150 (MQ=255)
                                                                                                          aTGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTgg                                         >  2:60425/1‑150 (MQ=255)
                                                                                                              aGATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCg                                     >  1:683014/1‑150 (MQ=255)
                                                                                                               gATTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGc                                    >  1:289440/1‑150 (MQ=255)
                                                                                                                aTTTAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCt                                   <  1:4604622/150‑1 (MQ=255)
                                                                                                                   tAGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCgg                                <  1:372242/150‑1 (MQ=255)
                                                                                                                    aGCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGa                               >  1:1990019/1‑150 (MQ=255)
                                                                                                                     gCCAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGAt                              >  2:4479860/1‑150 (MQ=255)
                                                                                                                       cAACATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGAtct                            <  1:2312827/150‑1 (MQ=255)
                                                                                                                        aaCATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGAtctc                           <  1:3224443/150‑1 (MQ=255)
                                                                                                                        aaCATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGAtctc                           >  1:476139/1‑150 (MQ=255)
                                                                                                                        aaCATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGAtctc                           <  2:2394357/150‑1 (MQ=255)
                                                                                                                        aaCATGACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGAtctc                           <  2:4212596/150‑1 (MQ=255)
                                                                                                                             gACACAGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAgg                      >  1:2584428/1‑150 (MQ=255)
                                                                                                                                acaGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGtt                   >  1:3551097/1‑150 (MQ=255)
                                                                                                                                acaGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGtt                   <  1:3350533/150‑1 (MQ=255)
                                                                                                                                  aGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAg                 <  1:1101135/150‑1 (MQ=255)
                                                                                                                                  aGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAg                 >  2:1824886/1‑150 (MQ=255)
                                                                                                                                  aGATCGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAg                 >  1:4062062/1‑150 (MQ=255)
                                                                                                                                      cGACTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCtt             >  1:148246/1‑150 (MQ=255)
                                                                                                                                         cTACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAgg          >  1:2697002/1‑150 (MQ=255)
                                                                                                                                          tACGATCGGTAACTGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGt         <  2:1996245/150‑1 (MQ=255)
                                                                                                                                          tACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGt         <  1:126596/150‑1 (MQ=255)
                                                                                                                                          tACGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGt         <  2:3764953/150‑1 (MQ=255)
                                                                                                                                           aCGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTc        >  1:2081344/1‑150 (MQ=255)
                                                                                                                                            cGATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCg       <  2:1521909/150‑1 (MQ=255)
                                                                                                                                             gATCGGTAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGa      >  2:2324874/1‑150 (MQ=255)
                                                                                                                                                 gggAAATGGCGTTTGTAGCGGTCCGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCt  >  2:3089675/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGCCCGATTTTGCCAGAAGATGAGCCACAACTTCAGCAATTCATTTCGCGAGTCACCAAAGAAGATCTTTATTACCGCTACTTTAGCGAGATCAACGAATTTACCCATGAAGATTTAGCCAACATGACACAGATCGACTACGATCGGGAAATGGCGTTTGTAGCGGTACGACGTATTGATCAAACGGAAGAGATCCTCGGCGTCACGCGTGCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCT  >  CP009273/2715496‑2715790

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: