Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,770,281 T→A 62.8% D146E (GAT→GAA ydiF → putative acetyl‑CoA:acetoacetyl‑CoA transferase: alpha subunit/beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,770,2810TA62.8% 230.5 / 282.9 218D146E (GAT→GAAydiFputative acetyl‑CoA:acetoacetyl‑CoA transferase: alpha subunit/beta subunit
Reads supporting (aligned to +/- strand):  ref base T (42/39);  new base A (59/78);  total (101/117)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 220 total aligned reads displayed.

GGACAATCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGT  >  CP009273/1770132‑1770422
                                                                                                                                                     |                                                                                                                                             
   cAATCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCa                                                                                                                                            >  1:574534/1‑150 (MQ=255)
   cAATCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACAGTTGTCGATCCa                                                                                                                                            >  2:4086189/1‑150 (MQ=255)
    aaTCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCAc                                                                                                                                           >  1:4234627/1‑150 (MQ=255)
     aTCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACg                                                                                                                                          <  2:3927200/150‑1 (MQ=255)
     aTCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACg                                                                                                                                          <  1:4776399/150‑1 (MQ=255)
     aTCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACg                                                                                                                                          <  2:3571646/150‑1 (MQ=255)
       cgccgcGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGcc                                                                                                                                        >  2:3719287/1‑150 (MQ=255)
         ccgcGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAg                                                                                                                                      >  2:3096185/1‑150 (MQ=255)
             gTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAg                                                                                                                                  >  2:4197557/1‑150 (MQ=255)
             gTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACAAAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAg                                                                                                                                  <  1:435107/150‑1 (MQ=255)
                    tGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCACGGCGGCaa                                                                                                                           >  2:1244949/1‑150 (MQ=255)
                    tGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCaa                                                                                                                           >  2:1838222/1‑150 (MQ=255)
                     gAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCaaa                                                                                                                          >  1:1132717/1‑150 (MQ=255)
                     gAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCaaa                                                                                                                          <  2:4287140/150‑1 (MQ=255)
                       acccGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACt                                                                                                                        <  1:858355/147‑1 (MQ=255)
                        cTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACtg                                                                                                                       >  2:1621948/1‑150 (MQ=255)
                           gCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACtgaat                                                                                                                    <  2:3438170/150‑1 (MQ=255)
                              gAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACtgaatgaa                                                                                                                 >  1:2994825/1‑150 (MQ=255)
                              gAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACtgaatgaa                                                                                                                 >  1:2994874/1‑150 (MQ=255)
                                aCAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGt                                                                                                               <  1:3959389/150‑1 (MQ=255)
                                  aaaataaaatTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCa                                                                                                             <  2:3343265/150‑1 (MQ=255)
                                    aataaaatTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACt                                                                                                           >  2:370640/1‑150 (MQ=255)
                                    aataaaatTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACt                                                                                                           >  2:1712176/1‑150 (MQ=255)
                                    aataaaatTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACt                                                                                                           <  1:4739358/150‑1 (MQ=255)
                                      taaaatTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTaa                                                                                                         <  2:4251003/150‑1 (MQ=255)
                                       aaaatTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTaaa                                                                                                        <  2:2188912/150‑1 (MQ=255)
                                       aaaatTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTaaa                                                                                                        <  1:4213468/150‑1 (MQ=255)
                                        aaatTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAaag                                                                                                       >  2:927972/1‑150 (MQ=255)
                                         aatTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAaaga                                                                                                      <  1:2921385/150‑1 (MQ=255)
                                         aatTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAaaga                                                                                                      >  1:4448707/1‑150 (MQ=255)
                                           ttattGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAaagaag                                                                                                    <  2:2621968/150‑1 (MQ=255)
                                           ttattGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAaagaag                                                                                                    <  2:4783723/150‑1 (MQ=255)
                                             attGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGAc                                                                                                  <  1:2703725/150‑1 (MQ=255)
                                               tGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCt                                                                                                <  2:3111374/150‑1 (MQ=255)
                                                 cTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGa                                                                                              >  1:127720/1‑150 (MQ=255)
                                                    ataACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTa                                                                                           >  2:1196835/1‑150 (MQ=255)
                                                     taACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTaa                                                                                          >  2:4482760/1‑150 (MQ=255)
                                                     taACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGGATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTaa                                                                                          >  2:2473075/1‑150 (MQ=255)
                                                      aaCTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTaaa                                                                                         >  2:2592730/1‑150 (MQ=255)
                                                       aCTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAAc                                                                                        >  2:4567748/1‑150 (MQ=255)
                                                       aCTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAAc                                                                                        >  2:2725163/1‑150 (MQ=255)
                                                           cccACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGt                                                                                    >  1:4514070/1‑150 (MQ=255)
                                                           cccACAAGGTGTACTGACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGt                                                                                    <  1:3877145/150‑1 (MQ=255)
                                                                aaGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGt                                                                               >  1:2122710/1‑150 (MQ=255)
                                                                   gtgtACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTg                                                                            <  1:4180400/150‑1 (MQ=255)
                                                                   gtgtACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTg                                                                            >  1:727910/1‑150 (MQ=255)
                                                                   gtgtACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTg                                                                            >  2:382531/1‑150 (MQ=255)
                                                                    tgtACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGa                                                                           >  2:1374119/1‑150 (MQ=255)
                                                                     gtACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGAt                                                                          <  2:2051074/150‑1 (MQ=255)
                                                                       aCTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATaa                                                                        <  1:1422681/150‑1 (MQ=255)
                                                                        cTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAAc                                                                       <  1:2211002/150‑1 (MQ=255)
                                                                         ttACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACa                                                                      <  1:1231784/150‑1 (MQ=255)
                                                                         ttACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACa                                                                      <  2:3312314/150‑1 (MQ=255)
                                                                           acacAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACaaa                                                                    <  2:4136805/150‑1 (MQ=255)
                                                                           acacAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACaaa                                                                    >  1:2556208/1‑150 (MQ=255)
                                                                            cacaAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAg                                                                   >  1:187409/1‑150 (MQ=255)
                                                                             acaAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGa                                                                  >  2:2547061/1‑150 (MQ=255)
                                                                                aaCCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAata                                                               <  2:1449381/150‑1 (MQ=255)
                                                                                aaCCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAata                                                               >  1:3429205/1‑150 (MQ=255)
                                                                                 aCCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAatat                                                              >  1:823330/1‑150 (MQ=255)
                                                                                 aCCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAatat                                                              <  2:1587841/150‑1 (MQ=255)
                                                                                 aCCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAatat                                                              <  1:2748048/150‑1 (MQ=255)
                                                                                   cTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATAtct                                                            >  2:2538281/1‑150 (MQ=255)
                                                                                   cTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATAtct                                                            <  2:1212830/150‑1 (MQ=255)
                                                                                                   gCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGAtt                                            >  1:7975/1‑150 (MQ=255)
                                                                                                   gCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGAtt                                            <  2:4192886/150‑1 (MQ=255)
                                                                                                   gCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGAtt                                            >  1:259752/1‑150 (MQ=255)
                                                                                                    cccACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTg                                           >  2:3364096/1‑150 (MQ=255)
                                                                                                      caccaGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTgcg                                         <  2:1021400/150‑1 (MQ=255)
                                                                                                          aGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAg                                     <  1:1293722/150‑1 (MQ=255)
                                                                                                          aGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAg                                     >  1:42938/1‑150 (MQ=255)
                                                                                                          aGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAg                                     >  2:4207631/1‑150 (MQ=255)
                                                                                                          aGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAg                                     >  2:4207695/1‑150 (MQ=255)
                                                                                                             cTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGata                                  <  2:4831667/150‑1 (MQ=255)
                                                                                                             cTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGata                                  <  2:2512499/150‑1 (MQ=255)
                                                                                                             cTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGata                                  <  1:1798892/150‑1 (MQ=255)
                                                                                                             cTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGata                                  <  1:3100785/150‑1 (MQ=255)
                                                                                                              tGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGatat                                 >  1:504730/1‑150 (MQ=255)
                                                                                                              tGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGatat                                 >  1:3727424/1‑150 (MQ=255)
                                                                                                               ggTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATAtt                                <  2:309239/150‑1 (MQ=255)
                                                                                                               ggTATTATTAGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATAtt                                >  1:2316422/1‑150 (MQ=255)
                                                                                                                 tattatTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGc                              >  2:2773896/1‑150 (MQ=255)
                                                                                                                    tattaGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCtt                           >  1:3182369/1‑150 (MQ=255)
                                                                                                                      ttaGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCa                         <  2:1469836/150‑1 (MQ=255)
                                                                                                                       taGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCAt                        >  1:1659176/1‑150 (MQ=255)
                                                                                                                       taGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCAt                        <  2:2258475/150‑1 (MQ=255)
                                                                                                                       taGTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCAt                        <  2:2801407/150‑1 (MQ=255)
                                                                                                                         gTGATATTGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTc                      <  2:3196594/150‑1 (MQ=255)
                                                                                                                               ttGGCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTa                <  2:574534/150‑1 (MQ=255)
                                                                                                                                tGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTac               <  1:1546773/150‑1 (MQ=255)
                                                                                                                                 ggCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTacc              >  2:4635803/1‑150 (MQ=255)
                                                                                                                                 ggCATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTacc              <  1:1196835/150‑1 (MQ=255)
                                                                                                                                   cATCGGGACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTaccac            >  1:3809775/1‑150 (MQ=255)
                                                                                                                                    aTCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTaccacc           >  2:353130/1‑150 (MQ=255)
                                                                                                                                       gggACATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGc        >  2:577206/1‑150 (MQ=255)
                                                                                                                                          aCATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGAc     >  1:3716012/1‑150 (MQ=255)
                                                                                                                                          aCATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGAc     <  1:4567016/150‑1 (MQ=255)
                                                                                                                                          aCATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGAc     >  2:2363713/1‑150 (MQ=255)
                                                                                                                                           cATTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACa    <  1:3096185/150‑1 (MQ=255)
                                                                                                                                            aTTTGTCGAACCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAg   >  1:3216616/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
GGACAATCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGT  >  CP009273/1770132‑1770422

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: